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ABSTRACT 

A  procedural  approach  to  the  design  and  Implementation 
of  a  computer  simulation  model  is  presented,  with  a  simula- 
tion model  description  and  computer  code.   The  discussion 
and  example  are  intended  as  reference  material  for  the 
computer  science  and  computer  war  gsunlng  courses  offered 
at  the  United  States  Naval  Postgraduate  School.   The 
sample  computer  simulation  was  designed  for  the  statisti- 
cal analysis  of  the  comparative  effectiveness  of  different 
ASV/  helicopter  search  tactics  in  a  variety  of  tactical  and 
physical  environments.   This  simulation  has  available  a 
wide  range  of  input  parameters  and  is  applicable  to  all 
AS'fl   helicopters  and  any  search  plan  employing  ten  helicop- 
ters or  less.   The  accompanying  FORTRAN  computer  code  was 
written  for  the  CDC  1604  computer,  and  is  adaptable  to  the 
IBM  7090-94  by  the  inclusion  of  the  appropriate  control 
cards  and  random  number  generator. 
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CHAPTER  I 
INTRODUCTION 

1.1     Ob.lective. 

The  objective  of  this  thesis  is  to  offer  the  beginning 
student  of  computer  war  game  simulation  an  easily  understood 
example  of  what  a  computer  simulation  is  and  how  one  might 
go  about  constructing  one.   The  first  year  Operations 
Analysis  students  of  the  U.  S«  Naval  Postgraduate  School 
constitute  the  particular  audience  for  which  the  thesis  is 
intended.   Accordingly j,  no  more  than  a  basic  understanding 
of  probability  and  the  Monte  Carlo  Method  is  presupposed. 
Fundamentals  of  War  Gaming:  by  P.  J.  McHugh  [l]  and  A  Survey 
of  Historical  Developments  in  War  GamJ^n^  by  John  Young  ["2 J 
are  recommended  to  the  reader  who  has  not  studied  manual 
war  gaming  methods  or  is  unfamiliar  with  the  historical 
development  of  war  games. 

The  particular  method  chosen  by  the  author  to  reach 
the  stated  objective  consists  of  a  sample  computer  Simula^ 
tion  and  a  proposed  plan  for  designing  and  implementing  a 
computer  simulation  model.   The  proposed  simulation  design 
methodology  is  essentially  one  solution  to  the  major  problems 
confronting  a  novice  computer  simulation  designer. 

The  documentation  and  computer  code  for  a  helicopter 
anti-submarine  warfare  search  simulation  are  included  in  Chapter 
VI  and  Appendix  A  respectively.   This  simulation,  subsequently 
referred  to  as  AHS-I9  is  cited  as  an  example  throughout  this 
thesis.     AHS  1  was  written  as  an  example  because 


It  treats  an  operation  of  general  Interest  in  the  U.  So  Navy, 
i.e.,  a  search  for  an  evading  submarine,  and  is  of  limited 
enough  scope  to  be  easily  understoodo   This  particular  sub- 
ject was  selected  as  a  result  of  the  author's  past  experience 
as  an  Anti-Submarine  Warfare  Officer  In  helicopter  squadrons- 
AHS-1  is  intended  as  an  analytic  tool  for  the  evaluation  of 
the  relative  effectiveness  of  different  helicopter  search 
plans  in  similar  environmental  and  tactical  circumstance So 
It  is  hoped  that  units  within  the  U,  S.  Navy,  whose  responsi- 
bility it  is  to  conduct  such  analysis,  will  consider  using 
the  simulation;  either  as  written  or  with  any  necessary  modi- 
fications.  With  this  application  in  mind  it  should  be  noted 
that  a  hypothetical  distribution  of  helicopter  figure  of  merit 
was  incorporated  in  the  simulation.   This  allowed  the  thesis 
to  be  unclassified.   Rewriting  this  particular  portion  of  the 
coding  is  a  relatively  minor  modif Icationo 
1.2  Definitions. 

Some  confusion  exists  throughout  the  literature  in  the 
definition  of  terms  used  in  war  gaming  and  simulation  in  gen- 
eral »   This  is  perhaps  partly  attributable  to  the  diverse 
backgrounds  of  the  individuals  engaged  in  the  field.   War 
gamers  include  people  from  nearly  every  branch  of  science  and 
engineering  as  well  as  professional  military  men.   This  con- 
fusion is  possibly  also  due  in  part  to  a  tendency  to  attach 
special  meaning  to  certain  terms  according  to  the  particular 
activity  in  which  the  individual  is  engaged.   The  following 
definitions  will  be  adhered  to  in  the  remainder  of  this  thesis 
and  were  chosen  for  simplicity  as  well  as  generality. 


V/ar  Game  has  been  applied  to  activities  ranging  in  com- 
plexity from  two  opponents  at  a  sand  table  to  strike  exer- 
cises involving  thousands  of  men  and  millions  of  dollars 
worth  of  equipment.   Ao  'iT..  Pennington  chooses  to  leave  the 
term  undefined  because  of  the  almost  universal  disac;:reenient 
on  a  definition.  [3  J   P.  Jo  McHugh  defines  a  war  game  as 
"a  simulationy  in  accordance  with  predetermined  rules>  data, 
and  procedures,  of  selected  aspects  of  a  conflict  situationc" 
l^lj  This  definition  is  favored  by  the  author  because  it  is 
clear  and  concise  and  does  not  tend  to  imply  any  special  type 
of  war  game,  size  or  mode  of  play.   The  words  "conflict  sit- 
uation" should  be  particularly  noted.   This  is  the  primary 
distinction  between  a  game  and  other  types  of  simulationo   In 
accordance  with  current  usage,  war  game  will  not  include  ex- 
ercises in  which  real  forces  or  equipment  are  employed o 

Because  of  its  size  and  complexity,  the  real  world,  or 
even  selected  aspects  of  the  real  worldv  is  usually  impossible 
to  study  in  detail  in  the  laboratory.   Consequently,  scientists 
must  usually  resort  to  the  study  of  a  model  of  a  real  life 
situation.   This  model  may  be  iconic »  such  as  a  model  airfoil 
constructed  by  the  aeronautical  engineer  for  study  in  a  wind 
tunnel,  or  it  may  be  symbolic »   A  stick  and  gumdrop  model  of 
an  atom  is  a  symbolic  model,  as  is  the  familiar  equation  P  =  ma. 
The  latter  is  an  example  of  a  mathematical  model  and  together 
with  several  other  equations,  or  models,  is  the  model  for 
Newtoniin  mechanicso   In  essence,  a  model  is  a  representation 
of  some  aspect  of  real  life  on  which  experiments  can  be  per- 
formed; and  from  the  behavior  of  the  model  certain  laws  or 


[reneral  rules  governing  the  real  world  situation  may  be  in- 
ferred.  It  1g  characteristic  of  every  model  that  the  result- 
ing conclusions  about  the  real  world  will  never  be  any  more 
valid  than  the  model  itselfo   This  thesis  will  be  concerned 
with  mathematical  models. 

Mathematical  models  can  be  further  classified  as  deter- 
ministic or  stochastic o   The  equation  P  =  ma  is  a  determin- 
istic model  in  that  if  mass  and  acceleration  are  Icnovm  the 
force  is  uniquely  determined o   If,  however,  the  measurement 
of  acceleration  is  subject  to  certain  random  errors  the  above 
equation  could  be  written  as 

P  =  m(a+e) 
where   a  is  the  measured  value  of  acceleration  and   e   re- 
presents a  random  measurement  error«   In  this  example,   P  is 
a  random  variable,  and  the  equation  is  a  stochastic  model 
representing  an  abstract  concept  which  is  called  force.   To 
Illustrate,  the  movement  of  the  tar^^et  submarine  in  .IHS-l 
during  any  time  interval  Is  represented  by  equations  involv- 
ing the  course  and  speed  of  the  submarine  durin/?;  this  in- 
crement of  timeo   Given  course,  speed,  and  time  of  travel j 
the  distance  traversed  and  far  direction  can  be  predicted  ex- 
actly.  The  course  Itself  though,,  is  assumed  to  be  a  random 
variable  distributed  uniformly  on  the  interval  I  0,360)  degreeso 

COURSE  =  (OOO+RN) 
where;  HN  is  a  random  niAmber  drawn  from  a  set  of  numbers  dis- 
tributed uniformly  between  zero  and  360.   Therefore,  the  di- 
rection traveled  during  the  time  interval  is  a  random  variable 
and  the  model  is  stochastic.   Once  the  random  number  has  been 


dravm,  COURSj^  "beooriies  the  observed  value  of  a  randooi  variable 
and  the  direction  of  travel  is  uniquely  determined. 

The  word  si-.nulation  is  used  interchant^eably  with  model 
by  many  writers  in  the  field-   We  choose  to  define  a  dis- 
tinction between  the  twoo   As  u?ed  in  this  paper,  a  simula- 
tion includes  both  a  model  and  its  lo'^ic.   Consider  the  model 
as  a  plan  which  is  essentially  static  until  put  into  effecto 
The  simulation  is  the  dynamic  Implementation  of  this  plan,  and 
includes  the  program  logic  and  the  completed  coniputer  instruc- 
tionso   When  the  term  simulation  is  used  in  this  thesis,  com- 
puter simulation  will  be  implied  unles<^  the  content  makes  it 
clear  that  this  is  not  the  case^   A  computer  v:ar  -^ame  is  one 
particular  type  of  simulationo 

A.  S.  Householder  defines  the  Monte  Carlo  method  as  "«o. 
the  device  of  studying  an  artificial  stochastic  model  of  a 
physical  or  mathematical  process. uo"  I  Aj   The  term  Monte  Carlo 
is  relatively  new  but  the  method  is  the  technique  of  model 
sampling  used  by  statisticians  for  at  least  sixty  years.   An 
alternative  and  perhaps  simpler  statement  of  the  above  defini- 
tion is  "the  use  of  sampling  to  estimate  the  answer  to  a  math- 
ematical problemo"  |5J   For  a  comprehensive  treatment  of 
Monte  Carlo  techniques  and  applications  the  reader  is  referred 
to  Svmpocium  on  Monte  Carlo  Methods o  edited  by  H,  A.  Meyer.  1 6] 

The  tenn"  real  world  and  system  will  be  u?ed  throuThout 
the  thesirs  and  in  this  particular  context  are  considered  to  be 
synonymous.   Both  terms  are  used  to  describe  the  particular 
part  of  the  universe  which  is  beinp;  ??tudied.   The  real  world, 
or  system,  could  comprise  the  particles  exislinir  in  the  nuclnur: 


of  an  ator,  or  it  mi.'^ht  consist  of  a  portion  of  the  ooeaji,  a 
submerged  submarine  and  a  flight  of  sonar  equipped  helicopters « 
In  the  latter  example  the  weather,  sonar  operators  and  pilots, 
or  anything  else  which  might  Interact  with  the  helicopters 
or  submarine  during  a  search  operation  would  be  part  of  the 
system. 


OH/iPT-^R  II 
aOKPUTSR  SliUI^ATION  DESIGN 

The  first  prlnoiple  of  conputer  simulation  de^i.-^n  Ig 
that  every  design  situation  is  unique.   Each  Simula  Lion  pre- 
sents problems  which  r.ust  be  solved  on  the  basis  of  the  ob- 
jectives of  the  particular  study  and  the  characteristics  of 
the  real  world  irysten  under  consideration.   It  is  believed, 
howeve'',  that  certain  factors  of  a  general  nature  nust  be 
considered  in  any  computer  simulation  de-^ign  and  that  th^re 
is  a  lor.ical  order  in  which  these  factors  should  ordinarily 
be  examined.   The  purpose  of  this  chapter  is  to  describe  a 
general  methodologic  approach  to  computer  simulation  design* 

Before  attemptins  to  describe  a  methodolo.^y  of  general 
applicability  to  computer  simulation  the  question  of  why  such 
a  procedure  is  needed  should  be  answered,   '.far  .raain.g;  has  a 
lone  history'-  of  application  to  military  problems  and  in  the 
course  of  this  history  many  principles  of  garae  desif:n  have 
evolved  and  been  recorded.   Unless  there  are  rome  basic  dif- 
ferences in  the  objectives  as  well  as  techniques  associated 
with  the  desisti  process,  these  principles  of  cla'^sical  war 
f^ane  deslcja  should  be  equally  applicable  to  computer  simula- 
tion.  In  an  attempt  to  determine  if  such  differences  do  ex- 
ist, some  of  the  factors  contributing  to  the  philosophy  of 
both  classical  war  game  aind  computer  simulation  design  will 
be  considered. 


2.1  I^ffect  of  the  t-ralnin~  role  on  classical  var  .-.ane    design* 
^far  games  began  to  achieve  the  stature  of  a  serious  busi- 
ness rather  than  a  form  of  amusement  in  the  early  19th  century. 
|2,7J  War  games  from  this  time  on  were  for  the  most  part  ex- 
tensions of  the  field  maneuvero   Field  maneuvers  provided  a 
means  of  practicing  for  war  without  the  high  costs  in  life  and 
resources  associated  with  actual  combat.   Similarly  war  games 
offered  militarists  as  well  as  prospective  military  officers 
the  opportunity  to  practice  the  business  of  war  without  field- 
ing troops  and  their  associated  equipment. 

In  both  field  maneuvers  and  manual  war  games,  the  primary 
purpose  ir;  generally  training  of  the  participants  rather  than 
analysis  of  the  effects  of  factor  interactions  on  the  outcome . 
Because  people  tend  to  learn  more  rapidly  when  placed  in  fami- 
liar surroundings,  the  game  designer  attempts  to  create  a  sim- 
ulated environment  which  will  be  familiar  to  the  players.   As 
a  result,  £;ame  factors  are  chosen  because  of  their  relation 
to  the  appearance  of  reality  rather  than  the  effect  these  fac- 
tors might  or  might  not  have  on  the  final  outcome  of  the  game. 

If  general  trends  or  cause-and-effect  relationships  be- 
tween numerical  or  qualitative  inputs  and  game  results  are 
noted,  they  are  often  considered  as  a  secondary  return  rather 
Ihsin  a  primary  game  objective.   Consequently  the  analysis  of 
results,  and  the  question  of  what  specific  factors  to  analyze, 
is  generally  not  considered  a  part  of  game  design,  and  is 
often  left  entirely  to  thf»  decision  maker  or  his  staff.  The 
nature  of  this  analysis  is  also  rt-lated  to  th*?  role  of  war 
gaming  as  a  training  device.   Throughout  the  play  of  the  came 
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a  detailed  battle  history  is  usually  ir.aintained.   This  history 
is  then  studied  from  the  standpoint  of  the  effect  of  each  in- 
dividual decision  on  the  final  outcome  of  the  game.   If  the 
game  is  repeated,  the  players  and  decision  rules  may  chan?3:ee 
In  this  way  the  participants  are  offered  the  opportunity  to 
react  to  a  wider  range  of  stimuli »   This  is  desirable  from  an 
educational  standpoint.   However,  when  war  games  are  used  for 
analytic  purposes  this  type  of  game  may  test  the  player  rather 
than  opprational  plans  or  doctrineso 
2.2  Effect  of  digital  computers  on  raame    da3i.":n. 

The  digital  computer  has  contributed  more  than  computa- 
tional speed  to  simulationo  With  this  means  of  rapid  calcula- 
tion the  analyst  has  gained  a  degree  of  control  and  repeat- 
ability not  attainable  in  manual  war  gaming.   Unlike  manual 
gaming  in  which  decision  rules  and  players  might  change  be- 
tween plays,  these  factors  do  not  vary  between  replications 
of  a  computer  simulationo   Neither  will  the  coiiputer  simulation 
results  be  affected  by  the  participants  learning  the  game  nor 
by  some  unusually  bold  stroke  of  genius  exhibited  by  one  of 
the  players.   Through  the  elimination  of  human  intervention, 
a  major  source  of  uncontrollable  experimental  error,  .and  the 
gain  In  repeatability,  computer  simulation  offers  the  analyst 
an  analytic  rather  than  a  primarily  educational  tool. 

The  automation  of  war  gaming  has  also  placed  certain  re- 
straints on  the  simulation  designer  and  these  restraints  must 
be  reflected  in  any  plan  for  simulation  design.   The  amount 
of  information  which  can  be  stored  in  the  memory  of  a  computer 


is  limltsd  and  computer  running  time  is  expensive,  on  the 
order  of  ten  dollars  a  minute  for  the  IBM  7094.   Purtherriore , 
before  the  sii;iulation  can  be  run  on  a  machine  it  must  be  re- 
duced to  a  set  of  computer  instructionso   This  list  of  in- 
structions, or  computer  program,  may  represent  a  considerable 
investment  in  time  and  money.  Lar^e  computer  simulations  such 
as  the  Air  Battle  Kodel  described  by  R.  H.  Ad.ams  and  J<.  L« 
Jenkins  I  3  I  may  be  tvro  years  or  more  in  the  construction.  [  9j 
Becaune  of  these  inherent  limitations,  it  is  ursually  aecessary 
to  exclude  certain  factors  which  do  not  contribute  substan- 
tially to  the  purpose  of  the  game  and  may  even  detract  from 
it  by  masking  the  game  response  to  factors  which  do.   The 
problem  confronting  the  simulation  designer  is  not  Just  how 
to  achieve  a  high  fidelity  representation  of  the  real  world, 
but  how  to  selectively  choose  those  components  of  the  real 
world  which  characterize  the  particular  questions  to  be 
answered,   A  major  effect  of  digital  computers  on  game  design 
is  that   a  need  has  arisen  for  criteria  which  v;ill  help 
the  designer  in  deciding  what  is  and  is  not  pertinent  to 
the  game  purpose. 
2.3  Ob.lectives  of  a  computer  desip:n  methodology » 

The  preceding  discussion  suggests  that  computer  war  gam- 
ing entails  something  more  than  adapting  manual  games  for  play 
on  high  speed  computing  machineso   It  seems  reasonable  then 
that  coiaputer  game  design  should  embody  more  than  a  straight- 
forward extension  of  mcinual  war  game  techniques.   The  fore- 
going also  indicates  some  of  the  requirements  a  plan  for  sim- 
ulation design  should  satisfy-   These  are? 
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(1)  The  simulation  should  be  designed  -rith  specific 
questions  in  mlnde 

(2)  The  de'^ign  should  capitalize  on  the  advantages  of 
repeatability  and  precise  parameter  control  offered 
by  high  speed  computers..   This  implies  utilizing; 
some  means  of  reducing  the  effect  of  experimental 
error. 

(3)  ^pha^is  chould  be  placed  on  diiiinishing  the  dis- 
advantages of  limited  space  ;ind  high  cost  associated 
with  digital  computers, 

(4)  Factors  should  be  chosen  on  the  basis  of  their  re- 
lation to  simulation  objectives. 

(5)  Planning  of  the  analysis  and  interpretation  of  re- 
sults should  be  concurrent  with  designing  the  sim- 
ulation. 

2.4     Sxp'srimental   desi-Tn  approach   to    siraulation  . 

It  is  considered  that  the  above  requirements  are  best  met 
by  the  methods  of  experimental  design.   This  is  not  surprising 
since  a  simulation  is  basically  a  sampling  experiment.   The 
principle  departure  of  simulation  from  the  agricultural  ex- 
periments, which  gave  experimental  design  its  impetus,  is 
founded  in  technique  and  degree  of  abstraction  rather  than 
theory.   In  simulation  the  experimenter  first  builds  a  re- 
presentation or  model  of  the  real  world  and  then  performs  his 
experimentation  on  this  model.   The  agricultural  worker  may 
use  real  plants  grown  in  real  soil  as  the  basis  of  his  experi- 
ment, but  not  just  any  plants  grown  on  Just  any  soil.   The 
plants  are  carefully  chosen  and  plots  are  arranged  in  a  se- 
lected pattern  in  order  to  control  the  effect  of  extraneous 
factors  such  as  soil  variability  or  differences  in  drainage. 
Similarly,  the  simulation  designer  selects  factors  according 
to  whether  they  contribute  to  the  essence  of  realism  in  a 
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specific  situation.   The  agricultural  experiraentor  eliminates 
weeds  which  might,  by  taking  water  from  the  plants,  result 
in  an  erroneous  estimate  of  the  characteristic  yield.   A 
problem  confronting  the  simulation  desir;ner  is  that  he  must 
first  determine  which  are  "plants"  and  which  "weeds".   Once 
this  determination  has  been  made  he  must  resist  the  tempta- 
tion to  simulate  the  "vreeds"  Just  because  they  are  part  of 
the  real  world  he  is  dealing  with. 

A  further  motivation  for  the  use  of  the  experimental 
design  approach  is  suggested  by  I .  P.  Wa'^ner. 

There  have  been  times  when  the  inexperienced  VTar  Gaming 
researcher  has  set  up  his  own  study,  planned  his  oirm  runs 
and  after  the  runs  have  been  made,  offered  his  data  for 
someone  to  interpret.   The  result  is  usually  an  inability 
to  place  proper  interpretation  upon  the  results,  nec- 
essary runs  have  been  omitted,  redundant  and  superfluous 
runs  have  been  conducted  and  the  time  and  money  expended 
may  exceed  the  value  received.  [  10  ] 

Implicit  in  this  approach  is  that  selection  of  factors 
to  Include  in  the  model  is  closely  related  to  data  collection. 
The  collection  of  data  is  as  important  a  part  of  statistical 
analysis  as  data  preparation  and  interpretation. 

For  those  who  may  not  be  familiar  with  basic  experiment- 
al design  the  following  from  Design  and  Analysis  of.  Industrial 
Expgrinents.  edited  by  0.  L.  Davie s,  should  further  clarify 
the  relation  between  experimental  design  and  the  discussion 
to  follovr. 

The  first  step  in  planning  any  experiment  is  to  form  a 
clear  picture  of  objects  of  the  experiment,  the  factors 
which  may  affect  the  results,  and  the  errors  v:hich  will 
inevitably  arise.   The  chief  objects  of  the  experimental 
deslfin  are: 

(1)  To  arrange  the  experiment  so  that  the  effect  of 
changiing  each  condition  can  be  readily  measured 
and  separated  from  the  effects  of  changin;^  the 
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other  conditions,  and  from  experimental  error. 

(il)   To  obtain  a  valid  estimate  of  error  appropriate 
for  assessing  the  statistical  significance  of 
the  effects  of  the  factors. 

(ill)   To  enable  the  effects  to  be  measured  with  the 
required  accuracy,  LnJ 

The  preceding  principle  of  experimental  design  vfill  play 
an  important  role  in  determining  the  course  to  be  followed  in 
designing  a  simulation.   The  first  step  will  be  to  form  a 
clear  picture  of  the  objects  of  the  experiment  and  a  first 
approximation  of  the  factors  affecting  the  results.   This  is 
treated  in  ::hapter  III,  and  vrill  take  the  form  of  a  prelimin- 
ary study  of  the  simulation  purpose  and  the  system  to  be  sim- 
ulated. 

Consideration  of  the  sources  of  error  and  the  first  ob- 
jective of  experimental  design  stated  above  will  shape  the 
main  body  of  the  proposed  simulation  design  method.   The  nec- 
essity of  designing  the  simulation  so  that  the  effects  of 
varying  parameters  can  be  measured  and  separated  from  the  ef- 
fects of  errors  and  each  other  was  seen  to  be  a  basic  departure 
from  classical  war  game  requirements.   This  objective  will  be 
realized  in  the  computer  simulation  by  designing  the  model  it- 
self as  a  factorial  experiment,  and  is  the  subject  of  Chapter 
IV.   The  second  and  third  objectives  of  experimental  design 
are  intrinsic  to  model  testing^  and  will  be  discussed  in  con- 
Junction  with  Implementation  of  the  model  in  Chapter  V. 

Before  proceeding  to  the  preliminary  study  phase  a  brief 
discussion  of  factorial  experiments  is  in  order.   This  class 
of  experimental  designs  wil  constitute  the  pattern  for  the 
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simulation  model. 

2.5   Factorial  experiments. 

Factorial  experiments  represent  a  class  of  experiments 
which  are  particularly  applicable  to  studying  variation  de- 
liberately introduced  by  the  experimenter.  Ill,12j  The  specific 
techniques  associated  with  constructing  these  experiments  is 
beyond  the  scope  of  this  thesis  and  will  not  be  considered. 
The  intention  here  is  to  consider  certain  very  basic  concepts 
of  experimental  design  rather  than  techniques.   The  purpose  of 
this  discussion  is  to  indicate  how  these  concepts  can  help  the 
simulation  model  designer  in  selecting  inputs  which  reflect 
the  objectives  of  the  system  as  they  pertain  to  the  simulation 
purpose . 

The  following  terms  will  be  used  throughout  the  remainder 
of  the  discussion  of  factorial  experiments  and  simulation 
design. 

Factor.   In  general  a  factor  is  any  quantity  which  con- 
tributes to  the  determination  of  the  state  of  the  simulated 
system,  and  can  be  controlled  by  the  experimenter.   Examples 
from  simulations  of  military  operations  are  number  of  aircraft, 
weather  conditions,  or  electronic  counter-measures. 

Factor  Levels.   The  Values  assigned  to  a  factor  are  term- 
ed levels.   A  level  may  represent  a  number  such  as  thirty 
bombers  or  it  may  be  qualitative.   Jamming  or  not  Jamming 
would  represent  two  levels  of  electronic  counter-measures, 
a  qualitative  factor. 

Treatment.   A  treatment  refers  to  the  level  of  the  factors 
which  describe  one  simulation  run.   A  convenient  notation  for 
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a  treataent  is  (A.BpO.).   This  indicates  a  simulation  run  with 
factor  A  and  C  at  their  lower  levels  and  factor  B  at  the  second 
level. 

Response.   This  is  synonymous  with  simulation  results  as 
we  have  used  this  term.   A  treatment  response  represents  the 
result  of  a  trial  of  an  experiment  corresponding:  to  that 
treatment. 

Main  effect.   A  main  effect  for  a  factor  at  two  levels 
is  the  average  difference  in  response  due  to  the  change  in 
levels. 

Interaction.  An  interaction  is  said  to  occur  when  the 
response  to  changes  in  one  factor  depends  on  the  particular 
level  of  another  factor.  An  example  might  be  artillery  and 
spotters.  Increasing  the  number  of  firing  piecer  might  not 
increase  the  damage  unless  additional  spotters  were  assigned 
to  direct  their  fire.  Similarly  increasing  the  number  of 
spotters  might  have  no  appreciable  effect  by  itself. 

Sample .   A  sample  consists  of  a  series  of  plays  which 
differ  only  as  a  result  of  the  series  of  random  numbers  used. 
Sainples  consisting  of  more  than  one  play  are  generally  used 
in  simulaions  whose  output  is  either  success  or  failure  rather 
than  a  detailed  battle  history,   \ihen  a  simulation  run  con- 
sists of  more  than  one  sample  the  additional  samples  represent 
replications  of  the  same  treatment  combination.   Replication 
provides  the  means  of  estimating  experimental  error. 

The  factorial  experiment  was  chosen  as  the  pattern  for 
simulation  design  because  it  is  generally  the  most  efficient 
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method  of  ntudyins  the  effect  on  the  results,  \rhcn   varyln^c' 
tvro  or  r.iorc  input  quanLiGie--.  Illl   .in  efficient  method  is 
one  which  yields  the  required  infor,;iation  with  r,he  uiiniinum 
effort.   In  essence  the  deci.Ti  ajid  analysis  of  a  factorial 
experiment  consists  of:   (1)  conductin/v;  trials  of  the  e::perl- 
:.ient  with  combinations  of  several  factors  at  more  than  one 
level,   (2)  obGorvin-  the  response  to  each  treatment  and 
determining  estimates  of  main  effects,  factor  interactions  and 
experimental  error,   (3)  te;^tlnr;  hypotheses  to  determine  \rhioh 
of  the  effects  are  real  and  which  are  attributable  to  p-'-'per- 
imental  error. 

\xi   individuil  experiment  on  a  completed  simulation  liiay 
be  conducted  ac  a  complete  factorial  experiment  in  vrhich  ob- 
servations are  nade  wilh  every  factor  at  every  one  of  its 
levels.   Thin  U'^ually  requires  a  jjreat  iLOny  computer  runs, 
however,  and  some  form  of  fractional  factorial  e:^perimentation 
ni:j;ht  be  e^iployed.   The  essential  point  is  that  by  building 
all  of  the  factorrj  vrhich  are  pertinent  to  the  stu:1y  purpose 
into  the  model  from  the  very  beginnlnr;,  the  capabili-^y  of 
efficiently  examining  anj-  combination  of  these  factors  and 
factor  levels  will  exist  in  the  completed  simulation.   In 
this  way  additional  assurance  that  the  analysis  will  reflect 
the  study  purpose  will  have  been  £;ained« 
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CHAPTBR  III 
PR^LIITIN-IRY  STUDY 

The  design  of  a  computer  game,  as  Tvith  essentially  all 
operations  research  studies,  must  begin  vrith  a  statement  of 
the  problem.  In  most  caces  the  initial  statement  of  the 
problem  will  come  from  the  person  or  persons  requesting  the 
study.  The  concern  here  is  more  with  the  interpretation  of 
the  problem  by  the  simulation  designer. 
3.1  Fornulation  of  the  problem. 

The  first  -^tep  toward  formulating  the  problem  is  to  be- 
come familiar  with  the  various  components  of  that  portion 
of  the  real  v/'orld  from  which  the  problem  comes.   Components 
in  this  context  may  include  hardware  r.uch  as  weapon  systems, 
operating  doctrine,  operational  environment  and  a  myriad  of 
other  items  which  make  up  the  particular  aspect  of  the  real 
vrorld  under  consideration.   It  is  suggested  that  a  block  flow 
diagram  showing  the  functional  position  of  each  component  in 
the  system  be  constructed  as  part  of  this  initial  research. 
A  flow  diagram  is  also  useful  for  showing  interactions  among 
components. 

The  determination  of  the  scope  of  the  proposed  design 
should  be  concurrent  with  the  preceding  step.   As  the  analyst 
familiarizes  himself  with  the  system  components,  the  size 
and  compl'^xity  of  the  problem  area  will  become  more  apparent. 
Conversely,  unless  he  has  some  feel  for  the  range  of  activity 
included  in  the  problem  it  is  impossible  to  know  vrhat  the  cot- 
ponentrj  of  the  system  are. 
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If  the  design  is  to  be  built  around  particul?.r  que  rations 
as  related  to  the  objjectives  of  the  systeri,  these  questions 
must  be  clear  from  this  point  on»   It  is  thereforo  prudent 
to  compare  the  interpretation  of  objectives  with  stated  ob- 
jectives.  It  is  possible  that  on  the  basis  of  the  above 
initial  research,  the  original  problem  may  be  modified  or 
expanded,  or  another  research  method  chosen. 
3 '2   Applicability  of  computer  simulation  to  the  nroblem. 

Once  the  scope  of  the  problem  is  kno^m  and  the  designer 
has  gained  some  understanding  of  the  system,  it  nust  be  de- 
cided whether  or  not  computer  simulation  techniques  fit  the 
situation.   This  is  a  step  which  is  perhaps  too  often  omitted 
in  simulation  design.   The  following  statement  by  0 .  J.  Thomas 
concerning  unfortunate  practices  in  military  gaming  is  equally 
applicable  to  computer  simulation. 

...there  is  an  unfortunately  common  tendency  to  think 
of  "military  gaming"  as  being  not  only  universally 
applicable  but  also  mechanically  applicable.   The  failure 
to  match  technique  with  problem  is  often  expensive  and 
almost  invariably  produces  disappointment.   [l3] 

Consideration  of  theoretical  model.   It  is  generally 
considered  that  simulation  should  be  used  only  if  it  is 
impossible  to  solve  the  problem  using  a  simple  theoretical 
model.  [  14,15,16  1   In  applying  this  generality  it  is  important 
that  \ie   do  not  lose  sight  of  the  given  problem  in  trying  to 
fit  the  system  to  an  analytic  model o   Any  situation  can  be 
treated  analytically  if  it  is  simplified  enough.   However 
oversimplification  of  the  model  and  the  subsequent  analysis 
may  invalidate  the  results  of  the  study.   Consider  the  heli- 
copter search  situation  treated  by  .U^S-l.   The  system  itself 
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is  not  overly  complicated,  and  analytic  models  of  a  search 
helicopter  versus  an  evading  submarine  conflict  have  been 
built.   One  of  the  helicopter  search  plans  currently  opera- 
tional with  the  U.  S.  Navy  is  based  on  a  game  theory  solution 
of  this  conflict  situation.   The  model  used  in  this  solution 
includes  the  assumption  that  submarine  maximum  speed  is  known, 
and  that  helicopter  sonar  ranges  and  dip  cycle  time  are  known, 
and  the  same  for  all  helicoptersc   The  simulation  study,  AHS-1, 
was  motivated  by  the  fact  that  these  factors  among  others  are 
subject  to  uncertainty.   The  author  wanted  to  know  the  rela- 
tive effectiveness  of  certain  search  plans,  including  the  so 
called  optimum  plan  found  analytically,  in  the  face  of  this 
uncertainty. 

Expected  number  of  plays.   A  second  consideration  in  de- 
ciding when  to  resort  to  computer  simulation,  is  the  number 
of  expected  repetitions  of  the  experiment.   The  validity  of 
the  Monte  Carlo  technique  is  dependent  on  repeated  trials.   It 
is  not  unheard  of  to  make  ten  or  more  consecutive  sevens  with 
two  fair  dice.   Inferring  that  this  is  routine  on  the  strength 
of  observing  the  phenomenon  once,  however,  could  be  financially 
disastrous.   Similarly,  confidence  in  simulation  results  is 
related  directly  to  the  number  of  replications. 

The  number  of  input  parameters  to  be  examined  and  the 
range  of  interesting  values  of  these  parameters  also  affects 
the  number  of  required  repetitions.   The  effect  on  experimental 
results  of  varying  parameter  values  is  usually  of  more  interest 
than  the  particular  outcome  corresponding  to  any  one  set  of 
input  values.   In  addition  certain  parameters  may  exhibit 
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interaction.   That  is  the  effect  of  changing  the  input  value 
of  tv;o  or  more  factors  simultaneously  is  considerably  dif- 
ferent from  that  of  varying  these  same  parameters  individ- 
ually. 

It  is  apparent  that  if  several  parameters  are  to  be  ex- 
amined at  more  than  one  level,  each  run  consisting  of  many 
replications,  the  number  of  simulation  runs  can  become  very 
large.   Even  vrith  a  game  of  relatively  limited  scope  such  as 
AHS-1,  several  months  might  be  required  to  examine  all  input 
parameters  at  only  a  few  levels  by  manual  methods.   On  the 
other  hand,  if  only  a  fevr  plays  of  a  game  are  contemplated, 
computer  simulation  is  probably  not  the  most  efficient  method. 
A  completed  computer  program  represents  a  sizeable  investment 
and  the  cost  in  time  and  money  of  running  this  program  on  a 
computer  one  time  represents  a  small  fraction  of  this  invest- 
ment. 

Player  trainin/^.   The  potential  value  of  training  derived 
from  player  participation  in  the  game  must  also  be  considered. 
If  the  primary  purpose  is  clearly  player  training  then  a 
manual  game  is  indicated,   l^niether  the  purpose  is  educational 
or  analytic  is  not  always  tacitly  obvious,  however,  and  a 
decision  to  satisfy  one  or  the  other  requirements  will  have 
to  be  made. 

Assuming  that  computer  simulation  is  theoretically  ap- 
plicable the  analyst  must  next  ask  if  it  is  physically  prac- 
ticable.  A  study  which  is  worth  doing  is  worth  doing  correct- 
ly.  If  the  limitations  of  time  and  available  phisical  faci- 
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lities  preclude  an  adequate  treatment  of  the  problem  by  sim- 
ulation methods  the  person  requesting  the  study  must  be  ap- 
prised of  this. 
3.3  Research  into  existinf;  simulations. 

Having  decided  to  use  simulation  methods,  the  designer 
night  next  ask  if  there  is  an  existing  model  vrhich  is  suited 
to  this  particular  situation.   However,  in  many  respects  try- 
ing to  fit  the  problem  to  an  existing  model  is  contradictory 
to  the  stated  objective  of  tailoring  the  simulation  to  the 
particular  study  purpose.   It  would  be  preferable  to  re- 
state the  question  more  generally  asj  what  is  the  existing 
level  of  simulation  knowledge  which  may  be  related  to  this 
system?   Although  it  might  be  unwise  to  modify  the  particular 
situation  to  fit  an  existing  model,  it  is  equally  foolish  not 
to  profit  from  the  effort  of  others  whenever  possible. 

To  illustrate  this  point  we  refer  to  the  scheme^  by  v;hich 
sonar  detection  ranges  are  computed  in  AHS-1.   The  particular 
objective  required  some  more  realistic  means  of  determining 
detection  ranges  than  the  usual  procedure  of  using  a  determ- 
inistic assured  sonar  range.   The  scheme  selected  was  design- 
ed by  R.  L.  Klinkner  as  a  general  method,  applicable  to  most 
existin::  Navy  r^onars.  |17|   Furthermore  an  algorithm  for  com- 
puting the  needed  sonar  ranges  by  this  method  had  been  "lirit- 
ten  and  coded  in  the  computer  language  previously  selected 
for  AHS-1.   The  availability  of  this  computer  coding  made  it 
possible  to  design  a  simulation  which  fit  the  requirements  of 
the  problon  within  the  time  allotted.   This  might  well  have 
been  impossible  othenrise. 
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OH;lPT'^R  IV 
DESIGN  OF  Tin  MODEL 

The  model  of  Inlerest  in  this  chapter  is  a  stochastic 
mathematical  model.   .Recall  that  a  mathematical  model  was 
defined  as  a  mathematical  representation  of  nome  aspect  of 
the  real  world.   This  model  can  usually  be  subdivided  into 
models  of  each  of  the  individual  conponents  or  activities 
which  characterize  the  irystem  of  interest.   It  is   partly 
because  of  this  subdivision  that  the  construction  of  a  block 
flou  diagram  wa:^  nugf^ented  as  part  of  the  for.nulation  of  the 
problem.   Such  a  flow  chart  is  often  a  useful  aid  to  recog- 
nition of  logical  functional  subdivision? . 
4.1   Functional  coiiipouents  of  the  model. 

There  are  several  reasons  for  subdividing  the  model  in- 
to functional  couponents,  the  degree  to  which  this  is  done 
bein-^  closely  related  to  the  scope  of  the  system  under  study. 
Any  system,  even  one  of  limited  complexity  such  as  the  helicop- 
ter-submarine conflict  of  our  example,  is  much  simpler  to 
visualize  if  it  ir:  broken  doim  into  functional  components  or 
activities. 

A  second  ar£;u].ient  for  subdividing  the  model  alon^  func- 
tional lines  is  that  for  more  complex  problems  several  persons 
may  be  simultaneously  engaged   in  modeling  the  system.   Con- 
sider a  simulation  of  a  large  area  anti-submarine  operation. 
Helicopter  operations  might  constitute  one  logical  subset 
to  be  allotted  to  one  individual.   One  person  could  develop 


this  part  of  the  sir.ulation,  code  it  and  complete  the  program 
check-out  independently  of  the  main  study  ^roup.   The  example, 
AHS-1,  could  with  only  minor  modification  represent  such  a 
subsection  of  the  larger  simulation. 

Dividing  the  model  into  such  subsets  also  facilitates 
parameter  validation  and  program  check-out.   This  applies 
equally  whether  the  simulation  is  being  built  by  one  person 
or  several. 

The  preceding  model  subdivisions  are  primarily  for  mech- 
anical or  operational  convenience  of  the  designer.   They  are 
aids  to  organizing  the  existing  knowledge  about  the  system  so 
that  the  model  will  become  physically  manageable.  Models  may 
also  be  subdivided  into  conceptual,  or  structural,  components. 
4.2   Structural  model  components. 

The  structural,  or  conceptual,  components  of  a  model  ares 

(1)  Assumptions 

(2)  Input  parameters 

(3)  Input  variables 

(4)  Algorithms 

(5)  Measures  of  effectiveness 

It  will  be  simpler  to  define  the  above  terras  and  explain 
why  they  are  considered  to  be  useful  to  the  model  designer  if 
the  following  two  definitions  are  agreed  upon.  These  defini- 
tions are  from  Van  Nostrand's  Scientific  Encyclopedia.  L  ■^'^J 

Parameter.   An  arbitrary  constant,  as  distingaish2d  frou 
a  fixed  or  absolute  constant.  Any   desired  numerical  value 
can  be  ':iven  to  a  parameter. 

V^irlablq .   A  quantity  as  distinguished  from  a  constant 
to  which  any  number  of  values  in  a  given  set  may   be  assign- 
ed.  If  the  set  comprises  a  doraain(a,b),  then  the  variable 
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is  only  defined  over  this  inte.^v-il  and  all  values  out- 
side of  the  interval  are  ignored. 

These  definitions  vrill  apply  vrhenever  parameter  or  variable 

is  used  in  the  succeedins  discussion. 

Assumptiors.  Any   conclusions  arrived  at  on  the  basis 
of  simulation  results  must  be  considered  invalid  unless  they 
were  made  with  a  complete  Icnowledse  of  the  major  model  assump- 
tions.  The  model  assumptions  reflect  more  than  any  other 
factor  the  designers  understanding  of  the  system  to  be  modeled 
and  his  interpretation  of  the  objective  of  the  simulation. 
The  assumptions  reflect  what  is  known  precisely  about  the  real 
world  as  well  as  vrhat  is  uncertain  and  the  degree  of  uncertain- 
ty.  To  illustrate,  the  following  are  model  assumptions  from 
AHS-1. 

The  ocean  is  represented  as  a  rectangular  coordinate 
system  in  the  simulation.   This  is  based  on  the  assumption 
that  the  area  encompassing  a  helicopter  search  is  small  enough 
that  the  earth's  curvature  need  not  be  considered.   Helicopters^ 
which  are  characterized  as  points  from  which  associated  sonar 
ranges  are  measured »  are  assumed  to  move  about  this  square 
grid  in  discrete  Jumps.   These  assumptions  reflect  certain 
well  defined  facts  about  helicopter  sonar  search  operations. 

It  is  further  assumed  that  the  above  sonar  ranges  vary 
among  helicopters  according  to  some  probability  distribution, 
and  that  this  distribution  Is  a  function  of  sonar  figure  of 
merit.   It  is  known  that  sound  propagation  conditions,  at 
any  one  point  in  the  ocean,  are  a  function  of  a  time  correlated 
random  variable.  [l7>19  J   This  time  dependence  of  sonar 
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ranges  has  not  been  Included  in  the  model  for  two  reasons'. 
First,  a  helicopter  search  occupies  a  relatively  short  time 
span  and  the  variation  due  to  the  above  random  variable  is 
expected  to  be  small  in  the  time  interval  of  interest. 
Secondly,  sonar  ranges  can  be  expected  to  decrease  or  in- 
crease with  the  same  relative  frequency  during  a  search.   If 
for  some  reason  sonar  ranges  were  biased  so  that  either  a 
decrease  or  an  increase  with  time  could  be  expected  this 
assumption  might  not  be  valid*. 

The  assumptions  concerning  sonar  range  all  have  a  bear- 
ing on  uncertainties  in  nature  and  the  effect  of  those  un- 
certainties on  the  simulation  purpose.   All  of  the  preced- 
ing assumptions  reflect  the  model  builder's  understanding 
of  basic  functional  relationships  of  the  system  components 
as  they  are  related  to  the  simulation  purpose. 

Input  parameters.   These  are  constants  which  may  be 
assigned  any  arbitrary  value  by  the  simulation  user.   Our 
interest  in  input  parameters  is  generally  in  the  way  changes 
in  parameter  values  affect  the  outcome  of  the  operation  sim- 
ulated.  The  factor  which  characterizes  these  quantities  as 
parameters  rather  than  variables  is  the  degree  of  human  con- 
trol in  the  real  world  situation  being  simulated.   For  example, 
the  particular  search  plan  designated  by  a  helicopter  flight 
leader  is  entirely  within  the  control  of  the  flight  leader 
in  questionv   In  assigning  values  to  this  particular  input 
parameter  the  analyst  is  concerned  with  finding  criteria  by 
which  a  certain  plan  should  be  selected »  given  a  specific 
state  of  nature o  Note  that  input  parameters  need  not  be 
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nu2ieric?.l  quantities  nor  'Till  they  nece'^narily  be  the  actual 
inputs  to  th'^  pr0  3ra'n.   Search  plan  is  an  inp^^"  '^arar.eter 
and  itn  value  nii^iht  be  Plan  Lisa  or  Spiral  Search  Three. 
This  quantity  will  ^ive  rise  to  input  nuantitieG  vhioh  do 
hav^"  numerical  values,  however.   In  this  e::a!rple  bearin-s 
and  distances  of  the  search  le^s  are  such  numerical  valued 
quantities. 

In^n-^  varinbles.   These  are  input  quantities  over  which 
we  cannot  expect  to  have  control  in  actual  practice,   -^or 
example  submarine  speed,  sonar  conditions  or  tirn.e  late  at  a 
submarine  datum.   Input  variables  are  characteristic  of  xhe 
systen  and  the  ran::;e  of  values  assif^ned  to  the  to  quantities- 
is  essentially  determined  by  nature 

As  noted  above  for  input  parai^eters,  input  variables  as 
initially  defined  by  the  nodel  desi~nor  may  not  be  the  actual 
pro^ra"-  inputs.   In  most  cas3s  these  vrill  be  determined  In 
conjunction  with  the  detailed  inve^^ti^^ation  and  data  collec- 
tion inrrin'iic  co  the  actual  nodel  conclruction. 

To  tie  the  above  definitions  to:;ether,  consider  the  func- 
tional relation  fro.n  AHS-lj 

where  R     =   Output  (Nur.iber  of  defections  during  a  run) 

3^  =  Submarine  Speed 

'  Input  variables 
T-  =  Tine  late  of  helicopter  at  datum 

p  =  Search  Plan  A 
^*-  )>  Input  naranetPTr; 

H„  =  Number  of  helicopters 

N 
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F  is  the  rule,  drteriUined  by  the  model  assU'Tiptions  and 
simulation  lo^ic ,  v:hich  determines  the  dependence  of  R^^  on 
the  input  variables  and  paranieterc.   Submarine  speed  and  time 
late  are  independent  variables  which  are  allowed  to  vary  over 
selected  portions  of  their  apsuned  ranges.   The  input  para- 
meters may  be  assigned  any  arbitrary  set  of  values.   V/hereas 
the  range  of  input  variables  depends  on  what  the  analyst  has 
assumed  to  be  the  possible  states  of  nature,  input  parameters 
reflect  particular  questions  about  the  system.   The  optinun 
values  of  these  parameters  corre spending  to  specific  states 
of  nature,  i.e.,  specific  values  or  ranges  of  values  of  the 
input  variables,  represent  part  of  what  we  hope  to  learn 
from  the  simulation. 

Input  parameters  and  variables  may  change  roles  accord- 
ing to  the  specific  circumstances  of  application.   For 
instance,  it  may  be  possible  to  control  time  la+eat  datum  to 
some  extent  but  at  the  expense  of  the  number  of  helicopters 
which  can  be  sent  to  the  datur..   It  might  happen  that  two 
helicopters  are  airborne  when  a  submarine  datum  is  established 
The  decision  must  be  made  to  launch  more  helicopters  and 
accept  a  gr^-ater  time  late,  or  proceed  imnediately  to  datum 
with  the  avT-ilable  aircraft.   Here  tirr.e  late  assumes  the  role 
of  a  parameter  and  number  of  helicopters  is  essentially  a 
variable. 

-U:orl  th  ■' .      "An  algorithjn   ir,  a   stated  rule   or  procedure 
that  r.ay   b.^    followed    to   arrivs   at   some    specified   goal.'Tso]     ;m 
algorittim  is  characterized  by   the   fact   that  it  may  be 
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mechanically  applied.   Ae  procedure  for  n:en?rat.in,2;  randoni 
niiEberf^  employed  in  .AHS-l  in  an  alrorithr. 

The  model  designer  must  consider  every  po?^sible  question 
which  mi^ht  arise  and  proprocram  an  answer  to  each  question* 
It  is  not  possible  to  defer  a  decision  until  such  time  during 
the  play  of  the  game  that  the  need  for  making  the  decision 
nay  arise.   Therefore,  in  order  to  build  a  model  which  can  be 
incorporated  in  a  computer  program  it  must  be  possible  to 
reduce  the  essence  of  the  real  world  situation  to  a  collection 
of  algorithms.   By  means  of  this  set  of  algorithms  the  program 
operates  on  input  variables  and  parameters  in  a  way  determined 
by  the  assumptions. 

Measures  of  effectiveness.   Saaty  defines  measure  of  ef- 
fectiveness as  "a  criterion  which  measures  the  extent  of  suc- 
cess of  a  solution,  as  related  to  the  objectives,  when  ap- 
plied to  a  problem  arising  in  an  operation."  [21I    The  prin- 
ciple measure  of  effectiveness  incorporated  into  .\HS-1  is 
number  of  submarine  detections.   In  this  example  the  objective, 
as  used  in  the  above  definition,  is  to  evaluate  the  comparative 
effectiveness  of  various  search  plans  in  similar  tactical 
and  physical  environmnents.   One  of  the  problems  in  this  ex- 
ample is  that  the  number  of  available  helicopters  is  generally 
limited. 
4,3  Determinin'^  measures  of  effectiveness. 

Military  simulation  studies  fall  into  two  broad  cate- 
gories according  to  purpose;  testing  plans  and  evaluation  of 
weapon  syst;^ms.   Regardless  of  whore  the  simulation  fits  in 
thio  general  clansirication,  however,  if  decisions  are  to  be 
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influenced  by  the  simulation  results  some  sieasur---;  of  the  suc- 
cess or  failuro  of  tho  operation  in  a  ^iven  instance  is  re- 
quired. Measures  of  effectiveness  nay  be  comparatively 
obvious,  or  quite  subtle.  In  the  case  of  AHS-1,  the  number 
of  submarine  detections  is  a  straightforward  measure  of  the 
comparative  effectiveness  of  two  search  plans  under  similar 
circumstances. 

Consistent  well  defined  measures  of  effectivenesr.  are 
essential  to  ^ood  simulation  design.   To  a  considerable  extent 
the  measure  of  v^hether  an  included  parsuneter  or  variable  is 
essential  to  the  simulation  purpose  is  the  sensitivity  of  the 
measures  of  effectiveness  to  changes  in  thir^  input.   Consequent- 
ly,  the  measures  of  effectiveness  will  represent  important 
criteria  by  which  the  analyst  must  decide  whether  a  factor 
is  to  be  included  or  omitted  from  the  simulation. 

Inconsistent  measures  of  effectiveness  indicate  incom- 
patability  of  the  simulation  objectives  as  well  as  the  lack 
of  a  clear  Interpretation  of  the  objectives  of  the  system 
bein:2  simulated.  In  essence  the  measures  of  effectiveness  are 
the  statements  of  the  objectives  of  the  system  beln^  studied. 
The  selected  measures  of  effectiveness  are  indicative  of  the 
desisners  understanding  of  the  system  objectives  Just  as 
model  assumptions  reflect  his  understanding^  of  basic  func- 
tional aspects  of  the  system. 

Output  data.   As  far  as  is  possible  the  output  of  the 
simulation  should  be  determined  in  conjunction  i^ith  selecting 
measures  of  effectiveness.   .0.1  of  the  measures  of  effective- 
ness will  be  output  quantities  and  these  may  also  cive  rise 
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to  associa-ted  nurr.bers  which  "irlll  be  of  interest  to  the 
analyst.   In  addition,  certain  control  variables  should  be 
Included  in  the  output.   These  controls  should  be  quantities 
whose  values  or  distributions  are  well  enough  known  that  they 
can  be  used  as  rough  checks  on  the  logic  and  numerical  inputs 
to  the  simulation. 
4.4  Selectin"^  inputs 

Every thine  that  has  been  accomplished  up  to  this  time 
has  been  preliminary  to  the  actual  design  of  the  model.   During 
the  formulation  of  the  problem  the  scope  of  the  simulation  was 
determined.   The  scope  in  turn  determined  the  acceptable 
level  of  aggregation  consistent  with  the  stated  problem.   Hav- 
ing chosen  measures  of  effectiveness  by  which  the  simulation 
results,  and  to  some  extent  the  simulation  itself,  will  be 
evaluated;  the  designer  is  prepared  to  select  the  inputs  and 
Incorporate  these  factors  along  with  certain  assumptions  and 
rules  into  a  representation  of  the  system. 

Inputs  were  classified  as  either  variables  or  par.araeters 
on  the  basis  of  the  extent  to  which  these  quantities  can  be 
controlled  in  the  real  world.   This  classification  V7ill  be 
useful  in  determining  the  ranges  of  input  variables  and  para- 
meters and  the  distinction  should  be  kept  in  mind  during  the 
follov/ing  discussion.   In  determining  which  inputs  to  include 
in  the  model,  however,  it  is  more  suitable  to  our  purpone  to 
combine  inputs  and  variables  under  the  single  title  of  factors. 

D.  R.  Cox  lists  four  items  to  be  considered  in  designing 
a  factorial  experiment .[  l^J    These  are,  the  factors  to 
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Include  In  the  design,  the  levels  of  each  factor,  the  number 
of  observations  to  be  made,  and  measures  to  reduce  the  effect 
of  uncontrolled  variation.   The  last  two  of  these  are  more 
closely  related  to  simulation  testing  and  analysis  of  the 
simulation  results.   The  factors  to  be  included  and  the  levels 
of  each  factor  determine  the  form  of  the  model  itself  and  will 
be  considered  at  this  time. 

It  should  be  noted  that  at  this  point  the  simulation 
designer  is  concerned  with  a  general  experimental  design.   One 
particular  experiment  might  involve  only  two  or  three  factors, 
each  at  one  or  many  levels.   However,  the  simulation  objec- 
tive often  involves  several  questions  about  the  system,  and 
the  analyst  must  design  the  simulation  to  encompass  a  class 
of  experiments.   Consequently,  all  of  the  factors  necessary  to 
conduct  any  experiment  of  this  class  must  be  built  into  the 
model . 

Choosing:  factors,  ^fhether  a  particular  factor  should  be 
included  in  the  simulation  is  primarily  a  matter  of  experience 
and  Judgment.   Once  a  working  version  of  the  sinulation  has 
been  completed,  statistical  methods  can  be  employed  to  test 
the  sensitivity  of  the  measures  of  effectiveness  to  the  select- 
ed inputs.   In  addition  a  certain  amount  of  parameter  testing 
can  be  accomplished  by  hand  or  desk  calculations  as  the  model 
is  being  built.   For  the  most  part,  however,  factors  irust  be 
selected  on  the  basis  of  experience,  that  of  both  designer 
and  user,  conditioned  by  what  was  learned  during  the  prelim- 
inary investigation  of  the  problem. 
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An   example  from  .'l9S-1  r^ay  lend  nore  meaning  to  this 
assertion.   Probability  of  ■z'^^inlri'Z   a  non-submarine  contac  c 
and  naxinum  time  to  evaluate  such  contacts  are  included  as 
input  parameterr.   Actually  these  represent  any  delay  exper- 
ienced by  a  helicopter  irhile  conducting  a  sonar  dip;  non-sub- 
marine contacts  are  one  of  raany  reasons  for  such  delays. 
That  delays  occur  is  \Tell   kno"5Tn  to  those  who  hav.^  parlijipated 
in  helicopter  search  operations.   However,  the  existence  or 
siGnlficance  of  variation  in  dip  tines  would  not  necessarily 
be  apparent  to  an  analyst  who  had  not  participated  in  an^.i- 
submarine  operations,  solely  from  a  review  of  the  literature. 

This  introduces  another  question  to  be  considered  oy   the 
model  desicTier.   To  what  extent  may  too  close  association  with 
the  day  to  day  operation  of  the  system  limit  the  porspec tlvc 
of  the  simulation  customer?   For  example,  a  catapult  officer 
may  be  a  valuable  source  of  information  concerning  aircraft 
carrier  air  operations.   However,  this  officer  is  involved 
with  many  factors  which  are  essential  to  air  operations  but 
may  not  be  pertinent  to  the  purpose  or  level  of  arji^recation 
of  a  simulation  of  air  operations.   Gaining  the  maxiraun  in- 
formation from  this  officer,  without  clutterinc  the  simulation 
with  unnecessary  detail  also  based  on  the  man's  experience 
is  a^ain  a  natter  of  Judf^ment. 

The  designer  must  also  appreciate  the  role  of  the  decision 
maker  when  considering  inputs.   It  is  the  uc-^er  who  inist  take 
the  final  rasponsibility  for  decisions  based  on  i.he  niiiiula- 
tion  results,  and  hence  these  decisions  uill  b^  influenced  to 
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a  lar^e  extent  'oy  his  confidence  in  the  simulation.   In  the 
above  example  if  exclusion  of  catapult  pressure  fron  the  sim- 
ulation will  result  in  a  complete  lack  of  faith  in  the  simu- 
lation by  the  decision  maker,  then  catapult  pressure  may  be 
a  valid  input  even  if  the  simulation  results  are  insensitive 
to  this  input.   It  may  be  that  only  by  including  this  factor 
and  conducting  a  more  detailed  analysis  of  the  sensitivity  of 
response  to  it  can  the  question  be  resolved  to  either  the 
designer's  or  user's  satisfaction. 

Factor  levels*   In  selecting  factor  levels,  Just  as  uith 
the  factors  themselves,  the  desif:;ner  will  usually  want  to  in- 
clude a  certain  aiiount  of  generality  in  the  model.   The  factor 
levels  chosen  at  this  time  vrill  in  some  cases  be  the  particu- 
lar levels  used  for  the  analysis  later  on.   In  most  cases, 
hovrever,  the  designer  will  be  interested  in  a  range  from  which 
specific  levels  may  be  chosen  at  the  appropriate  time.   In  this 
phase  of  simulation  design,  planning  the  experiment  around 
specific  questions  can  be  particularly  productive.   If  only 
certain  selected  ranges  of  factor  levels  have  a  bearing  on 
the  study  objectives,  time  spent  determining  the  entire  range 
of  possible  values  might  be  essentially  wasted.   Restricting 
the  levels  to  interesting  ranges  can  also  result  in  many  ef- 
ficiencies in  the  subsequent  computer  program.   As  the  range 
of  factor  levels  is  increased  so  is  the  computer  sLoragc 
space  necessary  to  contain  the  associated  numerical  values. 

Determining  factor  levels  will  also  give  rise  to  the 
need  for  determining  certain  distributions  of  associated 
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input  variables.   It  is  in  this  step  that  the  data  collection 
role  of  statistical  analynis  is  of  particular  interest  to 
the  simulation  desirner. 

Input  parameters  are  directly  related  to  the  particular 
questions  posed  by  the  person  or  activity  requesting  the  sim- 
ulation study.   Consequently  the  customer  is  often  the  best 
source  of  the  range  of  parameter  values.   Input  variables, 
however,  are  more  closely  associated  vrith  the  possible  states 
of  nature.   The  designer  is  free  to  include  or  reject  a  speci- 
fic variable  according  to  its  pertinence  to  the  simulation 
purpose.   Once  he  has  included  a  particular  variable  in  the 
model,  however,   the  analyst  must  determine  the  rang'^  of  values 
of  the  variable  according  to  what  actually  exists  in  the  real 
world.   This  determination  might  involve  a  comparatively 
simple  library  search  or  an  extensive  analysis.   Determination 
of  factor  levels  may  at  times  give  rise  to  analytic  or  simu- 
lation studies  which  are  more  extensive  than  the  original 
simulation. 

It  might  happen  that  the  required  information  cannot  be 
obtained  in  the  time  allowed  or  at  a  cost  consistent  vrith  the 
designer's  operating  budget.   Perhaps  the  information  is  not 
available  at  any  cost,  or  only  a  very  restricted  segment  of 
the  true  range  of  the  variable  can  be  determined  vrith  any  de- 
gree of  certainty.   The  designer  may  leave  this  variable  out 
or  he  may  resort  to  an  educated  guess.   'That  is  done  is  a 
matter  of  Judgment  and  depends  on  the  particular  situation. 
It  is  essential,  however,  that  the  simulation  roflect  the 
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nature  and  degree  of  uncertainty  relative  to  this  input. 
As  vrith  any  analysis  the  results  have  validity  only  with 
respect  to  the  assumptions  and  input  data  from  which  they  come. 
4.5  Corabinin.-^  factors  and  a3sur.r)tions  in  the  model. 

As  an  illustration  of  how  the  assumptions,  factors  and 
factor  levels,  and  algorithms  specify  the  model,  we  will  con- 
sider sonar  detection  ranges  in  AHS-1.   Pron  the  preliminary 
study  it  was  determined  that  the  simulation  purpose  dictated 
that  sonar  detection  range  was  an  essential  element  of  the 
helicopter  model.   This  quantity  is  essentially  dependent  on 
nature  and  was  therefore  defined  as  an  input  variable. 

Further  investigation  revealed  that  detection  ranges  de- 
pend on  certain  variables  such  as  figure  of  merit,  water  con- 
ditions and  sea  state.   These  were  tentatively  assigned  as 
factors  and  the  process  of  determining  general  factor  levels 
vras  commenced.   As  noted  above,  assigning  factor  levels  is 
closely  related  to  determining  the  range  of  input  variables, 
and  data  collection. 

The  first  step  in  assigning  general  factor  levels  was 
to  determine  the  functional  relationship  between  the  original 
input  variables,  sonar  detection  range,  and  the  tentatively 
selected  factors.   Investigation  in  this  area  revealed  the 
stochastic  nature  of  figure  of  merit  and  resulted  in  a 
modification  of  the  selected  inputs.   That  is  the  distribution 
parameters,  average  figure  of  merit  and  variance,  replaced 
figure  of  merit  as  inputs  to  the  simulation.   This  of  course 
involved  certain  assumptions  which  have  already  been  mentioned. 
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On  the  basis  of  the  above  research  and  assumptions  the 
valid  ranjre  of  factor  levels  was  determined.   It  was  then 
only  left  to  discover  al^orithris,  by  which  the  specific 
values  of  the  original  input  variable  could  be  computed,  and 
combine  these  vrith  other  al^^orithms  which  would  ciove  the 
helicopter  model  over  the  playing  area.   It  was  of  course 
necessary  that  these  algorithms  not  only  be  mathematically 
correct  but  that  they  be  applicable  over  the  range  of  both 
factors  and  the  original  input  variable.   The  end  product 
of  this  process  was  a  point  which  could  be  moved  about  a 
square  grid,  and  an  associated  detection  range.   This  is  one 
model  of  a  helicoDter. 
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CHAPTER  V 
BIPLEIiENTING  THE  MODZZL 

5.1   Adding  l0:~lc  to  the  model* 

The  model  is  the  heart  of  the  simulation,  but  to  imple- 
ment the  model  it  is  necessary  to  move  the  simulated  system 
components  throu^rh  time  as  vrell  as  space.   The  simulation  logic 
is  the  means  by  irhich  the  model,  a  static  representation  of 
the  system  components,  is  transformed  into  a  dynamic  time 
ordered  simulation  of  the  system.   In  practice  the  design 
of  model  and  loric  may  be  concurrent,  but  conceptually  there 
exists  a  distinct  interface  between  them. 

To  illustrate  this,  we  refer  to  the  helicopters  in  .IHS-l. 
As  stated  earlier,  a  helicopter  is  represented  by  a  point  on 
a  rec  tanciular  '^rid  torether  with  an  associated  sonar  detec- 
tion ranee.   These,  combined  with  a  suitable  algorithm  for 
chan^in'^  the  £rid  position  of  this  point,  complete  the  model 
of  a  helicopter.   The  implementation  of  this  model  in  a  sim- 
ulation of  a  iielicopter  search  operation  requires  some  pro- 
cedure for  causing  certain  actions  to  take  place  at  specified 
times.   Each  helicopter  must  move  to  successive  dip  stations 
at  prescribed  times.   Having  arrived  at  a  dip  station,  a 
search  for  the  submarine  and  a  determination  of  the  outcome 
of  that  search  must  be  made.   To  complicate  this  picture, 
other  helicopters  are  simultaneously  acting;  and  interacting 
with  the  target  submarine. 
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The  determina-':ion  of  vrhich  actions  epd  to  take  place 
and  their  relativr-  order  of  oocurren::^e  could  be  accomplished 
manually,  usin^^  the  above  model.   Similarly,  lo£:lcal  struc- 
tures other  than  the  one  ur.ed   in  AHS-1  could  be  employed  to 
implement  this  saine  model  in  a  computer  simulation. 

The  two  most  frequently  encountered  logical  structures 
are  the  ev'ini.  itore,  employed  in  .\HS-1,  and  the  time- f^ tap.  Both 
of  these  structures  are  essentially  techniques  for  approxi- 
raatint:;  time,  which  is  continuous  in  the  real  world,  by  i 
series  of  discrete  steps.   In  an  actual  operation  certain 
events  take  place  either  simultaneously  or  in  an  overlapping 
time  space,  but  a  computer  can  only  perform  one  operation 
at  a  time.   Consequently  when  an  action  occurs  ia  the  sim- 
ulation, time  must  be  arrested  for  all  units  not  participating 
in  that  action.   Some  adjustment  must  then  be  made  for  inter- 
acting: events  which  occur  simultaneously.   This  necessi- 
tates some  predetermined  logical  order  in  which  actions  and 
participating  units  are  to  be  considered  during  the  play. 

Time-step.   In  a  time-step  simulation  the  play  takes 
place  as  a  series  of  discrete  time  increments.   The  length  of 
each  interval  is  called  the  time  step  and  may  be  thought  of  as 
representing  the  length  of  one  rame  move. 

During  each  time  ':5tep  the  positions  and  capabilities 
of  all  uaits  are  noted  and  a  list  of  possible  actions  is 
consulted.   For  each  entry  in  the  list  of  possible  happenings, 
all  units  which  could  be  affected  by  this  action  are  then  con- 
sidered in  turn,   'fith  respect  to  each  of  these  units, 
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calculations  are  perfomed  to  see  if  the  action  in  question 
did  in  f.?.ct  oocur  and  to  what  extent  the  unit  vras  affected. 
After  every  ac  lion  in  i.he  list  has  been  considered,  the 
positions  and  capabilities  of  all  participants  are  reassessed 
and  the  play  moves  into  the  next  time  interval.   The  end  of  a 
play  may  be  sifj;nalled  by  the  occurrence  of  some  predetermined 
action  or  at  the  completion  of  a  specified  number  of  time 
steps. 

It  is  characteristic  of  this  method  that  all  events 
vrhich  occur  within  the  same  time  step  are  considered  to  have 
taken  place  either  simultaneously  or  in  the  order  in  which 
they  appear  in  the  list  of  possible  actions.   For  example, 
two  possible  events  might  be  weapon  detonations;  and  missile 
launches.   If  missile  launches  are  considered  ahead  of  de- 
tonations, then  it  would  be  possible  for  a  missile  to  launch 
vrhen  in  fact  the  missile  site  had  been  destroyed  by  a  vreapon 
detonation  prior  to  scheduled  launch  time.   The  probability 
of  this  happening;  can  be  decreased  by  dccreasin::  the  length 
of  the  time  step,  but  at  the  cost  of  increased  computer  run- 
ning time. 

Event  Store.  An.   event  store  simulation  is  played  as  a 
sequence  of  events  rather  than  discrete  time  intervals.   In 
contrast  vfith  the  time-step  method,  these  events  are  not  con- 
sidered until  they  are  scheduled  to  occur.   The  scheduling 
of  future  events  is  accomplished  by  entering  elements,  com- 
monly  referred  to  as  event  words,  in  a  table  called  the  event 
store ♦   ^ach  event  word  consists  of  the  type  of  event,  the 
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tine  the  event  is  to  take  place  and  any  other  information 
which  is  required  for  the  execution  of  the  event.   In  the 
sample  program,  .UiS-1,  the  identity  of  the  unit  involved  in 
the  event  is  part  of  this  word.   In  many  simulations,  events 
once  entered  in  the  list  may  be  cancelled  as  a  result  of 
subsequent  action,  necessitatinn;  the  inclusion  of  a  validation 
signal  in  each  event  word.   Events  which  have  been  invalidated 
may  then  be  discarded  when  they  reach  the  top  of  the  event 
store. 

It  is  characteristic  of  the  event  store  logic  that 
sections  of  the  program  representing  events  are  essentially 
independent  of  each  other.   Interactions  among  these  separate 
components  are  caused  to  take  place  by  two  events  usually 
referred  to  as  Take  Next  Event  (TNE)  and  Store  New  Event  (SNE) 
Store  New  Event  enters  event  words  representing  future  events, 
in  the  event  store  in  chronological  order,  and  Take  Nex: 
Event  causes  events  to  be  executed  at  the  appropriate  time. 

At  the  beginning  of  a  play  certain  events  which  can  be 
predicted  from  the  program  logic  and  the  input  data  are  enter- 
ed in  the  event  store.  Control  of  the  program  then  passes  to 
TNE  which  notes  the  event  type  of  the  earliest  valid  event 
word  and  calls  the  appropriate  event  subroutine  into  action. 
This  event  subroutine  carries  out  its  preassigned  function 
according  to  the  information  contained  in  the  event  word  and 
the  input  data.   As  a  result  of  this  action  other  events  may 
arise,  in  which  case  SNE  is  called  upon  to  store  the  necessary 
event  vrord.   At  the  completion  of  each  event,  control  is  re- 
turned to  TNE  and  the  cycle  repeated.   The  pln.y  terminates 
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when  some  particular  action  occurs,  such  as  a  submarine  de- 
tection in  .■\HS-1 ,  or  irhen  there  are  no  more  events  in  the 
event  store  to  be  processed. 

Time-step  vernus  TiVent  store*   Efficiency  of  both  pro- 
Gramminc  and  computer  operation  are  usually  prime  considera- 
tions in  choosing  the  lo3ic  for  a  simulation.   The  event  store 
method  is  somewhat  more  compatible  with  the  functional  sub- 
divisions of  the  model  previously  mentioned.   Each  event 
subroutine  usually  corresponds  to  some  particular  function 
of  a  system  component.   On  the  other  hand,  the  programmer  may 
have  difficulty  understanding  the  chronolocical  sequence  of 
events  at  first.  [22] 

Since  both  structures  are  methods  of  approximating  real 
continuous  time  by  discrete  steps,  the  number  and  hence  the 
length  of  each  step  is  closely  correlated  with  machine  run- 
ning time.   With  the  time-step  procedure  one  time  step  may 
be  considered  as  a  single  approximating  step.   Every  possible 
action  and  every  unit  is  considered  once  each  time  increment, 
so  as  the  number  of  time  steps  is  increased,  machine  running 
tine  becomes  greater.   Decreasing  the  number  of  operations 
by  lengthening  the  time  step  results  in  greater  time  aggre- 
gation which  is  usually  undesirable.   In  general  the  time-step 
method  results  in  greater  computer  running  time  when  the  exact 
order  of  occurrence  of  closely  spaced  events  is  a  critical 
consideration. 

In  the  event  store  method  each  event  may  be  considered 
as  one  approximation.   With  this  structure,  actions  begin  in 
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the  proper  time  sequence  but  in  the  machine   only  one  event 
can  occur  at  any  one  tine.  Consequently  interactions  among 
events  which  occur  in  overlapping  time  intervals  may  be  lost. 
An   ercample  of  this  is  the  relation  between  submarine  maneuvers 
and  submarine  detection  events  in  AHS-1.   A  detection  event 
consists  of  one  or  more  sonar  searches.   If  the  submarine  is 
due  to  maneuver  vrhile  a  detection  event  is  in  pro^recs,  then 
the  correct  submarine  position  will  not  be  available  from  the 
time  the  maneuver  is  scheduled  until  the  completion  of  the 
event  being  processed.   The  effect  of  this  could  be  minimized 
by  defining  each  individual  search  as  one  event.   This  would 
essentially  increase  the  number  of  events  and  would  result  in 
increased  machine  running  time.   The  procedure  that  is  used 
in  this  case  is  to  interrupt  any  detection  event  between 
searches  if  a  target  maneuver  is  to  occur.   This  is  determ- 
ined by  asking  if  the  submarine  maneuvered  at  the  end  of 
each  search.   The  submarine  is  then  maneuvered  and  the  detec- 
tion event  resumed  at  the  time  it  was  interrupted. 

The  above  is  an  example  of  superimposing  a  time-step  on 
one  portion  of  the  event  store  logic.   Similarly  an  event 
store  may  be  incorporated  into  a  basic  time-step  structure. 
This  technique  is  employed  in  STAGE,  a  global  atomic  exchange 
model  built  by  Technical  Operations  Incorporated  under  contract 
to  the  U.  S.  Air  Force.   The  STAG^  simulator,  which  is  played 
in  fifteen  minute  time  intervals,  uses  the  output  of  a  pre- 
processor for  its  input.   Part  of  the  STAGi:  preprccessor  out- 
put is  a  list  of  events  vrhich  have  positive  prooability  of 
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occuring  during  each  time  interval.   This  list  is  read  from 
magnetic  tape  into  the  main  simulation  program  as  it  is  need- 
ed and  essentially  represents  an  external  event  store.   The 
simulation  itself  then  "Steps  through  each  time  period,  con- 
sidering the  various  events  scheduled  to  occur  in  each  time 
period  and  determining  which  events  actually  occur."l  23J 

Flow  charts  of  the  lop:ic .   Block  flow  diagrams  can  be 
very  useful  in  organizing  the  logic  as  they  were  in  delineating 
the  functional  subdivisions  of  the  model.   Plow  charts  also 
provide  a  convenient  vehicle  for  communication  between  sim- 
ulation designer  and  programmer. 

Plow  charts  at  this  stage  of  the  design  will  usually  be 
constructed  on  at  least  two  levels  of  detail.   General  flow 
charts  consisting  of  boxes  enclosing  plain  English  statements 
are  one  means  of  describing  the  model  and  the  logical  way 
in  which  model  components  interact.   General  flow  charts, 
along  with  the  documentation  of  the  structural  model  compon- 
ents may  completely  specify  the  model  and  logic.   However, 
more  detailed  flovr  diagrams  must  usually  be  constructed  be- 
fore a  professional  programmer  can  reduce  the  model  and  logic 
to  a  useable  computer  program. 
5«2  Computer  lanp:ua,n:e3 

Having  reduced  the  model  and  logic  to  detailed  flow 
charts  the  designer  is  ready  to  turn  his  creation  over  to  a 
programmer,  or  as  is  sometimes  the  case  to  write  the  computer 
program  himself.   A  computer  program  is  essentially  a  list  of 
instructions  written  in  a  format,  or  language,  which  is  use- 
able by  the  computer.   The  prograun  is  the  vehicle  by  which  a 
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complicated  mathematical  model  is  reduced  to  the  very  limited 
set  of  operations  which  can  be  performed  hy   a  computer. 

A  complete  hierarchy  of  computer  languages  is  available 
and  the  selection  of  one  of  these  depends  on  several  factors. 
Some  of  these  considerations  ares   (1)  simulation  size, (2) 
expected  number  of  runs  of  the  completed  simulation, (3)  pro- 
gramiing  facilities,  (4)  available  time  for  programming  and 
pro::ram  testing,  and  (5)  the  available  computer  facilities. 

Computer  languages  may  be  classified  according  to  re- 
lative programming  ease.   In  general  the  simpler  the  program- 
ming the  more  inefficient  the  program  is  in  its  utilization 
of  operating  time  and  memory  space.   One  possible  classifica- 
tion is  as  follows? 

1.  Binary  code 

2.  Assembly  language 

3.  Compiler  language 

4.  Simulation  oriented  language 

Binary  code.   This  is  the  only  language  the  computer 
actually  "understands"  and  all  other  language  systems  must 
ultimately  reduce  the  programmers  instructions  to  the  binary 
number  system.   Writing  binary  instructions  is  extremely 
time  consuming  since  complicated  operations  must  be  reduced 
by  the  programmer  to  a  sequence  of  basic  operations  ;ihich  the 
computer  itself  can  accomplish.   These  basic  operations  con- 
sist of  storing  and  transferring  information,  addition,  sub- 
traction, multiplication  and  division. 

Assembly  or  syinbolio  lanfTjua.'Tjes.  These  languages,  such 
as  IBM  PAP  and  the  symbolic  language  used  with  the  CDC  1604 
employ  what  is  called  an  assembly  program.   The  programmer 
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writes  assembly  language  Instructions  which  are  then  trans- 
lated into  machine  instructions  by  the  assembly  program.   The 
programmer  might  write  a  statement  such  as  PAD  (R),  which  is 
symbolic  for  Floating  Point  Add  (Location  R).   This  state- 
ment instructs  the  computer  to  note  the  contents  of  a  loca- 
tion previously  designated  as  R  and  add  it  to  what  ever  has 
been  placed  in  a  section  of  the  machine  called  the  A  register. 
This  sum  will  then  be  left  in  the  A  register  and  may  be  stored 
in  another  location  by  additional  Instructions. 

This  is  a  much  simpler  instruction  to  write  than  the 
corresponding  binary  configuration.   In  addition,  the  number 
of  the  location  in  which  R  is  stored  need  not  be  known  to  the 
programmer,  as  the  assembler  will  keep  track  of  R  once  it  has 
been  defined.   VThen  binary  machine  language  is  used  the  program- 
mer must  keep  track  of  every  location  and  its  contents  by 
number. 

Compiler  languages.   The  use  of  languages  such  as  P0RTR.1N 
or  ALGOL  further  simplify  the  programmer's  task.   In  general, 
however,  the  step  from  assembly  to  compiler  language  involves 
a  considerable  loss  in  storage  space  and  flexibility,  and  an 
increase  in  operating  time.   The  FORTRAN  program  is  written 
as  a  sequence  of  arithmetic  statements  such  as  R  =  R+A.   This 
means  take  the  contents  of  the  location  assigned  to  R,  add 
this  to  the  location  assigned  to  A  and  place  the  answer  back 
in  R.   Here  the  format  is  again  floating  point  but  this  is 
signaled  to  the  compiler  by  the  particular  letters  used. 

The  FORTRAN  program,  consisting  of  statements  such  as 
the  one  above  and  many  other  instructions  which  are  a  great 
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deal  more  coniplicated,  must  first  be  compiled.   This  consists 
of  a  compiler  program  translating  the  FORTRAN  statements  in- 
to an  assembly  language  program.   This  assembly  language 
program  is  then  converted  to  machine  instructions.   This  gen- 
erally involves  inefficiencies  both  in  the  utilization  of 
time  and  storage  space.   The  compiler,  lacking  the  ability 
to  reason^  is  not  as  sophisticated  as  a  human  programmer  and 
tends  to  waste  instructions  to  some  extent. 

Compiler  languages  do  have  some  very  Important  advantages 
In  addition  to  ease  of  programming.   One  of  these  is  the  large 
library  of  function  sub-routines  and  programs  which  is  avail- 
able.  For  instance  the  simple  statement  R  =  SINP(THETA)  can 
be  used  to  compute  the  sine  of  an  angle  THSTA  and  store  this 
value  in  the  location  R.   To  accomplish  this  using  an  assembly 
language,  the  programmer  would  have  to  write  the  instructions 
for  computing  the  Taylor's  expansion  of  sine  THETA,  add  the 
terms  out  to  some  desired  accuracy  and  place  the  sum  in  lo- 
cation R. 

Simulation  oriented  languages.  Languages  designed  pri- 
marily for  simulation,  such  as  SIMSCRIPT  or  MILITRAN  are  cur- 
rently being  developed.   24,25]   These  languages,  or  program- 
ming systems,  are  designed  to  further  simplify  programming. 
In  general,  the  program  is  written  as  a  series  of  functions 
comparable  to  the  event  subroutines  in  AHS-1,  rather  than 
arithmetic  statements.   The  source  progrsun  consists  of  phrases 
similar  to  those  commonly  used  in  general  flow  diagrams. 
These  statements  are  then  translated  into  subroutines  written 
In  a  core  language  such  as  FORTRAN. 
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Choosing  a  computer  language o   Probably  the  first  con« 
sideration  In  selecting  a  language  is  size  of  the  resulting 
program.   If  the  problem  will  not  fit  into  an  available 
machine,  running  time  and  programming  considerations  become 
academic.   If  the  simulation  is  large,  however,  the  program 
size  should  be  considered  in  conjunction  with  available  data 
processing  facilities.   In  many  cases  the  effective  machine 
size  can  be  expanded  far  above  actual  size  through  the  use  of 
satellite  equipment  such  as  magnetic  tape  drives.   Parts  of 
the  program  and  input  data  may  be  stored  on  magnetic  tape 
external  to  the  machine.   The  taped  information  is  then 
read  into  the  computer  only  when  it  is  needed. 

If  space  requirements  indicate  that  a  compiler  language 
is  feasible  then  the  balance  between  machine  operating  cost 
and  programming  expense  is  very  important.   If  a  simulation 
is  to  be  run  only  a  few  times,  a  compiler  language  is  prob- 
ably indicated  because  of  the  lower  programming  cost.   Of 
course  if  programming  in  assembly  or  machine  language  would 
delay  the  completion  until  after  the  information  was  no  long- 
er of  interest,  a  compiler  program  might  be  the  only  altern- 
ative no  matter  how  many  simulation  runs  were  anticipated. 

The  availability  of  existing  subroutines  which  may  be 
useable  in  the  program  should  also  be  given  careful  consider- 
ation.  In  selecting  a  language  for  AHS-'l,  the  programmer  was 
oblli^ed  to  choose  between  two  compiler  languages,  PORTRAU 
and  NSLIAG.   This  restriction  arose  because  the  program  was 
intended  as  an  example  and  these  two  languages  are  the  more 
commonly  used  within   the  United  States  Navy.   The  final 
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decision  to  use  FORTRAN  was  greatly  influenced  by  the  avail- 
ability of  a  sonar  range  prediction  program  written  in  FORTRAN. 

One  of  the  most  important  drawbacks  to  FORTRAN  for  sim- 
ulation studies  is  the  inability  to  pack  more  than  one  quan- 
tity into  each  word.   For  example  the  event  store  table  used 
in  AHS-1  consists  of  three  arrays*   In  assembly  language  these 
could  be  combined  into  one  table  by  allotting  only  part  of  a 
word  to  each  of  the  three  items  of  information.   It  is  possible 
to  write  packing  and  unpacking  routines  in  assembly  language 
and  combine  them  with  a  FORTRAN  program.   Such  routines  take 
time  to  program  and  it  might  be  almost  as  convenient  to  write 
the  entire  code  in  an  assembly  language. 

The  various  dialects  of  ALGOL,  such  as  NELIAC  or  JOVIAL, 
do  retain  the  packing  feature  of  an  assembly  language,  but 
at  the  cost  of  an  increase  in  computer  operating  time. 
5.3  Proriram  check-out 

Regardless  of  the  lauiguage  in  which  the  computer  program 
is  written,  it  can  be  expected  to  contain  some  errors.   Care- 
ful planning  and  attention  to  detail  will  go  a  long  way  to- 
ward minimizing  these  errors ^  but  a  computer  code  is  an  ex- 
tremely complicated  structure p  smd  the  specifications  for 
computer  prograims  are  very  exacting.   The  omission  of  one  in- 
struction or  even  one  decimal  point  may  result  in  the  entire 
program  being  rejected  by  the  computer. 

It  is  probably  easier  to  test  the  prograjn  if  it  has  been 
constructed  in  Independent  sectionsj  and  it  was  partly  toward 
this  end  that  the  model  was  divided  Into  functional  components. 
Each  of  these  sections  provides  a  particularly  convenient  unit 

when  tenting  the  program  because  of  the  relation  between  input 

48 


data  and  Individual  functional  components.   If  the  coding  is 
tested  in  blocksj  input  data  which  would  normally  be  provided 
by  other  sections  of  the  program  is  usually  introduced  by  a 
check-out  control  program.   This  may  consist  of  a  few  cards 
containing  the  necessary  data  or  might  incorporate  previously 
tested  subroutines.   In  this  way,  data  which  is  Input  to  the 
subsection  being  tested  can  be  precisely  controlled.   Local- 
izing an  error  in  a  program  section  which  uses  self  generated 
input  data  is  complicated  by  the  fact  that  both  the  program 
logic  and  input  data  are  possible  sources  of  the  error. 

Although  coding  the  program  in  subsections  is  probably 
the  most  significant  single  aid  to  program  testing,  it  is  not 
a  complete  panacea.   Though  each  section  performs  properly 
over  its  designed  range  of  inputs,  this  does  not  insure  that 
all  of  these  components  will  interact  with  one  another  as  the 
designer  intended.   One  method  of  testing  the  completed  pro- 
gram is  to  compare  a  limited  number  of  runs  with  hand  calcu- 
lated results.   This  will  generally  involve  checking  the 
program  at  certain  critical  points.   The  procedure  used  in 
testing  AHS-1  was  to  have  the  values  of  key  variables  printed 
at  certain  points  throughout  the  program.   In  this  way  the 
output  from  one  operation  and  the  input  to  the  next  link  in 
the  chain  were  checked  simultaneously.   These  were  then  com- 
pared with  hand  calculationso 

The  above  procedure  is  a  tedious  and  time  consuming  task 
for  even  a  relatively  uncomplicated  program^  but  it  may  be 
the  only  way  the  simulation  designer  can  be  reasonably  con- 
fident of  the  correctness  of  the  program.   Even  when  such  a 
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testing  procedure  has  been  completed  there  is  generally  no 
way  of  knowing  with  one-hundred  percent  confidence  that  the 
program  is  free  of  errors.   Such  assurance  would  ordinarily 
Involve  testing  every  possible  combination  of  Input  values 9 
a  virtually  Impossible  tasko 
5.4  Auxiliary  program  components. 

The  questions  of  what  to  Include  as  input  and  output  data 
confronted  the  simulation  designer  throughout  the  model  de- 
sign phase.   The  question  of  how  to  get  this  information  in 
and  out  of  the  computer  must  also  be  answered.   With  the  cap- 
acity to  process  large  amounts  of  data  Inherent  in  a  digital 
computer,  data  processing  is  not  a  trivial  consideration. 

The  input  data  format  will  often  be  Influenced  to  a  con^ 
siderable  extent  by  the  size  of  the  program.   Space  limita- 
tions may  necessitate  reading  data  in  only  as  it  can  be  used. 
If  space  is  not  a  problem^  time  may  be  saved  by  reading  as 
much  data  at  one  time  as  is  practicable  and  storing  it  intern- 
ally. 

Space  and  time  considerations  will  affect  the  particular 
form  of  input  data  as  well.  In  a  small  simulations,  data  may 
be  Introduced  in  a  form  which  is  more  convenient  to  the  user 
and  any  necessary  conversions  accomplished  Internally.  Large 
simulations  might  require  considerable  pre-processing  of  data 
to  conserve  room  for  calculations  which  can  be  conducted  only 
in  conjunction  with  the  operating  section  of  the  program. 

Errors  in  the  input  data  are  also  an  important  considera- 
tion.  Data  input  to  AHS-1  is  printed  out  as  a  check  on  what 
values  were  actually  used  by  the  simulations,  but  an  important 
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section  of  the  program  was  omitted  because  of  time  limitations. 
This  omission  is  a  routine  for  checking  the  input  quantities 
to  see  that  they  fall  within  acceptable  limitSc   This  would 
not  eliminate  mistakes  entirely  but  could  reduce  lost  comput- 
ing time  due  to  misplaced  decimal  points  and  other  gross  errors 
in  the  input  data. 

The  output  section  of  the  program  will  be  influenced  to 
a  considerable  extent  by  the  simulation  objectives.   In  addi- 
tion to  selecting  outputs  which  will  yield  the  maximum  informa- 
tion, the  designer  should  specify  an  output  format  which  is 
compatible  with  the  particular  analysis  techniques  to  be 
employed.   For  example  when  space  permits^  computation  of 
sample  means  and  variances  might  save  the  analyst  a  great  deal 
of  time.   Similarly,  percentages  or  ratios  of  numbers  are  often 
more  meaningful  than  the  numbers  themselves. 
5'5   Testing  the  simulation. 

In  designing  the  modely  factors  were  included  because 
they  were  known,  or  In  some  cases  only  suspected,  to  affect 
the  system  In  a  way  which  was  important  to  the  study  purpose. 
Before  the  simulation  is  ready  for  use  as  an  analytic  tool, 
the  designer  should  determine  both  quantitatively  and  qualita- 
tively how  these  factors  affect  simulation  results.   This  will 
require  that  some  testing  program  be  conducted  and  representa- 
tive levels  of  all  or  part  of  the  factors  examined. 

Requirements  of  a  testing  program.  The  testing  program 
may  include  a  comparison  of  simulation  results  with  observa- 
tions of  the  system  itselfy  or  with  experimental  results  ob- 
tained by  other  analytic  methodSc   It  should  Include  some  an- 
alysis of  the  sensitivity  of  measures  of  effectiveness  to 
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changes  in  factor  levels.   Whatever  the  analytic   echniques 
employed,  the  following  appear  to  be  fairly  general  objectives 
of  the  testSo 

(1)  Verify  that  simulation  response  to  well  known  factors 
is  in  general  agreement  with  observation. 

(2)  Determine  the  sensitivity  of  response  to  all  factors- 

(3)  Determine  the  number  of  replications  necessary  to 
attain  certain  levels  of  statistical  validity. 

(4)  Determine  which  factors  interact  and  the  extent  of 
such  interaction. 

Conducting  the  testing  program.   The  number  of  runs  neces- 
sary to  attain  the  desired  accuracy  should  be  one  of  the  first 
considerations  during  simulation  testing.   In  testing  the  simu- 
lation, as  well  as  in  conducting  subsequent  experiments,  the 
analyst  is  usually  guided  by  two  major  considerations,  how  to 
gain  as  much  useful  information  as  possible  from  the  simulation, 
and  still  minimize  the  expenditure  of  time  and  resourcesc   Test- 
ing should  also  provide  some  general  guidelines  as  to  the  number 
of  replications  necessary  to  achieve  certain  levels  of  confidence 
in  the  simulation  results.   These  may  then  be  applied  during  both 
employment  and  testing  of  the  simulation.   Testing  the  simula- 
tion using  too  few  runs  may  essentially  invalidate  the  results 
of  the  tests. 

When  comparing  simulation  results  with  observations  of  the 
systemj  it  is  well  to  consider  the  source  of  such  observations. 
In  analyzing  military  operations  the  results  of  peace  time  ex- 
ercises may  be  the  nearest  the  analyst  can  attain  to  an  actual 
observation.   Field  maneuvers  and  fleet  exercises  may  at  times 
approach  reality,  but  the  fact  remains  that  they  are  another 
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form  of  simulatlorio   As  with  real  battles,  an  exercise  may  be 
conducted  only  a  few  times.   Consequently  the  analyst  must  con- 
sider the  dangers  associated  with  drawing  inferences  from  small 
samples.   The  generality  of  the  model  must  also  be  considered 
when  C(  mparing  simulation  response  with  observations  from  the 
real  world.   Agreement  between  simulation  results  and  one  parti- 
cular observation  may  reflect  that  the  designer  has  only  simu- 
lated that  particular  observationo 

During  the  Initial  testing  it  may  be  noted  that  certain 
combinations  of  factor  levels  tend  to  saturate  the  simulation 
output.   That  is,  some  factors  because  of  their  greater  effect 
on  the  simulation  response,  may  tend  to  completely  determine 
the  results  when  assigned  values  near  the  extremes  of  their 
ranges.   Knowledge  of  these  saturation  levels  should  be  use- 
ful whea  conducting  sensitivity  studies  on  the  remaining  input 
parameters  and  variables. 

Since  the  simulation  was  designed  as  a  factorial  experi- 
ment from  the  start  It  is  expected  that  this  method  of  analysis 
will  play  an  important  part  in  the  sensitivity  testing.   Factor- 
ial experiments  are  particularly  adaptable  to  sensitivity  anal- 
ysis as  well  as  to  the  determination  of  factor  interactions.  Il0,26j 

R.  J.  Matteis  and  W.  C.  Su.iler  have  described  a  testing 
program  employed  in  sensitivity  testing  of  the  Carmonette  Model. 
l26|  The  testing  program  was  conducted  in  two  phases?   The  first 
phase  consisted  of  examining  the  more  important  factors  and  gen- 
erating some  estimate  of  effects  and  variances.   In  this  phase, 
factors  which  exhibited  saturation  effects  were  treated 
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individually.   The  remaining  factors  were  then  examined  in  a 
complete  factorial  design.,   In  the  second  phase,  knowledge 
gained  during  the  initial  testing  was  used  in  designing  a  frac- 
tional factorial  experiment  to  test  the  overall  model  for  inter- 
actions. 
5.6   Documentation. 

Documentation  can  very  easily  be  the  defining  line  between 
a  partially  completed  computer  simulation  and  a  useful  analytic 
tool.  The  importance  of  assumptions  has  been  emphasized  through- 
out the  preceding  discussion  of  model  design.  Unless  these  as- 
sumptions are  written  down,  however,  they  may  not  be  considered 
during  the  analysis.  This  is  further  complicated  by  changes  in 
personnel  or  by  physical  separation  of  designer  and  ultimate 
user. 

Properly  documenting  program  logic  and  computer  code  is  no 
less  important  than  completely  describing  the  model.   The  real 
world  is  dynamic  and  the  simulation  must  often  reflect  change. 
Attempting  to  modify  a  program  without  complete  flow  charts  and 
a  dictionary  of  variable  names,  as  well  as  the  major  model  as- 
sumptions, can  be  a  formidable  task.   In  some  cases  it  may  be 
simpler  to  rewrite  the  coding  than  to  modify  an  existing  program. 

Standards  of  format  and  adequate  content  vary  with  each  or- 
ganization engaged  in  simulation  design^,  and  standardization 
among  simulation  groups  is  virtually  non-existant .   The  descrip- 
tion of  AHS-1  reflects  the  author's  own  idea  of  adequate  docu- 
mentation, with  one  exception;  the  simulation  has  not  been 
thoroughly  tested  and  no  mention  is  made  of  tests  results.   In 
summary,  the  following  list  is  considered  to  contain  the  minimum 
requirements  for  documentation  of  a  computer  simulation. 
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(1)  A  statement  of  all  major  assumptions  pertaining  to 
the  model  with  appropriate  references,  and  a  descrip- 
tion of  the  model. 

(2)  A  description  of  the  logic.   This  should  include  gen- 
eral flow  chart So 

(3)  Complete  rules  for  input  data  preparation,  and  the 
allowable  range  of  input  values* 

(4)  Statement  of  measures  of  effectiveness  and  a  descrip- 
tion of  the  output  datao 

(5)  A  list  of  all  variable  names  used  in  the  program  with 
their  definitions. 

(6)  A  description  of  the  statistical  test  procedures  used 
and  the  results  of  these  tests. 
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CHAPTER  VI 
AHS-1   A  COMPUTER  SIMULATION 

AHS-1  is  an  event  store  computer  simulation  of  an  antl° 
submarine  warfare  helicopter  small  area  search^   In  accordance 
with  the  definitions  of  Chapter  I^  it  may  properly  be  described 
as  an  analytic  computer  war  game.   A  maximum  of  ten  dipping 
aonar  equipped  helicopters  and  one  evading  submarine  are  the 
principles  in  the  gsune. 

The  simulation  is  intended  as  a  tool  for  the  statistical 
analysis  of  the  comparative  effectiveness  of  different  helicopter 
search  tactics  in  similar  tactical  and  physical  environments. 
The  model  for  AHS-1  and  the  underlying  assumptions  concerning 
sonar  parameters  and  submarine  capabilities  are  described  in 
detail  in  the  succeeding  sections. 
6.1   Tactical  situation  and  play  of  the  game« 

A  play  of  the  game  begins  with  the  submarine  at  a  known 
datum  and  a  flight  of  sonar  equipped  helicopters  enroute  to  or 
in  the  vicinity  of  the  datum.   The  submarine  leaves  the  surface 
at  game  time  zero  and  departs  the  datum  on  a  course  and  speed 
which  is  unknown  to  the  helicopters.   The  helicopters  arrive 
at  datum  at  some  time  determined  by  the  game  user  and  proceed 
to  dip  stations  in  accordance  with  the  assigned  search  plan. 

The  submarine  leaves  datum  on  a  course y  speed  and  depth 
selected  in  one  of  four  ways  according  to  the  game  inputs. 
Each  of  these  operating  modes  reflects  a  different  degree  of 
submarine  randomness,  ranging  from  a  completely  predetermined 
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track  to  random  selection  of  course y  depth  and  speed  within 
restrictions  imposed  by  the  game  planner* 

The  helicopters  proceed  to  their  designated  dip  stations 
at  datiun  time  plus  time  late  and  commence  the  search.   A.t  each 
dip  station  as  many  as  five  sonar  sweeps  at  various  depths  can 
be  executed.   Helicopters  are  assumed  to  be  equipped  witb  scan- 
ning type  variable  depth  active  sonar.   The  time  to  complete 
each  sonar  dip  is  determined  by  input  time  delays  representing 
the  followings   (1)   The  time  initially  required  to  establish 
a  hover  and  to  lower  the  transducer,  (2)  the  times  necessary 
to  change  the  transducer  depth  between  sweeps,  (3)  the  elapsed 
time  from  completion  of  the  last  sonar  sweep  until  the  helicopter 
has  transitioned  from  hovering  to  forward  flight j  and  (4)  time 
spent  in  the  prosecution  of  non-submarine  contacts  if  applic- 
able.  The  dip  cycle  time  is  the  sum  of  dip  time  and  time  en- 
route  between  dips.   Except  for  delays  resulting  from  non-sub- 
marine contactSj,  dip  time  is  the  same  for  each  helicopter.   All 
helicopters  fly  at  the  same  airspeed  between  dips. 

At  the  beginning  of  each  sonar  sweep  the  submarine  position 
is  determined  and  a  target  range  computed.   The  detection  area 
is  doughnut  shaped,  centered  about  the  helicopter" s  positiono 
Detection  will  occur  if  the  submarine  is  within  the  annulus 
determined  by  the  maximum  detection  circle  and  a  smaller  circle 
within  which  detection  cannot  occur  because  the  helicopter  is  too 
close  to  the  target.   This  minimum  detection  range  is  not  de- 
pendent on  sonar  conditions  and  may  be  input  as  zero  if  desired. 
The  maximum  detection  raoige  is  computed  for  each  sweep  and  each 
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helicopter  at  the  beginning  of  the  play  and  recomputed  whenever 
the  submarine  changes  depth.   Detection  ranges  are  based  on  a 
randomly  selected  figure  of  merit  for  each  unit  and  the  prevail- 
ing  sonar  conditions.   Once  selected^  individual  helicopter 
figure  of  merits  are  constant  throughout  each  play© 

The  scheme  for  computing  detection  ranges  was  designed  by 
R.  L.  Klinlcner  of  the  Applied  Physics  Laboratory  [17J  and  is 
considered  by  the  author  of  this  thesis  to  be  an  important  im- 
provement over  the  deterministic  "cookie-cutter"  detection  range 
employed  by  many  sonar  simulations-   The  actual  detection  range 
for  each  helicopter  depends  on  sonar  frequency,,  sea  state »  layer 
depth,  temperature  in  the  layer,  target  depth,  transducer  depth, 
and  the  quality  of  the  sonar  operators-equipment  comblnationo 
The  latter  is  assumed  to  vary  among  helicopters  in  a  prescribed 
random  manner . 

The  play  proceeds  until  all  helicopters  have  completed 
their  last  dip  or  a  detection  occurs.   At  the  completion  of  the 
play  the  running  tally  of  detections  Is  brought  up  to  date  and 
a  new  play  commenced o   This  sequence  continues  until  the  desired 
number  of  plays  have  been  completed.   Input  parameters  and  vari= 
ables  may  then  be  changed  and  another  sequence  repeated «   A 

» 

series  cf  plays  with  any  given  set  of  input  data  will  be  refer- 
red to  as  a  game  run» 
6.2  Submarine  movement q 

The  target  submarine  maneuvers  in  a  three  dimensional  play^ 
ing  area.   Horizontal  movement  Is  relative  to  an  X^Y  coordinate 
system,  with  yards  the  basic  unit  of  distance.   Depths  are 
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measured  in  feet  to  the  nearest  foot.  Courses  are  to  the  near- 
est degree  and  speeds  to  the  nearest  knot.   Position  and  speed 
vectors  are  computed  using  a  standard  polar  coordinate  system^ 
howeveTj,  courses  and  speeds  are  adjusted  so  that  the  Y  axis  cor- 
responds to  north  and  all  courses  are  measured  clockwise  from  ito 

Pour  different  modes  of  submarine  operation  are  available  to 
the  user  of  AHS-1  as  follows? 

Option  (1)  Submarine  track  is  predetermined  by  the  usero 
This  mode  might  be  used  when  studying  the  effect  of  evasive  sub- 
marine maneuvers  on  a  particular  helicopter  search  plan*   If  the 
probability  of  delay  due  to  non-submarine  contacts  (see  Section 
6o4)  is  input  as  zerOj,  the  location  of  the  helicopters  will  be 
known  at  all  times  to  the  game  plannero   Any  level  of  intelli- 
gence concerning  helicopter  movement  may  then  be  attributed  to 
the  submarine.   Similarly ^,  dip  times  may  be  made  essentially 
random  by  assigning  positive  probability  to  non-submarine  con- 
tacts.  In  this  way  the  submarine  may  be  placed  in  the  situation 
of  knowing  where  the  helicopters  are  at  any  one  time  but  not 
knowing  when  they  will  move,  or  where  they  will  Jump   to. 

Option  (2).  Speed,  depth  and  the  bearings  of  each  track  leg 
relative  to  the  first  course  are  predetermined  as  in  Option  (1). 
However^  at  the  beginning  of  each  play  a  pseudo-random  number 
is  selected  from  the  interval  1 0^,360)  degrees  and  added  to  each 
course.   Election  of  this  mode  is  equivalent  to  assuming  a  com- 
plete lack  of  knowledges,  by  the  submarine p  concerning  helicopter 
movement.   As  opposed  to  a  completely  random  submarine  (Option?)^ 
however^  the  capability  of  changing  course »  depth  and/or  speed 
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during  the  play  is  retained. 

Randomizing  the  submarine  track  serves  to  preclude  accident- 
al biasing  of  game  results  by  the  usere   This  bias  may  be  in- 
troduced by  a  particularly  fortuitous  selection  of  the  sub- 
marine track  with  respect  to  the  helicopter  positionso  X   sit- 
uation particularly  suited  to  this  mode  is  that  of  a  submarine 
departing  datum  on  a  straight  course  but  decreasing  speed  at  the 
end  of  preset  time  intervals.   This  tactic  might  be  used  by  a 
submarine  commander  who  desires  to  open  the  range  to  datura  as 
expenditiously  as  possible  without  exhausting  the  submarine 
batteries. 

Option  (3)*  Submarine  course ^,  speed  and  depth  are  uniformly 
distributed  random  variables.   Course  is  distributed  between  zero 
and  360  degrees.   Depth  and  speed  range  between  upper  and  lower 
limits  chosen  by  the  game  user.   Course^,  speedc  and  depth  are  de- 
termined by  generating  three  pseudo-random  numbers  at  the  begin- 
ning of  each  play,  and  remain  the  same  until  the  termination  of 
that  play. 

Option  (A),  Course  and  speed  are  chosen  randomly  as  in 
Option  (3)  but  depth  is  uniquely  determined  by  the  chosen  speed 
and  the  game  inputs.   After  selecting  a  speed,  the  minimum  depth 
is  taken  as  the  shallowest  depth  at  which  the  submarine  can  oper- 
ate without  cavitating.   It  should  be  noted  here  that  cavitation 
does  not  directly  affect  the  possibility  of  detection  in  the 
game,  as  the  helicopters  do  not  conduct  passive  sonar  operations. 
This  submarine  mode  is  introduced  as  a  means  of  restricting  the 
target  to  depths  which  are  realistic  for  the  speed  used.   Target 
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depth  is  a  parameter  which  is  considered  in  determining  active 
sonar  range. 

6.3  Helicopter  movement. 

Helicopter  movement  consists  of  Jumps  betweea  sonar  dip 
stations.   All  helicopters  Jump  at  the  same  speed,  and  as  the 
game  is  now  programmed  dip  stations  are  deterministic.   Provisions 
have  been  made  in  the  model  and  accompanying  computer  code  for 
the  introduction  of  random  bearing  and  course  errors.   However, 
determination  of  the  distribution  of  such  errors  would  require 
the  services  of  fleet  units  not  available  to  the  writer.   With 
the  increasing  sophistication  of  helicopter  navigation  equipment 
the  omission  of  these  errors  is  not  considered  to  be  critically 
detrimental  to  the  purpose  of  this  simulation.   The  provision 
for  including  navigation  errors  is  intended  primarily  as  an  area 
for  further  study  by  the  interested  reader. 

Helicopter  movement  over  the  playing  area  is  unrestricted 
in  azimuth,  and  dip  stations  are  determined  by  the  game  user. 
The  game  is  primarily  intended  to  simulate  close  area  search 
plans,  and  helicopter  dip  stations  are  computed  relative  to 
datum.   VThen  simulating  screening  or  other  support  operations 
the  datum  coordinates  may  represent  an  aircraft  carrier  or  other 
helicopter  take-off  point.   The  first  submarine  maneuver  may 
then  be  used  to  move  the  submarine  from  datum  to  the  desired 
starting  points   This  maneuver  can  be  executed  at  a  high  enough 
speed  so  that  it  will  be  essentially  instantaneous. 

6.4  Submarine  detection. 

Each  helicopter  has  the  capability  of  making  up  to  five  sonar 
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sweeps  at  each  dip  station.   Transducer  depth  and  the  time  re- 
quired to  adjust  transducer  depth  and  conduct  a  search  are  in- 
put by  the  user  and  may  be  different  for  each  sweep.   The  number 
of  sweeps,  sensor  depth  and  time  to  complete  corresponding  sonar 
sweeps  are  the  same  for  all  helicopters. 

No  provision  is  made  for  submarine  motion  during  a  sonar 
sweep.   This  is  of  very  little  consequence  if  scanning  type  sonar 
is  simulated,  but  could  have  some  effect  on  the  outcome  if  search- 
light sonar  is  being  used.   In  the  case  of  searchlight  type  sonar, 
the  time  to  train  the  transducer  may  be  included  in  each  sonar 
sweep  time;   however,  range  to  the  submarine  will  only  be  deter- 
mined  at  the  beginning  of  each  sweep,  regardless  of  the  time  re- 
quired to  step  train  the  transducer  through  a  complete  sweep. 

At  the  beginning  of  a  sonar  sweep  the  range  to  the  submarine 
is  computed  according  to 


Rs  =/(Xj^  -  X^)2  *  (Y^  .   Y^)^  (6.1) 


Where   X^^  and  Y,   are  the  dip  coordinates  of  the  helicopter 


» 


and  X   and  Y   are  the  submarine  position  coordinates.   Detec- 
s       s 

tlon  occurs  if   SR  -  Rs  ^  Rd  where  Rd  is  the  detection  range 
for  the  particular  helicopter  at  the  appropriate  transducer 
depth,  Rs  is  the  target  range,  and   SR  is  the  minimum  range  at 
which  detection  can  occur.   SR  may  be  input  as  zero  if  desired. 
The  course y  speed  and  depth  of  the  submarine  at  the  begin- 
ning of  a  dip  are  used  in  computing  target  range.   If  the  target 
is  scheduled  to  maneuver  during  any  sweep  this  information  is 
noted  and  at  the  completion  of  the  sweep  in  question  the  dip  is 
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discontinued.   The  submarine  maneuver  is  then  executed  and  the 
dip  completed  using  the  correct  submarine  track  information » 

In  addition  to  sonar  sweep  times  and  delays  representing 
the  time  required  to  lower  and  retrieve  the  sonar  transducer, 
one  other  factor  contributes  to  sonar  dip  time.   This  is  a  random 
variable  representing  time  engaged  in  pursuing  non-submarine  con- 
tacts. 

For  the  purposes  of  this  game  the  time  required  to  classify 
actual  submarine  targets  is  not  considered  pertinent.   If  a  sub- 
marine has  been  detected,  the  search  plan  has  been  effective  and 
classifying  the  target  is  another  problem.   It  is  conceivable, 
however,  that  the  effectiveness  of  different  search  plans  may  be 
more  or  less  sensitive  to  variation  in  dip  times  among  helicopters 
One  of  the  primary  factors  contributing  to  differences  in  dip 
time  is  the  classification  of  false  contacts.   In  the  model  it 
has  been  assumed  that  non-submarine  targets  such  as  fish  or 
whales  are  uniformly  distributed  throughout  the  ocean.   This  im- 
plies that  one  sonar  operator  has  about  the  same  chance  of  con- 
tacting such  a  target  as  any  other  sonarman.   Once  a  non-submarine 
contact  has  been  generated,  the  time  required  to  classify  it 
as  non- submarine  will  vary  among  sonar  operators.   These  two 
assumptions  have  been  incorporated  into  the  model  in  the  follow- 
ing way. 

The  actual  delay  for  each  helicopter  is  dependent  on  two 
input  parameters,  the  probability  of  obtaining  a  false  contact 
(P^  )  and  the  maximum  time  any  non-submarine  contact  will  be 
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prosecuted  (Tmax).   A  random  number  from  the  interval  |^0,  ij  is 
first  generated  to  determine  whether  the  helicopter  gained  a 
false  contact  durins  the  dip  in  question.   If  this  random  num- 
ber is  less  than  or  equal  to   P^   a  non-submarine  contact  was 
acquired  and  a  delay  time  must  be  added  to  dip  time.   The  length 
of  time  to  be  added  is  determined  by 

T  =  p^   .   Tmax  (6.2) 

^fc 

WTiere   T  is  the  actual  delay  and  RN  is  the  previously  generated 

random  number. 

6.5   Sonar  detection  ranges. 

Detection  ranges  for  each  helicopter  are  computed  at  the 
beginning  of  every  play  and  whenever  the  submarine  changes  depth 
during  the  play.   As  detection  range  depends  on  transducer  depth, 
every  helicopter  will  have  associated  with  it  as  many  different 
detection  ranges  as  the  number  of  sonar  sweeps  per  dip.   The 
actual  computation  of  these  detection  ranges  is  a  problem  in 
acoustics  rather  than  simulation  and  will  not  be  considered  in 
detail.   The  interested  reader  is  referred  to  Fundamentals  of 
50NAR.  by  J.  W.  Horton  [l9J  for  further  information  on  this  sub- 
ject. 

The  basic  equation  for  active  sonar  isj 

EE  =  M  +  TS  -  2PL  (6.3) 

where   EE  is  the  echo  excess  (signal  level  relative  to  that  re- 
quired for  a  50%   probability  of  detection) 

PM  =  Sonar  figure  of  merit 
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TS  =  Target  strength 
PL  =  One-way  propagation  loss 
All  quantities  in  equation  (6.3)  are  expressed  in  decibels. 
Figure  of  merit  and  target  strength  are  input  variables.  Figure 
of  merit  is  a  measure  of  the  quality  of  the  sonar  operator- 
equipment  combination  independent  of  water  conditions.   Target 
strength  is  a  characteristic  of  the  submarine  size  and  shape,  the 
external  surface  of  the  submarine  hull,  and  target  aspect.   An 
average  value  of  target  strength  is  used  in  the  model.  An   echo 
excess  of  zero  corresponds  to  a  50^  probability  of  detection 
and  this  value  is  used  in  the  computation. 

In  computing  sonar  ranges  it  is  assumed  that  sonar  figure 
of  merit  varies  among  helicopters  according  to  some  known  prob- 
ability distribution.   It  is  also  assumed  that  figure  of  merit 
does  not  vary  appreciably  between  dips  for  the  same  helicopter. 
Therefore  the  figure  of  merit  is  computed  for  each  helicopter 
only  once  for  each  play  of  the  game.   In  order  that  this  thesis 
be  unclassified  no  attempt  has  been  made  to  duplicate  the  actual 
distribution  of  helicopter  figure  of  merit.   For  purposes  of 
AHS-1  a  Normal  (gaussian)  distribution  is  assumed  with  mean  and 
variance  as  input  parameters. 

Rewriting  equation  (6.3)  with  echo  excess  as  zero: 

PL  =  MpM  +  TS)  (6.4) 

Propation  loss  (PL)  combined  with  transducer  depth  and  the  re- 
maining sonar  parameters,  is  used  by  the  routine  which  computes 
detection  ranges  for  each  individual  helicopter.   The  other 
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input  variables  entering  into  detection  range  are  target  depth, 
sonar  frequency,  layer  depth,  temperature  in  the  layer  and  sea 
state . 
6.6   Description  of  program  subroutines* 

The  following  brief  descriptions  of  the  individual  sub- 
routines are  intended  to  show  how  the  model  described  above 
has  been  implemented.   A  secondary  purpose  of  these  descriptions 
is  to  indicate  auxiliary  functions  performed  by  each  event  sub- 
routine.  Rules  for  input  data  card  preparation  are  included 
under  Subroutine  PRINT.   Plow  charts  showing  the  logical  struc- 
ture of  each  subroutine  are  included  at  the  end  of  the  individual 
descriptions.   the  Following  flow  chart  symbols  are  used  through- 
out. 

Flow  chart  symbols.   The  meaning  of  each  symbol  is  determined 
by  shape  and  letters  within  the  symbol.   Differences  in  size  re- 
flect only  space  consideration. 


ENTER 


-  Beginning  point  in  a  subroutine. 


JL 


-  Computation  to  be  carried  out  within 
the  subroutine, 


-  Arrows  show  direction  of  flow, 


Computation  to  be  carried  out  by  a  sub- 
routine other  than  the  one  in  which  the 
symbol  appears.   In  the  detailed  flow 
charts  the  information  needed  by  the 
routine  being  called  is  written  inside 
the  block. 
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(jTTy^ 


Decision  point  in 
to  be  followed  is 
or  NO. 


the  logic .   The  route 
indicated  by  the  YES 


A  logical  Jump  Is  to  be  made  to  some 
point  within  the  subroutine. 


Beginning  point  of  a  subsection  of  logic 
within  a  subroutine.   A  logical  Jump  will 
always  go  to  such  a  point.   The  numeral 
indicates  page  number. 


© 


Continuation  of  the  logical  flow  to  the 
next  page.   This  symbol  would  appear  at 
the  bottom  of  page  one  and  at  the  top  of 
page  two . 


3nd  of  event  subroutine, 
process  the  next  event. 


Program  is  to 


EXIT 


End  of  an  event  subroutine  when  for  any 
reason  the  normal  program  logic  is  to  be 
interrupted. 


RETURN 


Logical  end  of  any  subroutine  which  is 
not  an  event  routine.   Indicates  a  Jump 
back  to  the  section  of  the  program  that 
called  the  subroutine  into  action. 


STOP 


Instruction  to  stop  the  game.  Indicates 
that  all  the  data  has  been  processed  or 
a  mistake  was  made  in  the  input  informa- 
tion. 


5 

NAME 


( PIO . 2 ) 


Data  is  to  be  read  into  the  machine.   In 
the  detailed  flow  charts  the  top  numeral 
indicates  the  tape  unit  to  read  from  and 
the  characters  in  parentheses  indicate  the 
format  to  be  used.  NAME  is  the  variable 
name  to  be  read. 
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-  Write  output  on  tape  unit  six.   The  num- 
eral at  the  bottom  Is  a  code  Indicating 
the  format. 


2jTg       I  -   Store   New  Event.      Type   of  event  and 

event  time  will  appear  within  the 


block. 


Symbols  used  within  blocks  in  detailed  charts. 

RN  -  Random  Number 

♦  -  Multiplication 

/  -  Division 

**  -  Exponentiation 

A+B  — >  C       -  Add  A  to  B  and  store  the  answer  in  the 

variable  name  G. 

List  of  variable  names.  Following  is  a  list  of  the  variable 
names  with  their  definitions.   The  names  are  arranged  according 
to  the  subroutine  in  which  they  were  originally  defined.   This 
also  corresponds  to  the  order  of  the  COMMON  statements.   Input 
data  names  are  designated  by  an  asterisk.   Names  representing 
data  to  be  printed  are  identified  by  0  or  #  symbols.   Dimensioned 
variables  are  indicated  by  parentheses  following  the  nsune . 

(1)   Utility  Words 

I'^  Used  for  temporary  storage  or  con- 

Tl  venience  throughout  the  program. 

T2 


IDIP  Helicopter  dip  counters.   12  is  set 

12  equal  to  the  total  number  of  dips 

(number  of  helicopters  times  num- 
ber of  dips  per  helicopter)  at  the 
beginning  of  each  play.   IDIP  is 
incremented  each  time  a  helicopter 
dips.   If  a  detection  does  not  occur 
the  play  is  over  when  IDIP  =  12. 
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J2 


Counter  for  number  of  samples.  Used 
In  computing  average  number  of  de- 
tections per  sample. 


(2)  Main  Control  Program 


*NRSAM 


Number  of  samples  in  each  simula- 
tion run. 


♦NSSIZE  (I) 


Number  of  plays  constituting  the 

-.th     , 
I   saunple. 


$NHSDET  (I) 


Counter  for  keeping  track  of  the 
number  of  detections  made  during 
one  sample  by  the  ^th  helicopter. 


$TMIN 


Minimum  time  in  minutes  from  the  time 
the  helicopters  arrive  at  datum  until 
a  detection  occurred  for  current 
sample . 


$TMAX 


Maximum  time  required  for  a  detec- 
tion during  any  sample. 


8BART 


Average  time  elapsed  between  time 
helicopters  arrive  at  datum  and 
detection  when  detection  occurs. 


$VART 


Sample  variance  of  detection  time 
when  detection  occurs. 


#MIND 
^MAXD 
#BARD 


Minimum  number  of  plays  resulting  in 
detections  during  a  series  of  runs. 

Maximum  number  of  plays  resulting  in 
detections  during  a  series  of  runs. 

Average  number  of  plays  resulting 
in  detections  per  sample  during  a 
series  of  runs. 


$ 


Printed  after  each  sample  is  completed. 


# 


Printed  at  the  end  of  runs  designated  by  game  user. 
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♦  YARD 

$   PERO 

♦  DAT 

♦  XDAT 

♦  YDAT 
I  NDET 

NDTEMP 


PL(I) 


♦    IRNO 


Sample  variauace  of  number  of  detections  per  sample 

Percentage  of  plays  resulting  in  a  detection  dur- 
ing a  sample. 

Time  at  which  helicopters  arrive  at  datum  and 
begin  search. 


X  and  Y  coordinates  of  datum. 


Counter  which  is  set  to  zero  at  beginning  of  each 
sample  and  incremented  each  time  a  detection  occurs. 

Temporary  counter  which  is  set  equal  to  NDET  at  the 
beginning  of  each  play.   Comparison  of  NDET  and 
NDTEMP  signals  whether  the  play  ended  because  a 
detection  occurred  or  because  the  helicopters 
completed  all  dips. 


One  way  sonar  propagation  loss  of  the  I 
Computed  at  the  beginning  of  each  play. 


th 


helicopter. 


Number  of  pseudo-random  numbers  to  be  generated 

and  disposed  of  at  the  beginning  of  a  simulation  run. 


(3)   Subroutines  SNE  and  TNE 


The  following  three  words  are  used  throughout  the  program 
to  transfer  information  to  or  from  the  Event  Store  Table. 

TIMET    Time  at  which  an  event  is  to  be  stored  or  executed. 

NREVT    Number  of  the  event  to  be  stored  or  executed. 

NRUNT    Number  of  the  unit  affected  by  the  event. 

NTNE     Counter  which  keeps  track  of  the  location  of  the 
last  event  in  the  table. 

TIME(I)   These  three  arrays  comprise  the  event  store  table, 
NREV(l)   and  represent  time  of  execution,  number  of  the 
NRUN(I)   event  subroutine  to  be  called,  and  the  number  of 
the  unit  involved  in  the  event. 
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(4)   Subroutine  SME  (Submarine  Maneuver  i:vent) 

KT        Submarine  maneuver  counter. 
*  NRAl^SS 


«   MAN 


STIME 
SMNEXT 

DT 

VXS 

VYS 

XS 

YS 

TT 


Signals  vfhether  submarine  maneuvers  are  predeter- 
mined or  random. 

Signals  whether  the  first  course  is  randomly   se- 
lected when  submarine  maneuvers  are  predetermined. 

Time  the  submarine  last  maneuvered. 

Time  at  which  the  submarine  is  scheduled  to  make 
the  next  maneuver* 

Depth  at  STIME. 

V^  at  time  STIME. 

V  at  time  STIME. 

X  and  Y  coordinates  at  time  STIME. 

Temporary  storage  for  random  numbers  which  are  to 
be  added  to  submarine  courses. 


th 


Submarine  maneuver  table. 

*  SMTIME(I)   Time  to  execute  the  I ^"  submarine  maneuver. 

♦  DE?TH(I)    Depth  (feet)  at  time  SMTIME(I). 

Course  (degrees)  at  time  3MTIMS(I). 
Submarine  speed  (knots)  at  time  SMTIME(I). 


♦  scusd) 

*  SSPD(I) 

*  NR14AN 


Number  of  times  the  submarine  is  to  maneuver 
when  maneuvers  are  predetermined. 


Input  parameters  when  submarine  maneuvers  are  random. 


*  NHISPD 

»  LOSPD 

«  MAXSD 

*  MIN3D 


Upper  and  lower  limits  (knots)  of  submarine  speed. 
Upper  and  lower  limits  (feet) of  submarine  depth. 
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Cavitation  speed  versus  depth  table 


«  NCAVS(I) 

♦  NCAVD(I) 

«  NCAV 


Cavitation  speed  in  knots  for  the  I"^^  point 
on  the  curve. 

Minimum  depth  corresponding  to  NCAVS(I). 

Number  of  entries  in  the  cavitation  table. 


(5)   Subroutine  HSM  (Helicopter  Maneuver  Event). 


IDIPN(I) 

♦  NRDIPS 

♦  HTE(I,J) 

♦  HBRG(I,J) 

XH(I) 
YH(I) 

♦  NRHS 

♦  H3PD 

HSPDT 

♦  TDD 

ERB 

ERT 

XHT(I) 

YHT(I) 

♦  NTE 


Counter  which  keeps  track  of  the  number  of  dips 
made  by  the  I"^^  helicopter. 

Number  of  dips  to  be  made  by  each  helicopter. 

Jump  time  of  I"^^  helicopter  enroute  to  the  J^^ 
dip  station.   Input  as  yards  and  converted  to 
minutes  internally. 

Bearing  relative  to  the  preceding  leg  flown  by 
the  I^^  helicopter  in  transiting  to  the  J^^  dip 
station.   Input  as  degrees  and  converted  to 
radians  internally. 

X  and  Y  coordinates  of  the  I^h  helicopters  last 
dip  station. 

Number  of  helicopters. 

Speed  in  knots  at  which  helicopters  transit 
between  dip  stations. 

HSPD  converted  to  yards  per  minute. 

Time  (minutes)  required  to  establish  a  hover 
and  lower  the  transducer. 

Error  in  bearing  added  to  HBRG(I,J)  to  de- 
termine actual  flight  path  of  helicopters. 

Distance  error  added  to  HTE(I,J)a 

Table  for  temporary  storage  of  helicopter  X 
and  Y  coordinateso   Used  whenever  ERB  and  ERT 
are  not  computed. 

Signals  whether  ERB  and  ERT  are  to  be  computed. 
Value  is  one  if  errors  are  to  be  computed,  and 
zero  otherwise. 
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(6)   Subroutine  HDE  (Detection  Event) 


DRNG(I,J) 
»   NRSJ^ 

*  S'/fT(I) 

*  PD(I) 

*  TRD 

NOVER 

*  SHORTR 

NSWP(I) 

*  PRPO 

«    PGTMAX 

PCT 


Detection  range  of  the  I^^  helicopter  for  the 
Jth  sonar  sweep. 

Number  of  sonar  sweeps  to  be  made  by  each 
helicopter  at  each  dip  station. 

Time  (minutes)  required  to  complete  the  I 
sonar  sweep. 

Transducer  depth  for  the  I^   sonar  sweep. 

Time  (minutes)  required  to  retrieve  the  trans- 
ducer at  the  end  of  each  dip. 

Signals  that  a  play  has  been  completed.   Value 
is  one  if  play  is  over,  zero  otherwise. 

Range  (yards)  from  helicopter  within  which 
target  is  too  near  for  detection  to  occur. 

Signals  the  first  sweep  number  to  be  executed. 

Probability  of  detecting  a  non-submarine  contact. 

Maximum  time  for  evaluation  of  non-submarine 
contacts. 

Actual  time  to  evaluate  non-submarine  contact. 


(7)   Subroutine  COMPRG 


« 

NS 

•» 

T 

# 

P 

DP 

« 

DL 

# 

TS 

« 

POMMU 

« 

FOKVAR 

Sea  State,  (beaufort  scale). 

Temperature  in  the  layer,  (degrees  Fahrenheit). 

Acoustic  frequency  of  the  helicopter  sonar. 

Transducer  depth*  (feet). 

Layer  depth,  (feet). 

Target  strength,  (decibels). 

Average  figure  of  merit  of  the  helicopter  sonar, 
(decibels) . 

Variance  of  helicopter  sonar  figure  of  merit, 
(decibels) . 
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Main  Program*   The  main  AHS-1  program  performs  the  functions 
of  moderator  and  bookkeeper.   This  section  of  the  program  calls 
for  more  Input  data  when  It  Is  needed  and  records  the  results  of 
each  play.   The  main  program  also  sets  up  the  Initial  submarine 
and  helicopter  maneuver  events  so  that  each  replay  of  the  game 
will  begin  at  the  proper  point  in  time  and  space. 

The  following  terms  will  be  referred  to  throughout  the  de- 
scription of  the  program. 

Play.   A  play  of  the  game  begins  when  the  submarine  leaves  datum 
at  game  time  zero  and  ends  either  when  the  submarine  is  detected 
by  a  helicopter  or  when  all  helicopters  have  completed  their  last 
dip. 

Sample .   A  sample  consists  of  a  series  of  plays.   Sample  size  is 
determined  by  the  user  and  is  ordinarily  based  on  what  Is  consid- 
ered necessary  for  statistical  validity  of  the  results. 
Run.   A  run  consists  of  one  or  more  samples. 

At  the  completion  of  a  run,  any  number  of  input  parameters 
may  be  changed.   Inputs  do  not  vary  between  samples  of  the  same 
run,  additional  samples  representing  replications  of  the  same  ex- 
periment.  The  only  difference  between  two  samples  of  the  same 
run  is  the  series  of  random  numbers  used.   The  random  number  gen- 
erator is  reset  to  the  first  number  in  the  series  at  the  beginning 
of  each  run. 

The  various  functions  of  the  main  program  should  be  apparent 
from  the  accompanying  flow  diagram  once  the  reader  is  familiar 
with  each  individual  subroutine.   Only  those  parts  which  are  not 
immediately  clear  will  be  mentioned  here. 


74 


It  is  contemplated  that  the  usual  submarine  operating  mode 
will  be  a  completely  random  submarine,  option  (3).   The  datum 
coordinates  will  normally  be  the  origin  of  the  X,Y,  coordinate 
system.   For  this  reason  variable  names  associated  with  these 
Inputs,  referred  to  as  optional  variables  in  the  flow  chart, 
are  preset  prior  to  the  first  run.   This  obviates  the  necessity 
of  reading  these  inputs  in  except  when  other  than  this  standard 
operating  mode  is  desired. 

The  minimum  number  of  detections  ajid  maximum  time  to  detect 
are  initially  set  at  large  numbers  so  they  will  always  be  replaced 
by  the  actual  quantities  once  a  detection  occurs.   If  the  sub- 
marine is  not  detected  these  qauantities  will  be  printed  as  1000 
detections  and  600  minutes  respectively  and  should  be  disregarded. 
If  a  sample  size  of  1000  or  larger  is  used,  the  main  progrsun 
should  be  chsmged  to  reflect  this  by  setting  MIND  equal  to  a 
number  larger  thsm  the  sample  size. 

Bearing  and  course  inputs  to  the  program  are  measured  clock- 
wise, courses  being  relative  to  north, or  000.   For  computational 
convenience  the  zero  radial  is  congruent  with  the  X  axis  and  all 
bearings  are  measured  counter-clockwise.   For  purposes  of  the 
game  itself  this  presents  no  difficulty  since  only  relative  dis- 
tances and  bearings  affect  the  simulation  results.   With  no  ad- 
justment the  resulting  game  would  be  a  mirror  image  of  what  one 
would  plot  from  the  inputs,  and  relative  distances  would  be  pre- 
served.  In  anticipation  that  the  user  might  wish  to  have  unit 
positions  printed  out  at  some  time  during  the  play,  an  adjust- 
ment has  been  made  to  all  courses.   By  using  the  negative  of 
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courses  as  input  to  the  program  and  adding  ninety  degrees,  the 
coordinate  system  has  been  transformed  into  a  geographical  map. 
Consequently  submarine  and  helicopter  positions  will  agree  with 
the  usual  map  representation. 
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MAIN   PROGRAM 


(General) 


START 


SET   OPTIONAL   INPUT 
VARIABLES    TO    ZERO 


SET  CONSTANTS    FOR 
SONAR  RANGE 
COMPUTATIONS 


SET  FIRST  ENTRY 
IN  EVENT  STORE 
TABLE    TO    ZERO 


INITI-ILIZE   DETECTION 
STATISTICS 


O 


d) 
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© 


Main   Prog.      (G) 
2   of  6 


(is    A  NFrf   RUN   required)—^ 


-^ 


STOP 


YES 


READ   IN   DATA 
AND  PRINT 
NEW   RUN    INPUTS 


±. 


INITIALIZE  NEXT   SUBMARINE  MANEUVER 
TIME   SIGNAL    FOR   RANDOM   SELECTION 
OP   SUBMARINE-  COURSE    AND   SPEED 


SET   SUBMARINE  COURSE 
TO    ADD  90    DEGREES    TO 


ADJUSTMENT 
ALL  COURSES 


ARE   RANDOM    ERRORS 
TO   BE    ADDED   TO 
HELICOPTER  JUMP 
BEARINGS    AND  RANGES 


YES 


^0    TO" 
B3. 


NO 


COMPUTE 
DIP 


AND   TABULATE 
STATIONS 


© 


Main  Prog.    (G) 
3  of  6 


B3 


■^ 


J^ 


INCREMENT   SAMPLE  COUNTER 


:J^ 


RECORD   AVERAGE   FIGURE   OP 
MERIT    AND   POM   VARIANCE 


RECORD   TOTAL  NUMBER 

OP  DIPS    AND   SET    TOTAL 

DIP  COUNTER   TO    ZERO 


RESET    AND  CYCLE 

RANDOM  NUMBER 

GENERATOR 


BEGIN   NEW 
SAMPLE 


C3 


INITIALIZE 
HELO  DIP 


INDIVIDUAL 
COUNTERS 
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© 


Main  Prog.    (G) 
4  of  6 


INITIALIZE  MINIMUM 
MAXIMUM,    MEAN    AND 
VARIANCE   FOR 
DETECTION    TIME 


BEGIN 
NEW  PLAY 


<r 


€) 


^ 


INITIALIZE  SUBMARINE  DEPTH 


SET  UP  FIRST 

SUBMARINE  MANEUVER 

AND  STORE  IN  EVENT 

STORE  TABLE 


^ 


SET  UP  FIRST  DIP 
FOR  EACH  HELICOPTER 
^AND  COMPUTE  PROPAGATION 
LOSS  FOR  EACH  HELO 


J^ 


INITIALIZE   END  OF  GAME    FLAGS 


D 


© 
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Main  Prog.  (G) 
5  of  6 


TAKE   FIRST  EVENT. 
TNE  WILL  HAVE  CONTROL 

UNTIL  DETECTION  OR 
COMPLETION  OF  LAST  DIP 


NO 


^ 


-iL 


AS   SUBMARINE   DETECTED 


0 


YES 


INCREASE    INDIVIDUAL   HELO 
DETECTION  COUNTER 


COMPUTE    SQUARE   OP   TIME    TO 

DETECTION    AND  RECOMPUTE 

MINIMUM    AND  MAXIMUM 

TIME    TO    DETECTION 


IS    THIS    THE   LAST 
'LAY   OF   THIS   SAMPLE 


NO 


YES 


IJO    T(r 
.D4 
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Main  Prog.    (G) 
6  of  6 


COMPUTE  MEAN    AND  VARIANCE 
OP   TIME    TO   DETECTION 


COMPUTE   PERCENTAGE 

OP  PLAYS   RESULTING 

IN    SUBMARINE   DETECTION 


COMPUTE   SQUARE   OF  NUMBER 
OP   DETECTIONS    AND  RECOMPUTE 
MINIMUM    AND  MAXIMUl^   OP 
DETECTIONS    PER   SAMPLE 


IS    THIS    THE   LAST 
lAMPLE    POR    THIS   RUN 


NO 
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Subroutine  SNE.   Subroutine  SNE  (Store  New  Event)  enters 
future  events  in  chronological  order  In  the  event  store  table. 
The  event  store  table  Is  a  llstp  ordered  according  to  time,  of 
actions  which  are  to  take  place  in  the  course  of  a  play  of  the 
gajiie. 

VThen  an  event  is  to  be  stored,  SNE  is  entered  with  three 
Items  of  information;  the  time  at  which  the  event  Is  to  occur, 
the  type  of  event  to  be  executed,  and  the  number  of  the  unit  In- 
volved.  This  information  constitutes  one  event  word. 

Should  two  or  more  events  be  scheduled  to  occur  at  the  same 
time  there  is  no  designated  ordering  according  to  event  type. 
Two  or  more  such  events  will  take  place  in  the  order  in  which 
they  were  processed  by  SNE. 
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Subroutine  SNE 


(General) 


ENTER 


Ak 


COMPARE  EXECUTION  TIME 
OP  EVENT  TO  BE  STORED 
WITH  EVENT  TIMES  IN 
EVENT  STORE  TABLE 


nLI 


INSERT  EVENT  INFORMATION 
(TIME,  UNIT  NUMBER,  EVENT 

NUMBER)  IN  PROPER 
CHRONOLOGICAL  LOCATION 
IN  EVENT  STORE  TABLE 


JL 


INCREASE   EVENT 
COUNTER   TO   REFLECT 
ADDITION   OP  EVENT 
TO   EVENT    STORE   TABLE 


\/^ 


RETURN 
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Subroutine  3NE 


(Detailed) 


ENTER 


I  =  NTNE 
NTNE  =  NTNE+I 


NO 


friMET  <  TIME(I)V 


YES 


I  =■  I-l 


TIME(I+l) 
NRUN(I^-l) 
NREV(I+1) 


TIME(I) 
NRUN ( I ) 
NREV(I) 


TIMEd-t-l)  =  TIMET 
NRUN(i-»-1)  =  NRUNT 
NREV(I+1)  =  NREVT 


RETURN 
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Subroutine  TNE«   Subroutine  TNE  (Take  Next  Event)  examines 
the  event  type  of  the  earliest  entry  in  the  event  store  table 
and  calls  the  appropriate  event  subroutine  into  action.   Once  a 
play  has  commenced,  the  program  Is  controlled  by  TNE  until  the 
end  of  the  play.   The  end  of  a  play  is  signaled  by  the  detection 
event  subroutine.   VThen  this  signal  is  received,  TNE  trsinsfers 
control  back  to  the  main  program. 

When  an  event  is  executed  the  event  store  entry  referring  to 
It  is  removed  from  the  table.   An  alternative  method  would  be  to 
leave  the  entries  in  the  table  and  keep  track  of  the  earliest  un- 
executed event.   The  method  used  is  somewhat  more  time  consuming 
but  results  in  a  considerable  saving  of  storage  space. 
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Subroutine      THE 


(general) 


EUTER 


RETURN 
CONTROL 
TO  MAIN 
PROGRAM 


YES 


1<r 


HAS   A 
DETECTION 
OCCURRED 


NO 


STORE  EARLIEST 
EVENT  WORD  IN 
TEMPORARY  STORAGE 


V 


SHIFT  STORED  EVENT 
WORDS  DOWN  ONE  CELL 
IN  EVENT  STORE  TABLE 


iL 


CALL  NEXT  EVENT 
SUBROUTINE  CORRESPONDING 
TO  EVENT  WORD  IN 
TEMPORARY  STORAGE 
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Subroutine   TNS 


(Detailed) 


RETURN 


<r 


ENTER 


YES 


NOVER  ^  0 


NO 


TIME(2)  -^  TIMET 
NRTJN(2)  «^  NRUNT 
NREV(2)  -»  NREVT 


:iL 


3  ->    I 


TIME(I)  -^  TIME(I-l) 
NRUN(I)  -»  NRUN(I-I) 
NREV(I)    ->   NREV(I-l) 


^ 


NTE-1    -^   NTNS^^^^^    I   I  NTNE^ ^ >j  It-l 


-»    I 


(   NREVT  =   1  Jp:^^ 


GALL 


SME 


c 


NREVT  = 


STOP 


ny 


lES      ./    CALL 
HSM 


NREVT  =   3 
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Subroutine  SMS»   The  first  time  this  subroutine  (Sub- 
marine Maneuver  Event)  is  called  on  during  a  play,  the  sub- 
marine operating  mode  is  determined  from  the  information  in- 
put by  the  game  planner.   Once  this  determination  has  been 
made,  course,  speed  and  depth  are  determined  randomly  or  read 
directly  from  the  input  information o   Game  time  and  the  sub- 
marine position  are  noted  and  X  and  Y  components  of  velo- 
city are  computed.   From  these  four  pieces  of  information  the 
submarine's  coordinate  position  can  be  determined  at  any  later 
time. 

Additional  functions  performed  by  Subroutine  SME  are 
as  follows: 

(1)  On  the  first  submarine  maneuver  or  any  subsequent 
maneuver  involving  a  depth  change 9  helicopter  de- 
tection ranges  are  computed  by  calling  Subroutine 
COMPRG. 

(2)  If  another  submarine  maneuver  is  to  take  place,  the 
time  the  maneuver  is  to  occur  is  noted  and  subroutine 
SNE  is  instructed  to  place  a  submarine  maneuver  event 
In  the  event  store  table. 

(3)  The  time  of  the  next  submarine  maneuver  is  recorded 
for  use  by  the  detection  event  subroutine.   This  in- 
formation is  used  in  determining  whether  the  sub- 
marine maneuvered  during  any  helicopter  dip.   In  the 
event  of  a  submarine  maneuver  while  a  helicopter  is 
dipping,  submarine  track  and  depth  information  must 
be  updated  before  the  dip  is  completed. 
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Subroutine  3ME  -  Submarine  Maneuver  Event 


(General ) 


ENTER 


ML 


IS  THIS  THE 
FIRST  MANEUVER 


_^ 


YES 


^J^ 


Is   THIS    A  DEPTH 
CHANGE   ONLY 


NO 


ARE    SUBMARINE 
COURSES   RANDOM 


sTES 


GO    TO 
B2 


NO 


)iL 


-^TS    FIRST   COURSE' 
V  RANDOM 


YES 


±- 


^ 


ADD   RN   TO 
HLJj  COURSES 


^ 


COMPUTE   RN 
BETWEEN 
0-360 


\/_ 


COMPUTE   NEW   X, 
Y,V^,V^   and 

RECORD  TIME 


® 


^L 


-Nli. 


IS  THIS  A 
DEPTH  change; 


YES 


NO 


/  COMPUTE 
/HELICOPTER 
■^DETECTION 
\     RANGES 


^ 


->f     IS    THIS    THE 
LAST  MANEUVER 


NO 


SNE 
EVENT  1    AT 
vNEXT  MANEUVER; 
TIME 


YES 


90 


^ 


COMPUTE   ^N^ 
BETWEEN 
0-360 


Jl. 


3UBM/IRINE 
:  CURSE   =    RN1 


J^ 


IS    SUBM/iRINE 
NON-OAVITATING 


NO 


JiL 


COMPUTE  RN 
BETWEEN 
SPEED  LIMITS 


V- 


SPEED  =  MINDIUl'I 
SPEED   -f   RN 


J^ 


COMPUTE   RN 
BErv^EEN 
^DEPTH   LIMITS/ 


-V- 


DEPTH  =  MINIMUM 
DEPTH   +   RN 


C2 


€> 


SME      (G) 
2  of  2 


ICMPUxE   RN 
BETWEEN 
vSPEED  LIMITS/ 


-\k_ 


SPEED  =  MINIMUM 
SPEED  +  RN 


^ 


DETERMINE 

MINIMUM 

DEPTH  FOR 

NO  CAVITATION 


iOMPUTE  RN 
BETWEEN 
VDEPTH  LIMITi 


_y: 


DEPTH   =  MINIMU>: 
DEPTH   +   RN 


^ 


:^ 


COMPUTE   VELOCITY 
COMPONENTS 
V     and   V 


Jl 


COMPUTE 
HELICOPTER 
DETECTION 
RANGES 
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Subroutine  3ME- Submarine  Maneuver  Ivent  (Detailed) 


<  ^^•^rNRAN3S=0^ 


160   -»     RNC 
GALL   RAND 


RN 


TT 


<^ 


360  -»     RNO 
CALL      RA^TD 


i(L 


RN/57. 29573  -»  Tl 


:J^ 


NOAV  =  0 


YES 


^ 


'  NHISPD-L03PD 
+1  ->     RWO 

,     CALL   RAND 


^ 


(L0SPD+RN)*33.766 
->  T2 


±. 


MAXSD-MIN3D 
+1  ->•     RNC 
CALL   RAND 
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NO 


SME      (D) 
2   of   3 


-(^33US(KT)    =    30U3(KT«1)^ 


ye5 


S3PD(KT)    =    3SPD(KT-1) 


YES 


S3PD(KT)«33. 7666666  ->3PD 
(-3:US(KT)-^TT)/57. 29573   -»  0U3 
XS+(SMTIME(KT)-3TIME)«VX3   ->   XS 
YS+(SMTIME(KT)-3TIME)4*VTS   ->   YS 
SPD*  cos(CUS)-»  VXS 
SPD+f   sln(CU3)-»  VT3 
SMTIME(KT)->  3TIME 


YES 


DEPTH (KT)    =   DT 


NO 


Ad. 


DEPTH(KT)-^   DT<: 


^ 


CALL 
COMPRG 


->  KT  ^   NRM^N) 


YES 


NO 


>IL 


KT+1    -»  KT 
SMTIME(KT)->   TIMET 
l->  NREVT 
CALL   SNE 
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\L 


NC  AVS  ( NC  AV )  -NO  AVS'(  1 ) 
•►1    -^      RNC 
GALL  RAND 


N3AVS(1)+RN  -^  NSS 
NSS«33.7666  ->  T2 


^ 


(jISS  ^  NCAVS(li)-^I-»«l 


YES 


N3  AVL  ( NO  AV )  -  NO  AVD  ( I ) 
+1     -^      RNO 
CALL   RAND 


NOAVD(I)+RN   ->   DT 


Sl^IE      (D) 
3   of   3 


ID3 


-M- 


T2*cos(Tl)->  VXS 
T2«sln(Tl)-»  VYS 


/CALL 
^\^OMPRG 


Xgo  ton. 
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Subroutine  HSM.   Subroutine  HSM  (Helicopter  Maneuver  Event) 
updates  each  helicopter's  dip  station  coordinates  and  records 
this  information  for  use  by  the  detection  event  subroutine.   As 
the  program  is  presently  written,  dip  coordinates  are  computed 
by  the  main  program  and  stored  in  a  table.   This  table  is  avail- 
able to  Subroutine  HSM  whenever  a  helicopter  is  to  be  maneuvered. 
This  requires  that  the  computations  be  made  only  once  for  each 
game  run.   Should  the  necessary  code  for  computing  distance  and 
bearing  errors  be  >rritten,  new  dip  coordinates  will  be  computed 
by  Subroutine  HSM  every  play. 

Additional  functions  of  this  subroutine  are  as  follows: 

(1)  Records  the  number  of  the  last  dip  station  occupied  by 
each  helicopter.   After  the  first  dip,  the  order  of 
maneuvering  is  determined  by  dip  cycle  time  rather 
than  helicopter  unit  number.   Therefore  ar  individual 
dip  counter  must  be  maintained  for  each  helicopter. 

(2)  Stores  the  next  detection  event  for  the  helicopter 
which  has  maneuvered.   The  succeeding  detection  event 
will  occur  at  the  time  the  maneuver  is  completed  plus 
the  time  which  elapses  while  the  transducer  is  being 
lowered  into  the  water.   The  latter  time  Includes  normal 
delay  time  required  to  transition  from  forward  to  hover- 
ing flight  and  is  an  input  to  the  simulation. 
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Subroutine  H3M  -  Helicopter  Maneuver  E/ent 


(General) 


ENTER 


DETERMINE  NUMBER 
OP  HELICOPTER 
TO  BE  MOVED 


INCREMENT 
DIP  COUNTER 


ARE  ERRORS  IN 
:OURSE  AND  DISTANCE 
TO  BE  COMPUTED 


NO 


READ  NEW   X    AND 
Y   COORDINATES 
FROM    TABLE 


YES 


±- 


COMPUTE  COURSE 
AND  DISTANC 
ERRORS 


ICE 


GO    TO 
.A2. 


COMPUTE  NEVf   X    AND 
Y  COORDINATES 
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HSK  (G) 
2  of  2 


COMPUTE 
TIME  ENROUTE 


ADD  TIME  ENROUTE 
AND  TIME  TO 

LOWER  TRANSDUCER 
TO  GAME  TIME 


^ 


SNE 

DETECTION   EVENT 
FOR    THIS   HELO  TO 
0:CUR    AT    TIME 
COMPUTED   .IBOVE 


f    TKE     j 
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SubroutineHSI'.- Helicopter  Maneuver  ^vent     (Detailed) 


SN 


ER 


NRUNT   -♦   I 
IDlPN(I)fl    ->  IDIPN(I) 
IDTPNd)    ->  J 


\L 


(  NTE  ^0        J- 


YES 


NO 


COMPUTE 
ERB  and  ERT 


HBRG(I,J)fERB   ->   Tl 
HTE(I,J)+ERT   -►    T2 


XH(I)-»-T2#cos(Tl)    ->  XH(I) 
YH(I)-»-T2*sln(Tl)    -»  YH(I) 


XHT(I,J)    -»  XH(I) 
YHT(I,J)    -^  YH(I) 


TIMETfHTE(I,J) 
t-TDD   -»  TIMET 


GO^TO" 


'fc 


d) 
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T2/HSPDT   ->    Tl 


TIMET+TDDfTl 


TIMET 


A2 


3    -►  NREVT' 
1    -^  NRUNT 

CALL    SNE 


HSM    (D) 
2   of   2 
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Subroutine  HDS.   Subroutine  KDB   (Detection  Event)  Is  called 
Into  action  whenever  a  helicopter  arrives  at  a  new  dip  station. 
Using  the  information  previously  recorded  by  the  submarine  and 
helicopter  maneuver  events,  the  relative  positions  of  helicopter 
and  submarine  are  determined  at  the  beginning  of  each  sonar  sweep. 
On  the  basis  of  this  target  range  and  helicopter  detection  range, 
a  determination  is  made  as  to  whether  the  submarine  has  been  de- 
tected. 

Additional  functions  performed  by  Subroutine  HDE  are  as 
follows: 

(1)  Each  time  a  helicopter  dips  the  dip  is  recorded.   VThen 
all  helicopters  have  completed  their  assigned  number  of 
dips,  or  the  submarine  is  detected,  a  flag  is  set  which 
signals  the  end  of  the  play. 

(2)  If  the  play  terminates  with  a  submarine  detection,  the 
elapsed  time  from  when  the  helicopters  arrived  at  datum 
until  detection  time,  and  the  number  of  the  detecting 
helicopter  is  recorded.   This  Information  is  used  by 
the  main  progrsun  in  computing  detection  statistics. 

(3)  A  record  of  the  time  of  the  next  intended  submarine 
movement  is  maintained.   If  a  submarine  maneuver  is 

to  take  place  during  any  helicopter  dip,  the  detection 
event  is  discontinued  after  the  sonar  sweep  during 
which  this  maneuver  is  to  occur.   After  the  submarine 
maneuver  has  been  executed  the  helicopter  dip  Is  re- 
sumed at  the  same  game  time  at  which  the  interruption 
occurred. 

(^')   If  the  probability  of  contacting  a  non-submarine  target 
is  greater  than  zero  a  random  number  is  generated  and 
compared  with  this  probability.   On  the  basis  of  this 
comparison  and  the  maximum  time  to  pursue  non- sub- 
marine contacts  input  to  the  program,  appropriate  de- 
lays are  computed  and  added  to  dip  time. 

(5)   On  all  but  the  last  dip  for  each  helicopter,  a  maneuver 
event  is  stored  for  the  helicopter  currently  dipping. 
The  time  the  maneuver  is  to  take  place  is  determined 
by  adding  dip  time  to  the  game  time  at  which  the  dip 
commenced.   Dip  time  is  the  sum  of  all  sweep  times, 
non- submarine  contact  delay  time  and  the  time  required 
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for  the  helicopter  to  raise  the  transducer  and  make 
the  transition  from  hovering  to  fonrard  flight. 
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Subroutine   HDE  -   Detection  ^vent 


(General) 


Al 


-^ 


ENTER 

^f 

DETERMINE  WHICH 
HELO  INVOLVED 

_5C         ^ 

IS  THIS  HELO 
COMPLETING  A 
PREVIOUSLY 
INTERRUPTED  DIP 


YES 


DO   ONCE    FOR  EACH 
SONAR    SWEEP 


NO 


-> 


INCREMENT 
DIP  COUNTER] 


COMPUTE    TARGET 
POSITION   AND 
RANGE    TO    TARGET 


il- 


ls TARGET  RANGE 
^  SONAR  RANGE 


YES 


NO 


GO  t6^ 

.B2 
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TAJ^GET   RANGE 
>    MINIMUM 
DETECTION   R.INGE 


-                 ^ 

YES 

RS^^ORD   DETECTION 
AITD   TIME    THE 

DETECTION  oc::unRED 

\ 

f 

EXIT 

HDZ      (G) 
2   of   3 


ilDD   o'./EEP    TIME 
TO   GAME    TI>IE 


DID   SUBMARINE 
MANEUVER   DURING 
THIS    S^.VEEP 


NO 


YES 


IS    THI:) 

THE  LAST 
S'./EEP 


NO 


JiL. 


i.'j  i 


'DETECTION  EVSN1\ 
FOR   THIS   flELO 
AT   CURRENT 
GAME   TIME 


-^ 


SET  THIS  HELO'S 
FIRST  DIP  COUNTER 

TO  NUMBER  OF 
THIS  DIP  PLUS  ONE 


{    TNE  j 
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KDE      ( Ci ) 
3    of   3 


V 

IS   THIS    THE 

LAST    S'.fEEP 


YES 


IS    THIS    THE   LAST 
)I?   ?0R   THIS   HEL 


YES 


iL 


[AVE    .\LL   HELOS^ 
COMPLETED 
LAST    DIP 


NO 


■^   TNE     j 


YES 


RECORD  END 
OP   THIS    PLAY 


EXIT 


IS  PROBABILITY 

OF  NON- SUBMARINE 

CONTACT  >  ZERO 


YES 


T)ID  HELO   GAIN 
NON-SUBMARINI 

conta::t 


YES 


^ 


im. 


Ml 


COMPUTE    .ADDITIONAL 
DIP    TIME   REQUIRED 
FOR   PROSECUTION   OF 
NON- SUBMARINE  CONTACT 


3NE 

MANEUVER  EVENT  FOR 
'HIS  HELO  TO  OCCUR 
AT  GAME  TIME  PLUS 
TIME   TO    RAISE   DOME, 
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Subroutine      ED'2  -   Detection  Event 


(Detailed) 


ENTER 


V^ , 

:nur:r  ~»  i| 


V. 


NO 

-(n3WP(I) 

^0 

\ 

YES 

IDIFfl    -^  IDIP 

> 

^ 

s 

TIMET   ->    Tl 
T1-3TIME   ->   T2 

^ 

\ 

/" 

N3;7P(I)    ->   J 

XS+VXS<frT2    -^   XT 

Y3■^VYS^^T2    ->  YT 


© 


(XH(I)-XT)^+(YH(I)-YT)^ 
-^  RTT 


NO 


.GO    TO 


® 


RTT   r^  DRN( 


G(I.Jj) 


YES 


■^^^^-{riTT   >  '^HORTR    3 

YES 

> 

' 

1    -^  NOVER 
NDET-fl    -^  NDET 
0   -T>    IDIP 
T1^3'rfT(J)-D4T   ->  Tl 

> 

/ 

RETURN 
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Tl  +  3./T(J)-»  11 
T?  +  3:;t(J)->    T2 


^MN'XT  ^    Tl\ 


:ic 


YES 


/ ^^^^ N 

(j  -^  NR3:r     Y 


ye: 
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j+i  ->  iioM'Pd) 


3   ->  NREVT 
OilLL    3NE 


HD^.      (D) 
2   of   3 


^  J   ^  NR^W   ^ ^ 


J+1    ->   J 


YES 


-^PIPN(I)    ^   NRDIPs) 


YE.  3 


V- 


(iDIP  ^    12     ^ 


NO 


jiL 


1    ->  NOVEH 

0    ->»   IDIP 


JiL 


RETURN 


106 


hd:e:     (d) 
3  of  3 


±. 


FJPR    >  0 


NO 


YES 


CALL   RAN] 


rpCT  ^    POPR^ 


NO 


Tl+TRD  ->  TIMET 


YES 


(PGT*PGTMAJC)/P::!PR   -»►    PCT 
Tl  +  TRD+POT   ->   TIMET 


'2   -»   NREV: 
C:\LL    SNE 


1    ->  NS'.fPd) 
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Subroutine  PRINT,  This  subroutine  comprises  the  Input - 
Output  section  of  the  program.   It  is  called  at  the  beginning 
of  each  run  and  at  the  completion  of  the  last  run.   Input  data 
cards  are  arranged  in  groups  headed  by  a  card  having  a  Hollereth 
field  name  punched  in  the  first  six  columns.   This  name  signals 
the  information  to  be  read  from  the  succeeding  cards.   This  sys- 
tem allows  the  user  to  change  any  number  of  input  quantities  be- 
tween runs  without  preparing  data  cards  for  the  complete  set  of 
inputs.   After  the  inputs  for  each  run  have  been  read  into  the 
machine,  the  input  data  are  printed  out  and  output  data  headings 
are  printed. 

Input  data  and  rules  for  data  card  preparation.  Three  types 
of  fields  are  used;   External  Fixed  Point  (F),  Integer  (I),  and 
(A).   All  (F)  fields  are  F10.5  and  when  this  field  is  specified 
data  may  be  punched  using  either  one  of  two  methods;  (1)  with- 
out the  decimal  point,  requiring  that  the  last  digit  prior  to 
the  decimal  point  be  punched  in  the  fifth  column  of  the  field,  or 
(2)  with  the  decimal  point,  in  which  case  the  number  may  fall  any 
place  within  the  field. 

EXAMPLE:   To  read  in  the  numbers  1023.65  and  20.658  using 
a  2F10.5  field: 

(1)   Without  the  decimal  points 
Column 


1              5  6         10  11            15  16            20 

1|0|2   3|6|5   II     1      1    1    !2|0|6|5ia|   1    1 

(2)   With  the  decimal  point 
Column  1       5 


11012131.16151     I 


10   11  15 20 

I      I    I    I    \2\6rT^5m 


When  an  (I)  field  is  specified  the  last  digit  of  the  number  must 
fall  in  the  last  column  of  the  field.   The  (A)  field  is  used  for 
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reading  in  literal  characters  and  the  letters  or  numbers  to  be 
read  may  fall  any  place  within  the  six  columns  of  the  field. 

All  arrays  within  a  group  are  read  in  an  alternating  se- 
quence,  i.e.,   If  S'.fT  (I)  and  PD(I)  are  to  be  read  from  a 
6FIO.5  field,  I  =  1,4  the  following  order  is  used. 

Card  Columns  Name 

1  1-10  SWT(1J 

11-20  PD(1) 

21-30  SWT(2) 

31-40  PD(2) 

41-50  SiiT(3) 

51-60  PD(3) 

2  1-10  SVfT(4) 

11-20  PD(4) 

Two  dimensional  arrays  (Array  (I,J))  are  listed  on  the  card  in 
order  of  increasing  J  and  a  new  card  is  started  for  each  Increase 
in  I.   If  more  than  one  array  is  read  in  a  group p  the  names  are 
alternated  as  for  one  dimensional  tables.   In  the  following  in- 
structions only  the  first  elements  are  shownj  succeeding  elements 
will  be  in  the  order  shown. 
Data  Card  Format 

GROUP  I 

Card       Columns        Field        Name  Remarks 

1      1-6       A6    *DATE 

15-16     110     NDAY      Current  date. 
25-26  MONTH      Number  of  current  month. 

35-36  NYEAR      Last  two  digits  of  current 

year. 

This  Group  must  be  the  first  card  and  must  appear  only  once  In  a 
series  of  runs. 


Literal  characters. 
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GROUP  II 


Qard 
1 


Columns   Field 


1-5 
7-16 
26 


27-36 
2      1-6 
and  addl-7-12 
tlonal    . 
cards  if 
required 


A6 
110 


1216 


Name 

♦NRSAM 
NRSAM 
KD£T 


IRNC 
NSSIZE(l) 
NSSIZE(2) 


Remarks 


Number  of  samples. 
Always  equal  to  one  when 
included.   Signals  that 
detection  statistics  are 
to  be  printed  at  end  of 
run.   (S<=io  output  data) 


l*umber  ux  entries 
must  agree  with  NRSAM. 


GROUP  III 


Card       Columns        Field       Name 


Remarks 


1-6 

7-16 
17-26 
27-36 
37-46 


A6    *NHISPD     Submarine  speed  and  depth 

limits. 
no     NHISPD     Knots 
LOSPD      Knots 
MAXSD      Feet 
MINSD      Feet 


This  group  is  normally  not  Included  if  submarine  maneuvers  are 
predetermined.   However,  it  is  legal  to  include  both  groups  III 
and  V  in  the  same  set  of  data.   The  data  to  be  used  for  any  sub- 
sequent run  is  then  determined  by  the  value  of  NRANSS. 


Card 

Columns 

Field 

A6 
16 

GROUP  IV 
Name 

1 

1-6 
16 
22 

♦NRANSS 
NRANSS 
NRANC 

Remarks 


Equals  one  when  included. 
Equals  zero  if  random  number 
is  to  be  added  to  each  sub- 
marine course,  one  otherwise. 

This  group  need  not  be  included  unless  submarine  maneuvers  are 
predetermined.   The  one  exception  is  if  both  Groups  III  and  V 
have  been  read  on  the  same  run.   In  this  case  NRANSS  and  NRANC 
must  be  read  in  each  time  the  submarine  mode  of  operation  is  to 
be  changed.   A  zero  in  column  16  will  result  in  random  sub- 
marine maneuvers. 
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GROUP  V 


Card 


iolumris 

1-5 
7-16 


3 
and  as 
required 


11-20 

21-30 

31-40 

1-10 


Card 


Columns 

1-4 
7-16 

1-10 
11-20 
21-30 
31-40 
41-50 
51-60 


Field 
a6 

110 


Name 
*NRMAN 
NRMAN 


1-10  4F10.5        SMTIME(l) 


Field 

A6 
110 

6F10.5 


DEPTHd/ 

scusd; 

S3PD(1) 
3MTIME(2) 


GROUP  VI 
Name 

*NRS"»f 
NRS'.^ 

S'.fT(i) 
PD(1J 

s;;t(2) 

PD(2) 

3i/T(3) 

PD(3) 

GROUP   VII 


Remarks 


Number  of  submarine  maneu- 
vers. 

Minutes  (Time  the  maneuver 
is  to  take  place)  SMTIME(l) 
must  always  be  0.0. 
Depth,  (feet) 
Course,  (degrees) 
Speed,  (knots) 


Remarks 


Number  of  sonar  sweeps 

Minutes  to  complete  sweep. 
Transducer  depth,  (feet) ■ 


Card   Columns 


1-3 


Field   Name 


a6 


*TDD 


Remarks 


Card 


1-10 
11-20 
21-30 
31-40 


2F10.5 


Columns        Field 


1-4 
7-16 
17-26 

27-36 


A6 

110 


TDD 
TRD 
FCPR 
FCTMAX 


GROUP  VIII 

Name 

*NRHS 
NRHS 
NRDIPS 

HSPD 


transducero 
transducero 


Minutes  to  lower 
Minutes  to  raise 
0  ^  FCPR  r^  1.0 
Maximum  time  to  classify 
non- submarine  targets. 
(Minutes) . 


Remarks 


Number  of  helicopters. 

Number  of  dips  for  each 

helicopter. 

Helicopter  speed  (knots). 


Ill 


2      T-IO 

and  addi- 
tional cards 
as  needed. 

11-20 
21-30 

31-40 

41-50 

• 
• 

6P10.5 

HBRG(1,1) 

HTE(1,1) 
HBRG(1,2) 

HTE(1,2) 
HBRG(1,3) 

• 
• 

Bearing  of  first  hell- 
copter's  first  dip  sta- 
tion. (Degrees  relative 
to  000). 

Jump  distance  (Yards) 
Degrees  relative  to 
HBRG(1,1) 

Degrees  relative  to 
HBRG(1,2) 

New 
Card 

1-10 

11-20 

• 
• 
• 

6p10,5 

• 

H?RG(2,1) 
HTE(2,1) 

• 
• 
• 

GROUP  IX 

Card 

Columns 

Field 

Name 

Remarks 

1 

1-4 
7-16 
17-26 

A6 
110 

♦XDAT 
XDAT 
YDAT 

X  and  Y  coordinates 
of  datum. 

This  group  may  be  omitted  in  which  case  the  datum  coordinates 
will  be  (0,0). 


GROUP  X 

Card 

Columns 

Field 

Name 

1 

1-5 
7-16 

A6 

no 

*SONAR 

NS 

2 

1-10 

6F10.5 

FOMMU 

11-20 
21-30 
31-40 

FOMVAR 

DL 

T 

41-50 
51-60 

P 
TS 

GROUP  XI 

Card 

Columns 

Field 
A6 

Name 

1 

1-3 

♦DAT 

2 

1-10 
11-16 

F10.5 
A6 

DAT 
NAME 

Remarks 


Sea  State 

Average  sonar  figure  of 
merit.  (Decibels) 
Variance  of  FOM.  (db) 
Layer  depth.  (Feet) 
Temperature  in  the 
layer.  (^Fahrenheit) 
Sonar  frequency,  (kc.) 
Target  strength,  (db.) 


Remarks 


Time  late. (minutes) 
Alpha  —numeric  characters 
(Name  of  search  plan) 
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This  group  must  appear  in  every  run  and  must  not  be  placed 
ahead  of  Groups  I  or  II. 


Card 
1 


Columns 
1-2 


GROUP  XII 


Field 
A6 


Name 
♦GO 


Signals  program  to  print 
out  Inputs  and  commence 
new  run.  Must  appear  at 
the  end  of  each  set  of 
data. 


GROUP  XIII 

Card   Columns   Field   Name        Remarks 

1      1-5       A6    *FINIS      Signals  program  to  compute 

detection  statistics  at  the 
end  of  the  last  run  and 
stop  the  machine.   Need 
not  appear  if  detection 
statistics  are  not  desiredo 

Output  data.   The  output  statistics  reflect  the  writer's 
experience  in  ASW  Helicopter  squadrons.   No  assurance  is  given 
that  they  are  suited  to  any  other  potential  users  requirements. 
It  is  relatively  easy,  howeverp  to  write  the  coding  necessary 
to  compute  and  print  any  desired  output. 

The  following  output  data  is  printed  at  the  completion  of 
every  sample: 

(1)  Sample  size. 

(2)  Nuirber  of  detections.   This  is  the  primary  measure 
of  effectiveness  associated  with  this  simulation. 

(3)  Percentage  of  plays  resulting  in  detections.   This 
is  included  principally  as  an  aid  in  comparing 
samples  of  different  size. 


Percentage  = 


Number  of  detections 
Sample  size 


100 


(6.5) 


(4)   Average  time  to  detect  when  a  detection  occurs. 

The  time  to  complete  those  plays  which  do  not  end 
with  a  submarine  detection  are  not  included  in  this 
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^ 


average.   Time  to  detect  is  measured  from  the  time 
the  helicopters  arrive  at  datum.   A  problem  en- 
countered In  fleet  AS'.i  operations  is  when  to  aband- 
on a  search  once  begun.   Average  time  to  detect 
should  give  the  user  some  indication  of  the  time 
at  which  further  search  with  a  particular  search 
plan  is  apt  to  be  unproductive. 

(5)   Sample  variamce  of  time  to  detect. 

S^  =  (l/(n-l))  L^^-^    (X^  -  X)2  {6.6) 


where:   n  =  number  of  detections 

X,  =  time  to  detect  for  the  i   detection. 

X  ■  average  time  to  detect. 

2 
S  =  sample  variance  of  time  to  detect. 

(6)  Minimum  time  from  datum  time  to  detection 

(7)  Maximum  time  from  datum  time  to  detection.   Neither 
(6)  nor  (7)  include  plays  which  do  not  end  in  a 
submarine  detection. 

(3)  Number  of  detections  by  each  helicopter.  In  general 
considerable  discretion  must  be  used  in  drawirg  in- 
ferences from  individual  unit  performance.  Ireated 
with  due  caution,  however,  some  meaningful  informa- 
tion may  be  gained  from  this  knowledge.  This  part- 
icularly applies  to  asymmetrical  or  random  station 
search  plans. 

Printing  of  the  above  information  is  predetermined  by 
the  computer  coding  sind  the  user  has  no  control  over  whether 
the  data  is  printed  or  not.   The  succeeding  statistics  are 
computed  and  written  out  after  the  last  run  and  at  the  com- 
pletion of  any  run  determined  by  the  user. 

(1)   Average  number  of  detections  per  sample.   For  example 
If  this  is  printed  for  the  first  time  at  the  end  of 
run  two,  and  runs  one  and  two  each  consist  of  two 
samples. 

.4 


'^1=1  ^i 


D  =  (i)Z^,_T  D.  (6.7) 
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where:!)  =  average  number  of  detections  per  sample. 

D.  =  the  number  of  detections  for  the  i  sample. 

(2)  Sample  variance  of  number  of  detections  per  sample. 
Referring  to  the  above  examples 

^^ 
S^  =  (l/3)Z-i^i   (^i  "  ^^^  (^-^^ 

where:  S-.  =  sample  variance  and  D*  and  D  are  as  defined  above. 

(3)  Maximum  number  of  detections  taken  over  all  samples. 

(4)  Minimum  number  of  detections  taken  over  all  samples. 
The  above  quantities  have  meaning  only  if  all  samples 

are  of  the  same  size.   It  is  therefore 5  anticipated  that  the 
user  would  want  to  compute  and  print  this  data  whenever  sample 
size  is  changed.  Note  that  when  the  data  is  called  for  on  any 
run,  the  average  is  taken  over  all  runs  since  the  last  time 
these  quantities  were  printed. 
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Subroutine   PRINT 


General 


ENTER 


A 


NO 


IL 


COMPUTE 

DETECTION 
STATISTICS 


DETECTION 
STATISTICS 

_5^ 


IS    THIS   FIRST   RUN 


J^ 


<ARE   DETECTION^ 
STATISTICS 
TO   BE  COMPUTER 


YES 


W- 


NO 


LAST   RUN  COMPLETE 


NO  /  READ  NUMBER 


Ak 


YES 


OF  SAMPLES  & 
SAMPLE  SIZES 


STOP 


M 


V 


M. 


READ  NUMBER   OP 
SUBMARINE  MANEUVERS, 
TIME,    DEPTH,    COURSE    AND 
SPEED  OF  EACH  MANUVER 


SUBMARINE  COURSE 
SPEED  AND  DEPTH 
CHOSEN   RANDOMLY 


YES 


±. 


READ  HIGH  AND 
LOW  LIMITS  OF 
SUBMARINE  SPEED 
AND  DEPTH 


_.\l/. 


-^    READ  RANDOM 
SUBMARINE    AND 
RANDOM  COURSE 
INDICATORS 


-XT^    T0\ 
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PRINT   (G)   2  of  3 


0 


READ  NUMBER  OP  SONAR 
SWEEPS,  DEPTH  AND  TIME 
TO  COMPLETE  EA3H  SWEEP 


READ  TIMES  TO 
ESTABLISH  HOVER  AND 
RETRIEVE  TRANSDUCER 


READ  NUMBER  OP 
HELOS,  NUMBER  OP  DIPS, 
SPEED,  JUMP  BEARINGS 
AND  JUMP  DISTANCES 


^ 


CONVERT   HELICOPTER 
JUMP   DISTANCES   TO 
TIME   ENROUTE    AND 
BEARINGS   TO   RADIANS 


/■ 


READ   DATUM 
COORDINATES 


_^ 


READ   SONAR   INFO?MATION 
SEA   STATE,    AVERAGE   FOMp 
POM   VARIANCE,    LAYER   DEPTH « 
TEMPERATURE   IN   LAYER,    SONAR 
FREQUENCY    AND   TARGET   STRENGTH 


±. 


READ   DATUM 
TIME    AND   SEARCH 
PLAN  NAME 
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PRINT      (G)    3   of  3 


NEXT  OAJID 
FINIS 


NO 


NEXT  CARD 
LABELED  GO 


YES 


iL 


INCREASE  RUN 
NU14BER  COUNTER 


RETURN 
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Subroutine  PRINT 


(Detailed) 


ENTER 


Jl 


NAKE,N1,N?,N3,NA,N5   ^ 


(A6,5I10) 


Nl  ->  NDA.Y 
N2  ->  MONTH  1^ 
N3  ->  NrEAR 


y£5 


1  ->J1 

0  r^J2 

1  ->NRU 


6 
MONTH 
ND.VY 
NYEA.R 
TUifL 


Input   Data 
Title  '^ 


V 


0   ->  KDET 


<$^ 


jL 


NAME  =    DATE 


NO 


NO 


€> 


Al 


^AME  =    NRSAk)-^^^ ^KDET 


IZ3 


YES 


Nl   -»  NRSAM 
->N2   -»    KDET 


NSSIZE(I), 
I   =  1, NRSAM 


(1216) 


GO    TO 

ja. 
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PRINT   (D) 
2  of  9 


Jl  =  J2 


YES 


-^  Detectluu  Data 
Title 


NO 


6 
J1,J2 

Detection   Data 
Title 


JT   -►    T2 
B.IRD/T2    -*   BARD 
((VARD-T2«BilRD<t*2) 
/(T2-1. ))    -»  V4RD 


BARD,VARD, 
MIND,MAXD 


0.0 
0.0 
4321 
0 
0 


-> 
«► 

-» 
-» 
-> 


BARD 
YARD 
MIND 
MAXD 
JT 
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® 


PHINT      (D) 
3   of  9 


^ 


.NAME  =   NHI3PD 


yes 


Nl    -»   NR.1NSS 
N2   -»    NRA^'C 


NO 


Nl  ->NHISPD 
N2  -^LOSPD 
N3  ->MAXSD 
N4  -^'IIIITSD 


±. 


YES 
<r-^^ — (NAME  =  NRANSS 


NO 


<^0    T0> 


GO    TO 


NAME  =   NRM 


Nl    -^  NRSW  < 


YES 


^ 


YES 


Nl    -^  NRM  AN 


NO 


SMTIME(I),DEPTH(I), 
CCUS(I),3SPD(I), 

I   =   1, NRM  AN 


(4P10.5) 


ik- 


■(NAME  =   NRSV 


s-,;t(i),pd(i). 

I   =    l,NR3Vf 


(6F10.5) 


GO    TO 
Al 


NO 


k 


V. 


NAME  =    TDD 


yoo  Tb> 


y 


YES 


SRFO^PGTMAX 
Ti)D,TRD 


(4P10.5) 


NO 


® 
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PRINT      (D) 
4  of  9 


(^N^IE=NRHS    J^^^^ 


NO 


NAME=XDAT 


NO 


GO    TO 
Al 


NAME= SONAR 


YES 


Nl 


NCAV   ^ 


I£3. 


NO 


^OMMU,FOMVAR, 
DL,T,P,TS 


(6f10,5) 


NAME=NC  AV 


NCAV3(I),NCAVD(I) 
I    =   1,N3AV 


(1216) 


NO 


GO    TO^ 
.Al 


Jl 


ncavs{n:av) 

NCAVS(l)    - 

ncavd(n:av) 

NCAVD(l)    - 


->  NHISPD 

L03PD 
■^  MAXSD 

MINSD 


a" 
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0 


PRINT   (E) 
5  Of  9 


1  -»  I 


-> 


HBRG(I,J),HTE(I,J) 
J  =  1,NRDIPS 


(6F10.5) 


(90. -HBR5(I,l))/57. 29578 
->  HBRG(I,1) 


-  J 


2  -►  J 


HBRG(I,J-l)-(HBRG(I,J)/57.29578)  T 
-►  HBRG(I,J) 


c 


J  ^  NRDIPS 


YES 


^ 


li-l  -^I 


]c      ""    C 


I  ^  NRH 


D 


YES 


(i) 


NO 


Jfl  -»  J 
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\L 


HSPD*33. 76666 
-»  HSPDT 


:4^- 


C 


IITE  =  0 


m. 


YES 


V_ 


1    ~>  J 


^ 


hte(i,j)/hspdt 

-»  HTE(I,J) 


J+l->  J 


Ji^-'r    > 


<     ''"^      J   =   NRDIPS 


YES 


± 


I  £   NRHS 


"^     KQ 


YES 


PRINT      (D) 
6   of  9 


I+l-»  I 
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Q) 


C 


N-AME=DAT 


;^o 


GO    T0>< — ^^^^AME^GO^ 


^0 


Kl 


^^^  NAME=FINIS 


NO 


PRINT      (D) 
7   of   9 


YES 


DAT, NAME 
(F10.5.A6) 


NRU   ^   1 
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PRINT      (D) 
3   of   9 


HSPD,  TDD,  POiXMU ,  TRD, 

POK  VaJl ,  PR P3  ,  PC  TM  AX , 

NRDI PS , P , NRSW , NS , DL , T 

(I,PD(I),S/fT(I )  ,  I=l,NR3Vr) 


■^^^ — >(jJCAvTo 


>  (  S1»ITIME(  I )  ,  SOUS ( I )  ,  SSPD(  I ) , 
DEPTH(I),I   =   l.NRMAN) 


jS 
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PRINT      (D) 
9  of  9 


RETURN 
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Subroutine  OOMPRG*   This  subroutine  controls  the  compu- 
tation of  sonar  detection  ranges  for  each  helicopter.   The 
detection  ranges  are  computed  by  an  iterative  process  using 
the  applicable  zone  or  diffraction  region  equations.   These 
regions  are  determined  by  the  sonar  conditions  input  to  the 
program,  submarine  depth,  transducer  depth  and  the  helicopter 
figure  of  merit. 

COMPRG  is  a  modification  of  a  program  vn?itten  by 
R.  L.  Klinkner  of  the  Applied  Physics  Laboratory,  Johns  Hopkins 
University.   Subroutines  SGR(KP)  and  SSTCOP  together  with  sub- 
routine COMPRG  accomplish  the  sonar  detection  range  computa- 
tions.  The  first  two  of  these  subroutines  are  used  as  they 
were  written  by  Mr.  Klinkner. 
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Subroutine  COMPRG 


(General) 


l^NT^R 


® 


DO    FOR   EA'^H 
SONAR    SWEEP 


COMPUTE   COEPPICIENTS 
FOR   PROPAGATION 
LOSS    F.QUATIONS 


DO    FOR   EA:t 
HELICOPTER 


C 


ANY   LAYER 


) 


YES 


AVAILABLE    PROPAGATION 
LOSS    FALL    IN   ZONE   1. 


NO 


YES 


COMPUTE  RANGEN 
USING  ZONE  1 
EQUATION 


ARE  BOTH  THE 
PROJECTOR  AND 
TARGET  IN  THE 
LAYER 


YES 
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© 


OOMPRG    (G) 
2   of  3 


fAV  AIL  ABLE   PROPAGATION\__IES 
LOSS    PALL    IN   ZONE 


COMPUTE   RANGE^ 
USING   ZONE    II 
EQUATION 


NO 


COMPUTE   RANGED 
ZONE    III 
EQUATION 


,USING 


>/ 


COMPUTE   PROPAGATION   LOSS    FOR 
THE    ABOVE    RANGE    (IP    APPLICABLE) 
USING    PROPER   DIFFRACTION 
REGION    EQUATION   BASED   ON 
CRITICAL    PROPAGATION   LOSS 


C2 


IS    THIS    PROPAGATIO 
LOSS  LESS    THAN 
THE    AVAILABLE 

PROPAGATION   LOSS 


YES 


COMPUTE   RANGE   USING 
AVAIL /lBLE   PROPAGATION 

LOSS    AND   PROPER 
DIPPR.¥^TION   REGION 
EQUATION 


NO 


DETECTION   RANGE 

IS   LAST   RANGE 
COMPUTED    ABOVE 


<r 
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SUBROUTINE  TNE 
C 

C         •••  CIMENSION  AND  COMMON  STATEMENTS  ••• 
C 

8  IF(  NOVER)  9,g,lt» 

9  TIHET=Tir'E(2) 
NRUNT  =NRUN(2) 
NREVT  =NREV(2) 

C 

00  10  I   «  3,NTNE 
TIME(  I-n  *  T  IME(  I  ) 
NRUN(  I-l  )  «  NRUN(  I) 
K  =  I 

10  NREV(  I-l  )  «  NREV(  I  ) 

NINE  =  NINE  -  1 
C 

GO  TO  (  n,12,13)t  NREVT 
n  CALL  SME 

GO  TC  8 

12  CALL  HSM 
GO  TO  8 

13  CALL  l-DE 
GO  TO  8 

lU  RETURN 

END 
C 

SUBROUTINE  SNE 
C  ,   . 

C  •••  DIMENSION  AND  COMMON  STATEMENTS  ••• 

C 
C  SEARCH  EVENT  STORE  TABLE  FOR  PROPER  LOCATION  OF  THIS  EVENT* 

IT  =  NTNE 

NTNE  =NTNE  ♦  1 
1«4  IFCTI^'ET  -  TIKE(IT))  15,16,16 

15  TIME(  IT*1  )  »  TIME (  IT) 
NRUN(  IT* 1  )  =  NRUNI  IT) 
NREVI  IT4l  )  «  NREV (  IT) 
IT  »  IT  -  1 

GO  TO  lU  , 

16  TIME(  IT*1  )  «  TIMET 
NRUN(  IT  4l  )  «  NRJNT 
NREV(  IT  +1  )  =  NREVT 
RETURN 

END 
C 

SUBROUTINE  HSM 

c 

C  •••  DIMENSION  AND  COMMON  STATEMENTS  ••• 

C 

1  =  NRUNT 

IDIPN(  I  )  =  IDIPN(  I)  ^1 

200  J  =  ICIPNC I ) 
IF(NTE)2C1,201,260 

201  XH(  I  )=XHT(  I, J) 
YH(  I )=YHT(  I,J  ) 

TIMET  =  TI?'ET*HTE(I  ,  J)*TDD 

GO  TO  27C      ,   ■ 
C  ROUTINE  FCR  COMPUTING  BEARING  AND  TIMING  ERRORS  MAY  BE 
C    INSERTED  HERE. 

260  T1=HBRG( I,J)*ERB 

T2       «HTE(  I, J) ♦ERT 
C   COMPUTE  NEW  X  AND  Y  COORDINATES,  HTE  IS  JUMP  DISTANCE. 

XH{ I)  =  XH( I )  ♦  T2«C0SF(T1) 

YH( I )  =  YH( I)  ♦  T2»SINF(T1) 
C 

Tl  «  T2/HSPDT 
C   SET  UP  DETECTION  EVENT  FOR  THIS  HELICOPTER 

TIMET-TIMET  ♦TOO  ♦  Tl 
270  NREVT  =3 

NRUNT  =1 

CALL  SNE  / 

280  RETURN 

END 


141 


SUBROUTINE  SME 
C 

C  •••  CIMENSION  AND  COMMON  STATEMENTS  ••• 

C 

IF(KT-1  )  100, lOOt  108 

100  IF(NR/^NSS     )     101,101,102 

C   DETERMINE  SUBMARINE  COURSE  AND  CONVERT  DEGREES  TO  RADIANS 

101  RNC=3dO 
NXIN=PNC«RNG( 1 ) 
T1=FLCATF(NXIN)  /57. 29578 

C   DETERMINE  SUBMARINE  SPEED  AND  CONVERT  TO  YDS/MINUTE 

2  IF(NCAV)     3,3,M 

3  RNC  =  NHISPD-LOSPD*- 1 
NXIN=RNC«RNG( 1 ) 
T2=FLOATF(LOSPD>MXIN)»33. 766666    • 

C   DETERMINE  DEPTH 

RNC=MAXSC-MINSD*1 

NXIN^RNC«RNG( 1 ) 

DT=  MINSC+NXIN 

GO  TO  10 
C 

U  RNC=NCAVS(NCAV)-><CAVS(1  )♦! 

NXIN=RNC»RNG( 1 ) 

NSS=NCAVS(  n+NXn 

12=  FL0ATF(NSS)«33. 766666 
C  DETERMINE  MIN.  DEPTH  FROM  CAVITATION  DEPTH  VS.  SPEED  TABLE 

DO  6  I=1,NCAV  • 

IF(NSS-NCAVS( I ) )  5,5,6 

5  RNC=NCAVL(NCAV)-NCAVD( I)*l 
NXIN=RNC»RNG( 1 ) 

DT=  NCAVC(I)*  NXIN 
GO  TO  10 

6  CONTINUE 
C 

C   COMPUTE  VELOCITY  COMPONENTS  OF  SUBMARINE  AND  HELICOPTER 
C    DETECTICN  RANGES. 
10  VXS=T2«CCSF(T1 ) 

VYS=T2»SINF(T1) 

CALL  CCMPRG 

GO  TO  115 
C   IF  TRACK  IS  PREDETERMINED  ASK  IF  FIRST  COURSE  IS  RANDOM 

102  STIME  =  CO 

IF  (NRANC)  103,  103, 1 10 

103  RNC=360 
NXIN=RNC«RNG( 1) 
TT=NXIN 

GO  TO  lie 

108  IF(SCUS(KT)-SCUS( KT-1  ) 1110,109,110 

109  IF(SSPO(KT)-SSPD(  KT-1  ) )1 10,1 12,  110 

110  SPD  =  SSPD(KT)»33. 7666666 
CUS=(-SCUS(KT  )+TT 1/57.29578 
XS=XS*(SNTIME(KT)-STIME)«VXS 
YS'YS*( SMTIME(KT)-STIME)»VYS 

vxs=SPC»cnsF(cusi 

VYS=SPD»SINF(CUS) 

STIME  ^  SMTIME(KT ) 

IF(DEPTH(KT)-UT  1  112,12,112 
112  DT  =  CEPTH(KT) 

CALL  CCMPRG 
12  IFCKT-NRMAN) 1 13,1  lU,  1  lU 
C   IF  THIS  IS  NOT  THE  LAST  MANEUVER  STORE  MEW  EVENT  1  AT 
C    NEXT  MANEUVER  TIME. 
1  13  KT  =  KT4l 

TIMET»SMTIME(KT) 

NREVT  «  1 

SMNEXT  ■  TIMET 

CALL  SNE 

GO  TO  115 
IIU  TT^'O.O 

SMNEXT-  123U5.0 
115  RETURN  . 

END  .  ^ 
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SUBROUTINE  HOE 
C  .  • 

C   HELICOPTER  CETECTION  EVENT 
C 

C  •••  DIMENSION  AND  COMMON  STATEMENTS  ♦•• 

C 

I  ■  NRUNT 

IF(NSV^P(  I)-l)  300,300,302 
300  IOIP=iniP+l 
302  T1=TIMET 

T2  =  Tl-STIME 
C 

J1=NSKP(I) 

DO  32M  J=J1,NRSW 

XT  =  XS  ♦  VXS«  T2 

YT  =  YS  ♦  VYS»  T2 

RTT  =  SQRTF((  XH(  I  )  -  XT)»»2  ♦(YH(I)-YT  )••  2  I 

IF(RTT-DRNG(I, J))  310,310,319 

310  IF(RTT  -  SHORTR  )  319,319,311 

311  MOVER  =1 
IDIP=0 

NOET  =NCET  ♦  1 

Tl       =  Tl  ♦  SWT ( J)  -  OAT 

GO  TO  3UC 

319  T1=TUSWT(J) 
T2=T2*SWT(J) 

C  DID  SUBMARINE  MANEUVER  DURING  THIS  SWEEP. IF  ANSWER  IS  YES 
C   ASK  IF  THIS  IS  THE  LAST  SWEEP. 

IFCSMNCXT  -Tl)  320,320,3?»4     ,' 

320  IF(J-NRSh)    321,32U,32U 

C   RESTORE  THIS  EVENT  TO  ALLOW  SUBMARINE  TO  MANEUVER  BEFORE 
C    COMPLETING  SONAR  SEARCH. 

321  NSWP(  I  )=J+1 
TIMET=T1 
NREVT=3 
CALL  SNE 

GO  TO  3UC 
C 

32U  CONTINUE 
I       IF(  lOIPNC  I  )-NRDIPS)329.325.325 

325  IF(    IDIP-  12  )  3«40,356,356 

326  NOVER  =1 
IDIP=C 

GO  TO  3UC 

329  IF(FCPR)  33U,33U,330 

330  FCT=RNG( 1 ) 
IF(FCT-FCPR)  332,332, 33U 

C 

332  FCTx( FCT/FCPRI»F:TMAX 

Tl  «  Tl  ♦FCT 
33U  TIMET  =  TUTRD 
C 

NREVT  5=  2 
CALL  SNE 
NSWP(  I)=  1 
3U0  RETURN 
END 
C 

SUBROUTINE  PRINT  • 

C  . 

DIMENSION   NRCIO) 
C  •■•  CIMENSIOM  AND  COMMON  STATEMENTS  ••• 

C 

DO  50C0  1  «  1,10 
5000  NR(  I  I  «  1 
C 

10  READ  INPUT  TAPE  5,  1 000, NAME , N 1 , N2 , NS ,NU .NS 
1000  FORMAT(A6,5nO) 

19  IF(NAME  -  UHDATEI  25,20,25 

20  NDAY  «  Nl    . 
MONTH  »N2 
NYEAR  -N3 
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Jl  =  l 

J2  =  C 
NRU   =1 
KOET  '  0 
11  WRITE  OUTPUT  TAPE  6 , 3000, MONTH, NDAY.N YEAR 

3000  FORMAT!  1H1,3UX,36H  AHS-1  AN  ASW  HELICOPTER  SIMULATION  / 
1/  10X,5H  DATE  12,  IH/  12. IH/  12  ///I 

WRITE  OUTPUT  TAPE  6, 3001, NRU 

3001  FORMAT  (38X,26H  INPUT  DATA  FOR  RUN  NUMBER  tI3  ///) 
GO  TO  10 

25  IF(NAHE  -  5HNRSAM  )  27,26,27 

26  IF(KDET)  206,206,207 

206  NRSAM  =  Nl  .     , 

KDET  =  N2  I 

IRNC  =  N3  ', 

READ  INPUT  TAPE  5,1003,  (NSSIZE ( I ) , I« 1 tNRSAM  ) 
1003  FORMAT  ( 1216) 
GO  TO  10 

207  IF(J1-J2  )  209.208,208 

208  WRITE  OUTPUT  TAPE  6,20U0,J1 

20U0  FORMATC  //   9X,36H  COMPARATIVE  DETECTION  DATA  FOR  RUN 
1 ,6HNUMBER, 13///) 
Jl  =  NRU 
KDET  »  0 

GO  TO  225 

209  WRITE  OUTPUT  TAPE  6,20U1.J1,J2 

20U1  FORMAT(   //  9X,  36H  COMPARATIVE  DETECTION  DATA  FOR  RUNS 
1 ,I3,8H  THROUGH, 13  //) 
Jl  =  NRU  ' 

KDET  =  0 
GO  TO  225 

225    T2=JT 

BARD    =    BARD/T2 

VARD    =     (VARD    -    T2 •    BARD»« 2 ) / ( T2- 1 . ) 

WRITE    OUTPUT     TAPE     6 , 20U2 , BARD , VARD, MI  ND , MAXD 
20U2    FORMATl///    9X,32H    NUMBER    OF    DETECTIONS    PER    SAMPLE    // 

1  9X,ieH    SAMPLE    MEAN  =    F13.U    /9X,18H    SAMPLE    VARIANCE 

2  F13.U    /9X,  18H   MINIMUM  «    18/9X,     SH    MAXIMUM, 
39X,1H=     18     ) 

BARD=0.0 
VARD=C.O 
MIND*M321 
MAXO^C 
JT  =  0 
229     IF(K1  )    206, 206. Ul  1 

27  IF(NAME    -    6HNH1SPD     )     30,28,30 

28  NHISPC    «    Nl 

LOSPD  =^  N2  . 

MAXSD  =  N3 
MINSD  =  NU 
GO  TO  10 

30  IFCNAME  -  6HNRANS S )  UO, 3 1 ,U0 

31  NRANSS  «  Nl 
NRANC  =  N2 
GO  TO  10 

UO  IF(NANE  -  5HNRMAN  )  60,U1,60 
Ul  NRMAN  ■  Nl 

READ  INPUT  TAPES,  1001, (   SMT I  ME C I ),  DEPTH ( I »  .  ICUI  ( 11  . 


1SSPD( I ) ,  1=1. NRMAN  ) 
1001  FORMAT(UFl0.b ) 
GO  TO  10 


144 


60  IF(NAKE  -  UHNRSW  )  65,61,65 

61  NRSW  =  Nl  '  . 
READ  INPUT  TAPE  5,  1002, ( SWTt I ) ,P0( 11 ,  I  «  1,NRSW  I 
GO  TO  10 

C 

65  IF(NAHE  -  3HTD0)  68.66,68 

66  READ  INPUT  TAPE  5 , 1 006,TOD, TRO, FCPR,FCTMAX 
1006  F0RMAT(UF10.5) 

C 

GO  TO  10 
C 

68  IF(NAKE  -  MHNRHSI  70,69,70 

69  NRHS  =  Nl 
NRDIPS  »N2 
HSPD=N3 
NTE=NM 

C 

DO  5  1=1, NRHS 

52  READ  INPUT  TAPES,  1002,  (  HBRGd  ,  J  )  ,HT  E  (I ,  J  )  ,  J«l  ,NROI  PS  ) 
1002  FORMAT  (6F10.5) 

C 

CONVERT  HELC  DISTANCES  TO  TIME  ENROUTE  AND  DEGREES  TO  RADIANS 

HBRG(I,1)«  (-HBR5(  I,  l)+90,0)/57.29578 

DO  5  J=2, NRDIPS 
5  HBRG(  I,J)»(-HBRG(  I,J)/57.29578)*HBRG«  ItJ-ll 

HSPDT  =  HSPD«  33.76666666 

IF(NTE)  53,53.55 

53  00  5U  I»1,NRHS 
DO  5U  J=l, NRDIPS 

5U  HTE(I,J)  «  HTE( I,  JI/HSPDT 
55  GO  TO  10 
C 

70  IF(NAKE  -  MHXDAT  )  80,71,80 

71  XDAT=N1 
YDAT  «  N2 
GO  TO  10 

C 

80  IF  (NAME  -5HS0NAR  )  90,81,90 

81  NS=N1 

READ  INPUT  TAPE  5, 1005,FOMMU,FOMVAR,OL,T,F,TS,PCL 
1005  FORMAT(  7F10.5) 
GO  TO  10 

c 

90  IF(NAME  -  3HDAT)  100,91,100 

91  READ  INPUT  TAPE  5,  1 008, DAT , NAME 
1008  FORMAT  (   F10.5,A6  ) 

IF  (NRU   -  1)  93,93,92 

92  WRITE    OUTPUT    TAPE     6,     2001, NRU 

2001  FORMAT    ( 1H1,////U5X,26H    INPUT    DATA    FOR    RUN    NUMBER    13    ///) 

93  WRITE    OUTPUT    TAPE     6 , 2002 , NAME , NRHS, DA T 

2002  FORMAT(26X, 12HSEARCH    PLAN    ,A6,6H    WITH     12,    9H    HELICOPT 
1.8HERS    AND    ,F5.0, 18H    MINUTES    TIME    LATE    //    ) 

GO    TO    10 


C 
C 


t 


c 
c 


100    IF(NAME-UHNCAVI    110,102tll0 
102    NCAV»M 

READ    INPUT    TAPE    5,999,     ( NCAVSI I) , NCAVOf I) , 1-1 ,NCAV 1 
999    FORMAT( 1216) 

NHISPD  =  ncavs(n:av) 

LOSPO  =NCAVS( 1 ) 
MAXSD  '    NCAVDINCAV) 
MINSD  =  NCAVD(  1  ) 
GO  TO  10 

110    IF(NAME-2HG0)    311,330,311 

311  IF     (NAME    -    5HFINIS)    312.313,312 

312  WRITE    OUTPUT    TAPE     6, 2000, NAME 

2000    FORMAT    ( 1 H 1 , 3 1 X ,3 6HF0LL0 W I NG    DATA    CARD    CONTAINS    ILLEGAL 
1.5H    NAME    //12X,A6,U2X, 12HJ0B    DELETED.     ) 
§TOP 


l^O 


313  Kl=l 

GO  TO  207 
C  PRINT  OUT  HEADING  ANO  INPUTS 
C 

330    WRITE    OUTPUT    TAPE     6,2010,HSPD, TDD , FOMMU, TRD.FOMVAR, 
lNRniPS,F,NRSW,NS 

2010  FORMAT    (/    25X , 1 5HHEL I  COPTER    DATA.35X,  llH    SONAR    DATA    //    . 
19X,nH    JLMP    SPEED,  17X,3H    =    F5.  0,  5HKN3  TS,25X  , 

217H    FIGURE    OF    ME^IT    // 1  OX, 25HT IME    TO    LOWER    TRANSDUCER 
3.3H    «    F5.1,8H    MnuTES,2UX,12HAVERAGE  =    F8.3    /ICfX, 

U29HTIME    TO    RAISE    TRANSDUCER  *    F5.1,8H    MINUTES, 2MX, 

512HVARIANCE         =    F8.3    // 1  OX, 2UHNUMBER    DIPS    EACH    HELICOP 
6,5HTER    *    IU,31X,23H    ACOUSTIC    FREQUENCY         -    F6.2,6H       KG. 
7//10X,29HNUMBER   SWEEPS    EACH    DIP  =    IU,31X,5H    SEA    » 

816HSTATE  =     IM/     ) 

WRITE    OUTPUT     TAPE     6 , 20 1 2.01 , T , ( I , PD ( I ) , SWT ( I ) , I «1 , NRSW I 
2012    FORMAT( 1CX,5HS0NAR,7X,5HDEPTH,6X,6H   SWEEP. 35X, 

122H    LAYER    DEPTH  =    F6.2,6H       FEET    /! OX, 5HSWEEP , 

220X,     HTIMES.35X,23H    TEMPERATURE    IN    LAYER    -    F6.1,UH      F. 
3//( )0X, I2,F1U.0,F 12.1     )//     ) 
C 

WRITE    OUTPUT    TAPE     6, 201 1 , FCPR , FCTMAX 

2011  F0RMAT(/9X, 31H    PROBABILITY    OF    FALSE   CONTACT    -    ,F6.5    /9X, 
131H    MAX.    DELAY    FDR    FALSE    CONTACT    «    F5.1,8H    MINUTES    /    I 

380     IF(NRANSS)382,382,M01 

382  IF(NCAV)  UOO, UOO,  38U 
C 

38U  WRITE  OUTPUT  TAPE  6,2018 
2018  F0RMAT(//U7X,  25H  NON-CAV I  TAT  I  NO  SUBMARINE) 
C 

«400  WRITE  OUTPUT  TAPE  6,  2020  ,  LOSPD  ,NH  I  SPO  ,  M  INSD  .M  AXSD  ' 

2020  FORMAT(/25X,UOHSJBMARINE  COURSE,  SPEED  AND  DEPTH  CHOSEN 
1,  9H  RANDOMLY// lOX, 16HMINIMUM  SPEED   =  lU  /lOX, 
216HMAXIMUM  SPEED   =  lU  // 1 /X , 1 6HM INIM UM  DEPTH   «  15  / 
310X, 16HMAXIMUM  DEPTH   «  15  I 

GO  TO  M05 
C 

UOl  IF(NRANC)  U02,»402,U03 
U02  WRITE  OUTPUT  TAPE  6,2021 

2021  FORMAT(  3UX,  37HSUBMARINE  TRACK  PRESET(FIRST  COURSE 
l,15HChOSEN  RANDOMLY   //  ) 

GO  TO  UOU 

c 

U03  WRITE  OUTPUT  TAPE  6,2022 

2022  FORMAT!  M9X , 22HSJ BMAR INE  TRACK  PRESET  //  ) 

UOU  WRITE  CUTPUT  TAPE  6,2023,  (  SMT IME( I) , SCUS ( I ) ,SSPD( I ) t 
iDEPTHd  )  ,  I«1,NRMAN  ) 

2023  FORMAT( ICX.UHTIME  ,6X,  7H  COURSE, 6X,1 UH  SPEED  CKNOTSIf 
16X,13H  DEPTH  (FEET)  / / ( 1  OX , F6. 1 ,F 1 1 .0 ,F 1 6.0,F22.0  )) 

C 

U05  WRITE  OUTPUT  TAPE  6,2030,NRU 

2030  FORMATC     ///USX , 26H0UTPUT  DATA  FOR  RUN  NUMBER  I3  //) 
WRITE  CUTPUT  TAPE  6, 203 1 , (NR ( I ) . I = 1 ,N RHS) 

2031  FORMAT!///   30X,10H  SUBMAR I NE , 1 OX ,  1  7H   SEARCH  TIME  TO 
1 , lOHDETECTION  , 7X , 20HNUMRER  OF  DETECTIONS  /30X, 
210HDETECTIONS, lOX ,25H   WHEN  SUBMARINE  DETECTED  lOX, 
320H  BY  EACH  HELICOPTER   /59X,9H( MINUT ES )  /lOX, 
U15HSAMPLE    S  AMPL  E  ,  2  1  X  ,  5H  ME  AN,  Ux  ,  17  HM  INIMUM    MAXIMUM 
5/10X,i4UHNUMBER  .    SIZE     NUMBER    PERCENT    TIME   t 
628HVARIANCE     TIME        TIME   tlOU  T 


C 


WRITE  OUTPUT  TAPE  6,2032 
2032  FORMAT!//) 

NRU   -  NRU   ♦  1 
J2  «  J2  ♦  1 

UIO  RETURN 

un  STOP 

END 
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SUBROUTINE  COMPRG 
C 
C  •••  CIMENSIOM  AND  COMMON  STATEMENTS  ••• 


C 


00  22C  J=1,NRSW  . 
DP  =  PD(J) 
CALL  SETCOF 

DO    220    I    =     1,NRHS 

PLG=PL(I) 

IF(DL)    2C>»,20U,202 

202  IF(PLG    -    P(l))    203,203,210 

203  CALL       SCR( 1 ) 

DRNG( I, J)    '    RG«1000.0 

IF(OT    -    CD    22.22  ,20U 
22    IF(DP    -    CD    220,20«*,20U 
20U     IF    (RC)    27,27,2U 
2U    IF(PLG-    FL»    LOGF(RG)    -    B(5)»RG    -    C(5))    220,220,206 

206  CALL    SCR(5) 

DRNG( I.J )    =    RG»1000.0 
GO    TO    220 
27    IF(DD    208,208,207 

207  IF(PLG-    A(i4)»FLN»LOGF(RG)-B(«4)»RG    -  C  (H)  )    220,220,208 

208  CALL    SCR(U) 

DRNG( I. J )    »    RG-1000.0 

GO    TO    220  ^  ■  ' 

210  IF(PLG-    P(2n    ^11  ,211,212 

211  CALL    SCR(2) 

DRNG( I.J)     »    RG»1000.0 
GO    TO    20U 

212  CALL    SCR(3) 

DRNG( I.J)    -    RG'1000.0 
GO    TO    ^0*1 
220    CONTINUE 
RETURN 
END 

SUBROUTINE    SCRIKP ) 

•••    CIMENSION    AND    COMMON    STATEMENTS    ••• 

IF     (KP    -    3)  110,11  1,110 

110  C2    «    6.0206 
GO    TO     10 

111  C2    '    3.0103 

10  CO    «    e(KP) 
CI    »    C(KP) 

IF     (PLG-CO-Cll       lU, 11,11 

11  CU=    CO 

C5«    C2    ♦    2.«CU    ♦CI 
C6«    0. 
C7»    1. 

12  IF    (C5-PLG)       13,17,17 

13  CU    »    2.«CU  • 
C5    »    C2    4    CU    ♦    C5 

C6    =    C6    4    C7 
GO    TO    12 
lU    CU    =    CO*. 5 

C5    -    CI    4    CM    -    C2 
C6    «    -1. 
C7    =-1. 

15  IF     (PLG-C5)       16,17,17 

16  CU    =    CU».5 

C5    «    C5    -    C2    -    CU 
C6    «    C6    ♦    C7 
GO    TO     15 

17  TOP    »    PLG    -    CI    -   C2«(C6-. 95693) 
DENOM    «    CO    ♦    C2»(  .5««C6) 

C12    »    TOP/DENOM 
C15    -     I  1.UU269X2 

18  Tl    -      C15»    LOGF       (C12)    ♦    C0«C12    ♦    01    -    PLG 
T2    «    IC15/C12)    ♦   CO 


1  ■'♦7 


C13    =    C12-(T1/T2) 
IF    (C13)       25,25.19 

19  IF    {AESF{C13-C12I-.001)  21ff21«20 
21    RG    =    C13 

GO    TO    50  ; 

20  C12    «    C13 
GO    TO    18 

25  WRITE    OUTPUT    TAPE     6,26 

26  FORMAT( IX.UOH    PR3GRAM    DETECTED    AN    ERROR    -    NEGATIVE    ORt 
118HZER0    RANGES    IH    SCR) 

CALL  EXIT 
50  RETURN 
END 

SUBROUTINE  SETCQF  ,. 

•••  CIMENSIO^  AND  COMMON  STATEMENTS  ••• 

111  IBIT=2 

A6=(F/25.)» •0.333 3333333 

A7=(  1  .23)«(  10.«»(  6.-(2100.)/(T  +  i*59.6»  )) 

A8=F»F 

A10=SCRTF( ABSF(DP-DL) )  ♦  SQRTF ( ABSF (D T-DL) ) 

A16=(0.651»A7»A8|  / ( A8*A7» A7) ♦ I 0.0269* AS  I /A7 

IF  (DL)  133,133,112 

1 12  AO=SQRTF(DL) 
A1=SQRTF(DP) 
A2=A1/A0 
A3xSQRTF(DT) 
AU=A3/A0 
A5=(U.5/A0»«SQRTF (F) 

IF  (F-8. )     Hi,  iiu,  nu 

113  A9=l. 

GO  TO  ns 

nu  A9=1.M62«A6 

115  IF  M  .-A2)   1  16,1  19, 119 

116  IF  (1  .-AM)   1  17.1  18,1  18 

117  Al 1'0.2»(SQRTF(A2»A2-1.)  ♦  SQRTFI AU^AU- 1 . ) ) 
GO  TO  122 

118  All:=(  1.-AU)».25  ♦  .2»ISQRTF{A2»A2-1.)  ) 
GO  TO  122 

119  IF  ( l.-AM)    121,  120,  120 

120  Al 1«( 2.-A2-AU)«0.25 

IF  ( A2-AU)  122,200, 122 
200   A12  «  FL»(-6.9) 
GO  TO  12 

121  All     =  (  l.-A2)»0.25>0.2»(SQRTF(AU»Al|-l.)  ) 

122  A12=»  FL«LOGF(  AO»M  1  ) 

12    A13»    FL»LOGF(  (AO)  •(A1U.5)) 
Al 8=A8SF ( A2-AU ) 

123  A15=     .U»A9»( {  10.»»A2)>( 10.»»AU)4(10.4«A18)I 
IF     (A18-1 . )        12U,  125, 125 

12U    MU    =    .  1«A6«(  (  10.  )«»(2.3»A18)  ) 
GO    TO    126 

125  AlU    =    20. "Ad 

126  IF    (NS-3)       127,128,128 

127  A17    =    A5 
GO    TO    129 

128  A17    =    2.«A5 

129  A19  =  2.«(A15-A1U ) 
A20  «  A1UA19 

A27  =  A16  ♦  (A19/AC) 
A28  «  Al*«  ♦  60. 

A26    -    A2e    t    A13   t     ASTlAoIdiw.sf    «   A20 
IF       (CD       133,133,130 

130  A( 1 )»20. 

B(  1  )-    A16    ♦    (Ali»/(AO»An  )  ) 
C(  1  )»    60. 
P(1)«    A25 


i4ij 


i 


131  A(2)=  20. 
e(2)=AU4(A19/A0) 
C(2)=  A2e  -  A20 
P(2)=  A26 

132  A(3)=  10. 

B( 3)=  A16  ♦  A17 

C(3)»  A15  -  A17«A0«(A1U,5)  ♦  .5«A13  4  60, 

P(3)=  A26 

133  A(U)=  20. 

B(U)=  A16  ♦  5.»A6 

C(U)«  (25.-Al0)»A6  ♦  60. 

RC  =  (AlC-25. )«.2 

IF  (RC)   33,33,34 
33  PLC=123U56. 

P(U)=PLC 

GO    TO    13M 
3U    PLC=    FL«LOGF(RC)    ♦A16«RC    ♦    60, 

P(U)     =    PLC 
13U    A(5)     =    2C. 

B(5)    »    A16 

C(5)     =«    6C. 
P(5)    =    PLC 

RZl    '       AC'All 

RZ2    =       AC»( Al  14.5 ) 

IF    (RZl-.Ol)       135,135,136 

135  P( 1  )«-2^68. 

136  RETURN 
END 

SUBROUTINE  NORM 

•  ••    CIMENSIO>J     AND    COMMON    STATEMENTS    ••• 


YNORM  =  CO 

DO  10  IN  =1,12 
10  YNORMxYNCRM^RNGd  ) 

YNORM  =.  XMU*(  YNORM  -  6.0  )  •  SQRTF  (  SIG2) 

RETURN 

END 

FUNCTION  RNG(NR) 
C 

C  UNIFORM  (0,1)  FLOATI><G  POINT  NUMBER  GENERATOR.  (NR.O  -RESET  • 
CNR'l  -  GENERATE  NEXT  NUMBER  IN  THE  SEQUENCE  AND  PLACE  IN  A 
C  REGISTER. NR=2  OR  MORE-  CYCLE  GENERATOR  NR  TIMES! 
C  NUMBERS  GENERATED  BY  POWER  RESIDUAL  METHOD. 


C 


CON    (    Kl=12207031 25       ,    K2«2000000000000000B,       R«77777777Bt 
1    Rl=77777777    B)  i      , 

IF(NR)  10,10,20 
10  R=R1 

GO  TO  UO 
20  on  30  I»1,NR 

LDA(R),MUF(Kl ),0VF(K2),STQ(R) 
30  CONTINUE 

LLS(i48),CVF(K2),^RS(  n),ADD(K2),FAO(K2) 
kO    RETURN 

END 

END 


SAMPLE  PROGRAM  OUTPUT 

coo: 


• 

o     ►-    • 

. 

^       LULL 

LU 

U. 

o        o 

o           • 

•     oo 

O      Ou-j 

rj<\j   • 

c 

oo 

II          O  !■ 

oo 

»— 

Ui 

CJO 

II    II  OL 

►- 

•    • 

Oi 

< 

o 

>•        > 

.1 

♦-     o 

LJ           < 

«1 

mm        ^— 

Z            -J 

UJ 

►- 

oc 

LU 

at 

«i 

UJ  II    II 

_J          z 

»4 

(-) 

X 

CJ            ^ 

►- 

LJ 

a 

u. 

OC        XLU 

t/> 

« 

U       U) 

LL       »-Qt 

u 

^ 

UJO 

LUa_3 

^ 

■    u 

\JLJ-Ji^ 

O^LU  — 

O 

t>0 

ftr<j< 

—  <CJ< 

ij 

_ja  — 

►->-     a 

KM 

OLD  or 

oni/ict:iu 

a 

2. 
• 

—  >< 

3     Lua. 
002. 
OaJ«luj 

OL 

^  1 

<C/)_J>— 

JJ 

-- 

cu 

2. 

Z) 

A. 

^ 

<• 

^ 

K/~: 

3 

a. 

or 

L.- 

oc 

a 

o 

u 

u. 

< 

_i 

^- 

LU 

« 

X 

o 

-J 

v>0i/1 

h- 

OJOJ 

:d 

(•)►—►— 

CL 

X 

•-ZJ3 

z 

K- 

oz?- 

a 
a 

< 
a. 

z 

UJ 

u 

X 

o 


o>- 

-iQ. 

»-(_) 

^ 

OO 

J-f- 

UJ"— 1 

h-_» 

LUUJ 

CM 

OX 

ror- 

*     U-i 

uo 

^- 

< 

Osl'- 

QiUJ 

UJ 

aj> 

o* 

2.aj 

r' 

-^ 

Z 

ac    fo 

2.2.     • 

XH-CM 

< 

z 

2. 

OQ 

•* 

•—LU 

►-I— 

OO 

a:     o 

K 

LULJ 

_>    u^ 

UI 

^-^— 

2.-J    » 

OJ 

LUUJ 

—  2.0J 

a. 

uu— 

•Z  — 

J> 

uo— .►- 

z 

OUJJJ2. 

•-z»- 

z 

— •_) 

^ 

UJQCZ 

a 

>:<i- 

I       UJO» 

^2.3: 

0'*^ 

«v 

H-X  — 

■      Z(N 

l3 

ID 

<^ 

LL 

Iv/1 

»—    • 

O 

Qcr- 

< 

az 

<lCNi 

►- 

<LU 

> 

^ 

OJX 

o 

UO^ 

»- 

Zluoj 

^ 

<i3:f^ 

& 

LU— •    • 

^ 

a:t-o 

z 


UJ 
Q. 

</■" 

LU 

t-  UJ 

13  4/) 

Z  K 

—  ^ 

«              z—  ac  o                                 ;3 

X               ^^.j:  o  o                                   O 

•  oo 

oiDin  o  •  UJ 

fO  — >—  f<rirj  •— 

^  •  oc 

«         «  II  II  < 

►-  II    II    II         II    II  O-U,      O  <"  ii. 

«i  UJ2:      •  •    >—h-     CJ 

Z        a  cc     Lui—    — —     '^<^     3 

<  _J  i»-  <3<        LO 

-J  OC  I-  00  l-h- 

(LUJCX.  (XLL  ZZ  ^ 

^  UJLU        O  OO 

X  Q.  OO  oa.  oo 

o  o  _j_j  — —  oo 

tf  o  oo  _jo  UJLU  irnr>     oo 

<  •—  l/lOO  LU  oolo  ^      '-vn  Z 
UJ  -i  zz  XX  -JO 

\A           UJ           <<            OX  •    •       <X< 

I            ^ty  X«l>—  OO      XU.  MM        M    II 

>— k—  OUJO.  iTIlP  UJ 

<        ^  r-        u.a  ^LUO 

CKLU  UJi/)0  OO  OO       XX                                                                      Q.^s|0 

uJ4/1              i.  X  UJLU       h->—  2.  ~^•— 

:X—'  UOUJ  >-  UJLU        0.0.                                                                                <lOO 

aLJ<  Q-uJ  I— >•  O.Q.       LULU                                                                          (/) 

uj_ja     •-':x  -^<  i/)oo    OQ 

UJ  OOO  -J_) 

a.oo  --LU  acs:     acs: 

ooKi-     oca  cto  ziji     uu  Luoc 

UJUJQ^a.  <  2.2.       2.2.  OLU 

a..ijxc_iij'<uj  'J3»  ►_-^i--.~.  iajr— 

aL.2.>.       2.J-ZUJ      •— Oi      ox  ZXZX  2.2. 

:_).-.—     _3_)Oi  a.«5  — <i     — '<i  O 

-)►-•—      ZZoOoO  Q.2.  2.2.       2.2.  oOZ 


zo 

UJO 

oo 

LUOO 

oc   • 

zz 

UU'O 

«o 

Cl-» 

oc—' 

<»- 

2.0 

ojlu 

oc 

u^- 

LU>0 

OOLU 

OJJ- 

LJ 

2. 

3 

£ 


© 


lOMPRG  (G) 
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RETURN 
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Subroutine  NORM.   The  function  of  Subroutine  NORM  is  to 
compute  the  value  of  a  Normal  random  variable,  employing  the 
classical  central  limit  theorem.  Given  a  sequence  of  inde- 
pendent, identically  distributed  random  variables  X.   with 
finite  mean  E(X)  and  standard  deviation  (XCx) ,   then  the  se- 
quence Y^  defined  by 


n 


_  (X^  +  X2  + 


•  9  0 


+  Xn)  -  nE(X) 


V^  cru) 


(6.9) 


converges  in  distribution  to  a  random  variable  which  is  norm- 
ally distributed  with  mean  zero  emd  variemce  one.  I27J 

In  actual  practice  the  value  of  a  pseudo-normally  dis- 
tributed random  variable  is  computed  since  pseudo-random 
numbers  generated  by  RNG  are  used  for  the  sequence  Xj.   Under 
the  assumption  that  the  X.   are  uniformly  distributed  in  the 
interval  [o,l]  then  E(X)  =  h   and 

crw   =/i= 


Talcing  n  as  twelve  we  compute 


Y  = 


VT7 


12 


VT2 


[(X, 


+  X2  +  ••••»•  X-j^2^ 


(6.10) 


which  is  approximately  normal  with  mean  zero  ajid  variance  one. 
To  obtain  the  value  of  a  normal  random  variable   Z,  with  mean// 
and  variance  (T         we  use  the  transformation 


Y  = 


-M 


or 


cr 


(6.11) 
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The  follovlns  calling  sequence  Is  used  for  Subroutine  NORM: 

/^-»  XMU 

O'^-^  SIG2 
CALL  NORM 
The  resulting  normal  varlate  will  be  placed  Into  the  var- 
iable YNORM  by  the  subroutine. 
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Subroutine   NORM 


(General) 


ENTER 


GENERATE    .0 

SEQUENCE   OP 
.2    RANDOM 
NUMBERS 


COMPUTE  Y'^N(0,1) 
Y'=  ^    (RN)^-6 
1=1 


COMPUTE  Y'^N(;i,tf^) 
^-wrY'    -»  Y 


RETURN 
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Subroutine  NORM 


(Detailed) 


ENTER 


0.  -^  YNORM 


\/ 


1  ->  I 


-^ 


GALL 
RAND 


Y+RN  ->  Y 


C 


I  ^  12 


> 


NO 


> 

YES 

X14U+(Y-6.)^^  l/siG2' 
->»•  YN0H14 

N 

/ 

RETURN 

I+l  -»  I 
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APPENDIX  A 
LISTING  OF  AHS-1  FORTRAN  COMPUTER  CODE 

PROGRAM    AHS    1 
C 

DIMENSION       NSSIZE  (  1  CO ) .NHSDET ( 1 0) .HTE ( 10f20),HBRG(  10f20) 
1 ,XH( 1C).YH( 10>,XHT( 10,20) .YHT( 10.20). IDIPN( 10),SMTIME(20) 
2,DEPTH(20),SCUS{2  0)  ,  SSPD  (  20  )  ,  PD  (  5)  ,  TI  ME  (  31  )  ,NREV(31  )  ,. 
3NRUN( 31 ) ,DRNG{ 10,  5 ) , S WT( 5  )  ,N SWP ( 1 0  I ,NCAVD ( 30 ) ,NCAVS(30) 
U.A(10),B(10),C(1D),P(10),PL(10) 

COMf'CN    IT,  JT,  Tl,r  2,  IDIP,  I2,TT 

COMMON    NRSAM,NSSI ZE, NHSDET, TMIN.TMAX,  RART,VART,MINOt      r 
IMAXO, eARC,VARO,PERC,DAT,XDAT,YDAT,NOET,NOTEMP,PL 

COMMON    TIMET, NREVT,NRUNT,NTNE, TIME, NREV,NRUN 

COMMON    KT,NRANSS,NRANC,STIME,DT,VXS,VYS,YS,XS,SMTIMEf 
ISCUS, SSPC,NRMAN,^HISPD,LOSPD,MAXSD,MI  NSD ,NCA V , NC AVS t 
2NCAV0, DEPTH 

COMMON    ICIPN,NRDI PS , HTE , HBRG , XH , YH , NR HS , HSPO , HSPDT , TOO 
1  ,XHT,YHT,NTE 

COMMON    DRNG,NRSW,  S WT , PD, TRD ,NOVER , SHORTR , SMNEXT, NSHP 

COMMON    FCPR,FCTMAX,FCT  .     .» 

COMMCN    NS,T,F,DL,  OP,TS,FOMMU,FOMVAR 
'  COMMON    A,B.C,P,PLG,RG. AR,FL,FLN,FL3 

COMMON  IBIT, IFS,IR,PLC,RC,RS,RT,RZ1,RZ2.SG 

COMMON  AO.aI,  A2,A3,AU.A5,  A6,A7,A8.A^,  AlO,An.Al2,A13, 
1A17,A18. Al9,A20,A2l,A22,A23,A2U,A25,A26,A27,A28,A29,A30 
2,A1U, A15,A16 

COMMON  Ci,Cl,C2,CU,C5,C6,C7,C12.C13,C 15,DEN0M,T0P 

COMMON  YN0RM,SIG2 ,XMU 

XDAT=0.0 

YDAT=C.O  . 

IRNC=C 

NCAV  =0 

FCT=0.0 

FCPR'C.O 

JT»0 

BARD  =  0.0 

VAPD  *  0.0  .     '     ' 

MIND  =  ICOO 

MAXD  =  0 

FLN  =.U3»429U»*819 

FL  =  20. •FLN 

FL3    =     10. •FLN 
1     STIME    =    CO 

NRANC  =  C 

NRANSS  =  0 

CALL  PRINT 

SMNEXT  =  123U5.0 

TT=90. 
C 

C   TEST  NTE  TO  DETERMINE  IF  HELD  BEARING  AND  TIMING  ERRORS 
C   ARE  TO  BE  COMPUTED.  IF  NOT, COMPUTE  DIP  STATIONS. 

IFCNTE)  11,11,12 

11  DO  12  I=1,NRHS 
T3=XDAT 
TU'YDAT 

00  12   J=1,NR0IPS 
Tl  »  HBRG(  I,J  ) 
T2  »  HTE( I, J) 
VXH  »  HSPDT  •COSF(Tl) 
VYH  «  HSPDT  -SINF (Tl  ) 
T3  »  T3+VXH»T2 
TU  »  T»4  +  VYH«T2 
XHT( I ,J)=T3 
YHT(I,J)=TU 

12  CONTINUE 

NOTE.  DIMENSION  AND  COMMON  STATEMENTS  ARE  IDENTICAL  FOR  ALL 
SUBROUTINES  AND  ARE  LISTED  ONLY  IN  THE  MAIN  PROORAM 
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JT=JT*NRSAM  ' 

XMU  =  FONMU 

SIG2=F0MVAR 

TIME( 1)  =  0.0 

12=  NRHS«NRDIPS 

IDIP=0 
C   CYCLE  RANDOM  NUMBER  GENERATOR 

T3=RNG(0) 

T3=RNG( IRNC) 
C 

DO  30  N=1,NRSAM 

DO  6  I  =1,NRHS 

6  NHSOEK  11=0 
BART    =0.C 
VART    »0.C 
TMIN    =    6C0.0 
TMAX=0.0 
NOET=C 

N5  =  NSSIZE(N» 
C 

DO  20  P=1,N5 

DT=-123U5. 
C   SET  UP  FIRST  SUBMARINE  MANEUVER 

NTNE=  1 

KT=1 

XS  =  XDAT 

YS  =  YDAT 

TIMET  =  CO 

NREVT  =1 

CALL  SNE 
C  SET  UP  FIRST  DIP  FOR  €ACH  HEL ICOPTER, INIT I AL IZE  DIP  COUNTER 

DO  7  1=1 ,NRHS 

IF ( NT E)  16, 16,  15 

15  XH( I )=XDAT 
YH( I )=YDAT 

16  IDIPN(n=0 
NSWP(  n«  1 
TIMET  =  CAT 
NRUNT  =  I 

NREVT  =2  • 

CALL  SNE 

CALL  NORM 

PL( I)=  (YNORM  ♦  TS  )/2.0 

7  CONTINUE 
NDTEMP  sNDET 
NOVER  =  C 

C   TAKE  FIRST  EVENT. TNE  WILL  HAVE  CONTROL  UNTIL  DETECTION  OR 
C   LAST  DIP 

CALL  TNE 

IF(NDTEMP-NDEr )  10,20,20 
10  I  =  NRUNT 

NHSDET( I )=NHSDET( I )  ♦  1 

BART  =  BART  ♦  Tl 

TMIN  =  MINIFCTMH  ,T1  ) 

TMAX  =  MAX1F( TMAX .Tl) 

VART  =  VART  ♦  Tl  ••2 
20  CONTINUE 

M=NSSIZE(N) 

T2  =  NDET 

BART  =  BART/  T2 

VART  =  (VART  -  T2«BART»»2)  /  (T2-1.0) 

PERC=(T2/FL0ATF  (NSSIZE(N)))«100.0 

WRITE  OUTPUT  TAPE  6 , 1 00, N, M, NDET , PERC , BART , VART , TMIN , TMAX t 
1  (NHSCET(I),  I=1,NRHS) 
100  FORMAT  {  IIX.  I3,7X,  I3,6X,  I3,3X,F7.3fF8.2fFn.U,2Fl0.2t2X,(9IU) 

MAxo  «  xmax6f(maxd,noetj 

MIND    ■    XM  IN0F(MIND,NDET) 
BARD    =    BARD    ♦    FLDATFCNDET) 
VARO    «    VARO    ♦    FLOATFINDET    ••2) 
30    CONTINUE 
GO    TO    I 
END 
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